We first of all derive a method for solving this system when the given differential equation is linear,

i.e.

is of the form

where and are given functions. Defining the diagonal matrix

and the vector

we may write

The system (10.27) now reduces to the system of linear equations,

(10.28)

where

(10.29)