How to write an ODE solver in Matlab?
Hi, I would like to write another ode solver, similar to ode45 and ode23. I have the code that takes the timestep, but what I have no idea on how to do is to write a function that behaves (in the i/o sense) like ode**. I looked at ode45.m and if I were to start modifying those it’ll take me several days. What I am wondering if there’s a template that does all this stuff for me and all *I* have to do is program the time step. Thanks!
