# Why do NDSolve and OutputResponse not evaluate non-analytic functions numerically?

test = OutputResponse[TransferFunctionModel[1/(1 + s), s], Exp[-(1/t)], {t, 0, 10}]


does not evaluate. If Exp[-(1/t)] is replaced by Exp[-t], it's fine.

Starting the numerical evaluation at t>0 does not help, either. (v9.0.1.0 on Mac OS). Thanks.

The same issue occurs with NDSolve:

NDSolve[{y'[t] + y[t] == Exp[-1/t], y[0] == 0}, y, {t, 0, 10}]


does not generate an interpolating function. Replacing Exp[-(1/t)] by Exp[-t] does.

-
This behavior seems rather strange, and some suspicious messages are produced. Starting from t == 1 doesn't help; changing Exp[-1/t] into Exp[-1/(t + 1*^-10)] does. Bug? –  Oleksandr R. May 4 at 19:53
Thanks! The t to t + eps trick does seem to work for both OutputResponse and NDSolve. –  John Bechhoefer May 4 at 23:01
It probably should be noted that it does help to start from a value larger than zero for NDSolve (at least if one adopts the initial condition), but not for OutputResponse... –  Albert Retey May 6 at 15:44