how to reset/change the initial data?
9 views (last 30 days)
Hi there, I am working on a non linear system of 4 ordinary differential equations with time varying coefficients. I solved them numerically using ode 15s. It worked well. But now i want to stop my solver every 2 or 3 years and change/reset the initial data and run it again. I dont know how to do it. Any guidance/help is greatly appreciated.
Star Strider on 26 Jul 2012
Edited: Star Strider on 26 Jul 2012
You can give the ODE solvers a vector of time points at which to evaluate your ODE function. I've done this often, especially because it constrains the integration time without significantly affecting the accuracy of the integration.
From the documentation for the ODE solvers http://www.mathworks.com/help/techdoc/ref/ode23.html :
tspan — A vector specifying the interval of integration, [t0,tf]. The solver imposes the initial conditions at tspan(1), and integrates from tspan(1) to tspan(end). To obtain solutions at specific times (all increasing or all decreasing), use tspan = [t0,t1,...,tf].
For tspan vectors with two elements [t0 tf], the solver returns the solution evaluated at every integration step. For tspan vectors with more than two elements, the solver returns solutions evaluated at the given time points. The time values must be in order, either increasing or decreasing.
I also refer you to http://www.mathworks.com/help/techdoc/ref/odeset.html.