Solving Initial-Value Problem for a System of ODE’s
Extended from the previous section Solving Initial-Value Problem for ODE, all the numerical methods are the same as described before except that they are in vector forms and the absolute value involved with the estimate of local truncation error is changed to a vector norm.
In addition to the numerical feature as the previous section, a new function of illustrating the dynamic feature of a autonomous ODE system through phase portrait is added and shown in the figure of phase plot.
In one of menu items, examples, there are several systems of ODE’s
Listed as examples. Among them, Lotka-Volterra, pendulum, Van der Pol, Duffing, Lorenz equations are selected to illustrate their dynamic features; linear1, linear2, linear3 are selected to illustrate the numerical phenomenon of stiff ODE systems.
Lotka-Volterra equation:
, and .
Pendulum equation:
.
Van der Pol equation:
.
Duffing equation:
.
Lorenz equation:
, , and .
Linear 1 equation:
, and .
with exact solution , and .
Linear 2 equation:
, and .
with exact solution
, and .
Linear 3 equation:
, , and .
with exact solution
, , and .
References:
【1】
R. L. Burden and J. D. Faires, Numerical Analysis, PWS,
【2】
S. H. Strogatz, Nonlinear
Dynamics and Chaos,