2.5.20 ode: - ordinary differential equation solver CALLING SEQUENCE y=ode(y0,t0,t,f) [y,w,iw]=ode([type],y0,t0,t [,rtol [,atol]],f [,jac] [,w,iw]) [y,rd,w,iw]=ode("root",y0,t0,t [,rtol [,atol]],f [,jac],ng,g [,w,iw]) y=ode("discrete",y0,k0,kvect,f) PARAMETERS - : real vector or matrix (initial conditions).
- : real scalar (initial time).
- t : real vector (times at which the solution is computed).
- f : external (function or character string or list).
- type : one of the following character string: "adams" "stiff" "rk" "rkf" "fix" "discrete" "roots"
- rtol,atol : real constants or real vectors of the same size as y
- jac : external (function or character string or list).
- w,iw : real vectors.
- ng : integer.
- g : external (function or character string or list).
- : integer (initial time). kvect : integer vector.
DESCRIPTION ode is the standard function for solving explicit ODE systems defined by: dy/dt=f(t,y) , y(t0)=y0. It is an interface to various solvers, in particular to ODEPACK. The type of problem solved and the method used depend on the value of the first optional argument type which can be one of the following strings: - not given lsoda solver of package ODEPACK is called by default. It automatically selects between nonstiff predictor-corrector Adams method and stiff Backward Differentiation Formula (BDF) method. It uses nonstiff method initially and dynamically monitors data in order to decide which method to use.
| |
|