I would like to solve ODE and get a function. not just a final value. so i'd like to know a value of y(t) at all steps, for example incremented by 0.01. from my understanding the dodesol returns y(tmax) only. is there an efficient way to get the intermediate values? the routine that i was using before from numerical recipes used function evaluation during solving ode with adaptive step-size to interpolate the function so i can get an output for each increment of 0.01.
For more complete information about compiler optimizations, see our Optimization Notice.