I'm working on a project that involves a solid particle passing a supersonic shock. I'm trying to find the temporal response of the temperature numerically.
It all comes down to a final DE of the form: T'(t)=f(t)*T(t)
f(t) is a function that I retrieve from interpolating between several functions fitted through measurement points.
Evaluating the code returns the following error: NDSolve::ndnum: Encountered non-numerical value for a derivative at t == 0.`. >>
When I replace f(t) with a more simple function it works without issues. (Hint: the actual f(t) function is defined and non-zero for every t!)
Can someone spot what I'm missing here?
Attached you will find the complete notbeook. The interesting part is the last two cells: first one will result in an error, second one will work.
Attachment: forum.nb, URL: ,