To solve the system of differential equations

(1)

where is a matrix and and are vectors, first consider the homogeneous case with . The solutions to

(2)

are given by

(3)

But, by the matrix decomposition theorem, the matrix exponential can be written as

(4)

where the eigenvector matrix is

(5)

and the eigenvalue matrix is

(6)

Now consider

			(7)
			(8)
			(9)
			(10)

The individual solutions are then

(11)

so the homogeneous solution is

(12)

where the s are arbitrary constants.

The general procedure is therefore

1. Find the eigenvalues of the matrix (, ..., ) by solving the characteristic equation.

2. Determine the corresponding eigenvectors , ..., .

3. Compute

(13)

for , ..., . Then the vectors which are real are solutions to the homogeneous equation. If is a matrix, the complex vectors correspond to real solutions to the homogeneous equation given by and .

4. If the equation is nonhomogeneous, find the particular solution given by

(14)

where the matrix is defined by

(15)

If the equation is homogeneous so that , then look for a solution of the form

(16)

This leads to an equation

(17)

so is an eigenvector and an eigenvalue.

5. The general solution is

(18)

CITE THIS AS:

Weisstein, Eric W. "Ordinary Differential Equation--System with Constant Coefficients." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/OrdinaryDifferentialEquationSystemwithConstantCoefficients.html