MATLAB: Numerical solution of linear ODE involving function defined point-wise

Question

I need to solve a differential equation of the form y'(t)=1/z(t)-y^2 z(t) where z(t) does not have an analytical form, but is defined point-wise as a list of values. Is there any integrator in Matlab that can handle this direcly?

Answer 1

See Example 3 (link) in the ode45 (link) documentation.