General Approach

Most of the methods seen above can be applied to system of ODE's, or, indeed, higher order ODE's. The following approach can easily be extended to an ODE of any degree or to a system of more than 2 ODE's.

If the initial value problem is given by an ODE of higher order, one can rewrite it as a system of first order ODE's. For example, the initial value problem

[Maple Math] ,

can be written as a system of linear ODE's by the introduction of the variable [Maple Math] :

[Maple Math] ,

[Maple Math] ,

[Maple Math] .

>

The coefficients [Maple Math] may be functions of [Maple Math] ! This system can be rewritten in the us more familiar form

[Maple Math] ,

[Maple Math] ,

[Maple Math] .

>

We can generalise the previous methods to solve a system of ODE's of the form above. For example, the family of Taylor series methods and Runge-Kutta methods can be generalised as

>

[Maple Math] ,

[Maple Math]

[Maple Math] , [Maple Math]

where the [Maple Math] are approximations for the values [Maple Math] and [Maple Math] respectively.

In case of the Runge-Kutta method of order four, this translates as:

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math] ,

[Maple Math]

>