Differential Equations Project

Computer Lab Modules

 

Systems of D. E.'s

 

 

Background. We wish to study the system of D. E.'s
x' = a x + b y
y' = c x + d y

 

Computer Lab Work.

 

Exercise 1. Find the general solution to the system of D. E.'s
[Graphics:e12.txtgr1.gif]
Plot the solution curves where the starting points are (1,1), (-1,1), (-1,-1), (1,-1)
and the parameter t is in the interval 0 < = t < = 12.566

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr2.gif]

Put the D.E.'s in operator form and eliminate y to obtain a higher order D.E. for x, and find the roots of its characteristic equation.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr4.gif]

The roots are complex, [Graphics:e12.txtgr5.gif], so the general solution is formed as follows.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr6.gif]

Now construct the vector form of the solution R[t,c1,c2] and look at the initial value R[0,c1,c2].

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr7.gif]

Plot the solution curves where the starting points are (1,1), (-1,1), (-1,-1), (1,-1).

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr8.gif]
 

Exercise 2. Find the general solution to the system of D. E.'s
[Graphics:e12.txtgr9.gif]
Plot the solution curves where the starting points are (1,0), (0,1), (0.5,1), (1,0.5),
and the parameter t is in the interval 0 < = t < = 2.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr10.gif]

Put the D.E.'s in operator form and eliminate y to obtain a higher order D.E. for x, and find the roots of its characteristic equation.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr11.gif]

The roots are real and distinct, [Graphics:e12.txtgr12.gif], so the general solution is formed as follows.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr13.gif]

Now construct the vector form of the solution R[t,c1,c2] and look at the initial value R[0,c1,c2].

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr14.gif]

It is desirable to have the general solution in a form where (c1,c2) is an initial point that the solution curve passes through. So we adjust the above solution to fit this requirement

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr15.gif]

We will construct the vector function P so that the I.C.'s are easier to input.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr16.gif]

Plot the solution curves where the starting points are (1,0), (0,1), (0.5,1), (1,0.5).

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr17.gif]
 

Exercise 3. Find the general solution to the system of D. E.'s
[Graphics:e12.txtgr18.gif]
Plot the solution curves where the starting points are (1,1), (2,2), (3,3), (4,4).
and the parameter t is in the interval 0 < = t < = 5.9

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr19.gif]

Put the D.E.'s in operator form and eliminate y to obtain a higher order D.E. for x, and find the roots of its characteristic equation.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr20.gif]

The roots are pure complex, [Graphics:e12.txtgr21.gif], so the general solution is formed as follows.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr22.gif]

Now construct the vector form of the solution R[t,c1,c2] and look at the initial value R[0,c1,c2].

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr23.gif]

It is desirable to have the general solution in a form where (c1,c2) is an initial point that the solution curve passes through. So we adjust the above solution to fit this requirement

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr24.gif]

We will construct the vector function P so that the I.C.'s are easier to input.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr25.gif]

Plot the solution curves where the starting points are (1,1), (2,2), (3,3), (4,4).

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr26.gif]
 

Exercise 4. Find the general solution to the system of D. E.'s
[Graphics:e12.txtgr27.gif]
Plot the solution curves where the starting points are (1,0), (0,0.6), (0.6,0), (0,0.3).
and the parameter t is in the interval -0.2 < = t < = 1

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr28.gif]

Put the D.E.'s in operator form and eliminate y to obtain a higher order D.E. for x, and find the roots of its characteristic equation.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr29.gif]

The roots are real and distinct, [Graphics:e12.txtgr30.gif], so the general solution is formed as follows.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr31.gif]

Now construct the vector form of the solution R[t,c1,c2] and look at the initial value R[0,c1,c2].

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr32.gif]

It is desirable to have the general solution in a form where (c1,c2) is an initial point that the solution curve passes through. So we adjust the above solution to fit this requirement

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr33.gif]

We will construct the vector function P so that the I.C.'s are easier to input.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr34.gif]

Plot the solution curves where the starting points are (1,0), (0,0.6), (0.6,0), (0,0.3).

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr35.gif]
 

Exercise 5. Find the general solution to the system of D. E.'s
[Graphics:e12.txtgr36.gif]
Plot the solution curves where the starting points are (0.2,0), (0.6,0), (0,0.4), (0,1),
and the parameter t lies in the interval 2 < = t < = 0.8

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr37.gif]

Put the D.E.'s in operator form and eliminate y to obtain a higher order D.E. for x, and find the roots of its characteristic equation.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr38.gif]

The roots are real and equal, [Graphics:e12.txtgr39.gif], so the general solution is formed as follows.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr40.gif]

Now construct the vector form of the solution R[t,c1,c2] and look at the initial value R[0,c1,c2].

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr41.gif]

It is desirable to have the general solution in a form where (c1,c2) is an initial point that the solution curve passes through. So we adjust the above solution to fit this requirement

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr42.gif]

We will construct the vector function P so that the I.C.'s are easier to input.

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr43.gif]

Plot the solution curves where the starting points are (0.2,0), (0.6,0), (0,0.4), (0,1).

[Graphics:e12.txtgr3.gif][Graphics:e12.txtgr44.gif]
 

Solutions.

 



Return to the Differential Equations Project

Return to the Numerical Analysis Project

Return to the Complex Analysis Project

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews, 1998