Tutorial Problems #6
MAT 292 – Calculus III – Fall 2015
Solutions
3.2 - #31 Consider two interconnected tanks similar to those in Figure 3.2.9. Tank 1 initially contains 60 gal
of water and Q01 oz of salt, and Tank 2 initially contains 100 gal of water and Q02 oz of salt. Water containing q1
oz/gal of salt flows into Tank 1 at a rate of 3 gal/min. The mixture in Tank 1 flows out at a rate of 4 gal/min,
which half flows into Tank 2, while the remainder leaves the system. Water containing q2 oz/gal of salt also
flows into Tank 2 from the outside at the rate of 1 gal/min. The mixture in Tank 2 leaves the tank at a rate of
3 gal/min, of which some flows back into Tank 1 at a rate of 1 gal/min, while the rest leaves the system.
(a) Draw a diagram that depicts the flow process described above. Let Q1 t) and Q2 (t), respectively, be the
amount of salt in each tank at time t. Write down di↵erential equations and initial conditions for Q1 and
Q2 that model the flow process
(b) Find the equilibrium values QE E
1 and Q2 in terms of the concentrations q1 and q2 .
(c) Is it possible (by adjusting q1 and q2 ) to obtain QE E
1 = 60 and Q2 = 50 as an equilibrium state?
(d) Describe which equilibrium states are possible for this system for various values of q1 and q2
Solution
(a) We see the diagram to be given by
1
Tutorial #6 – Fall 2015 MAT 292
The DE’s are easily seen to be
4 1
Q01 = Q1 + Q2 + 3q1 Q1 (0) = Q01
60 100
2 3
Q02 = Q1 Q2 + q2 Q2 (0) = Q02
60 100
(b) The equilibrium is easily found setting Q01 = Q02 = 0, we see
! ! ! ! !
1 4/3 1/5 3q1 18 6 3q1 54q1 + 6q2
QE = =) QE = =
20 2/3 3/5 q2 20 40 q2 60q1 + 40q2
(c) Is it isn’t possible to obtain QE E
1 = 60 and Q2 = 50 as an equilibrium state since it’d require a negative
q2 . Since
2 3 1
⇤ 60 ⇤ 50 + q2 = 0 =) q2 =
60 100 2
(d) The states that we’d allow must satisfy
2 E 3 E 4 E 1 E
Q Q 60 & Q1 + Q 60
60 1 100 2 180 300 2
3 E 9 E
=) Q 6 QE
1 6 Q
20 2 10 2
3.3 - # 17-20 Consider x0 = Ax. If given the eigenvectors and eigenvalues:
(a) Sketch a phase portrait of the system.
(b) Sketch the trajectory passing through the initial point (2,3)
(c) For the trajectory in part b), sketch the component plots of x1 versus t and of x2 versus t on the same
set of axes.
#17
! ! ! !
~1 = 1 ~2 = 1 1 t 1 2t
1 = 1, & 2 = 2, =) x(t) = C1 e + C2 e
2 2 2 2
with a portrait like
2
Tutorial #6 – Fall 2015 MAT 292
#18 ! ! ! !
~1 = 1 ~2 = 1 1 1
1 = 1, & 2 = 2, =) x(t) = C1 et + C2 e 2t
2 2 2 2
with a portrait like
#19 ! ! ! !
~1 = 1 ~2 = 1 1 t 1
1 = 1, & 2 = 2, =) x(t) = C1 e + C2 e2t
2 2 2 2
with a portrait like
#20 ! ! ! !
~1 = 1 ~2 = 1 1 t 1
1 = 1, & 2 = 2, =) x(t) = C1 e + C2 e2t
2 2 2 2
with a portrait like
3
Tutorial #6 – Fall 2015 MAT 292
3.3 - #28 Consider ! ! !
d I R1 /L 1/L I
=
dt V 1/C 1/CR2 V
(a) Find a condition on R1 , R2 , C and L that must be satisfied if the eigenvalues of the efficient matrix are
to be real and di↵erent.
(b) If the condition found in part a) is satisfied, show that both eigenvalues are negative. Then show that
I(t) ! 0 and V (t) ! 0 as t ! 1, regardless of the initial conditions.
Solution
(a) Let’s compute the characteristic equation. We see
✓ ◆
+ R1 /L 1/L 2 CR1 R2 + L R 1 + R2
P ( ) = det(1 A) = = + +
1/C + 1/CR2 CLR2 CLR2
Now the eigenvalues are given by the roots, but we know the discriminate needs to be positive and non-zero
for the roots to be real and distinct. i.e
✓ ◆2 ✓ ◆2
CR1 R2 + L R1 + R 2 R2 R1 L 1 1 R1 1 4
b 2
4ac = 4 = 21 2 2
+ 2 2 4 = >0
CLR2 CLR2 L CL R2 C R2 CL L CR2 CL
(b) The actual eigenvalues are given by
✓ ◆ s✓ ◆2
1 R1 1 1 R1 1 4
± = + ±
2 L CR2 2 L CR2 CL
Clearly is negative, and + is negative since
✓ ◆ s✓ ◆2
1 R1 1 1 R1 1 4
+ = + +
2 L CR2 2 L CR2 CL
✓ ◆ s✓ ◆2
1 R1 1 1 R1 1
6 + +
2 L CR2 2 L CR2
✓ ◆
1 R1 1 1 R1 1
= + +
2 L CR2 2 L CR2
<0
4
Tutorial #6 – Fall 2015 MAT 292
We know the solution takes the form
!
I
= C1 ~ + e +t
+ C2 ~ e t
V
Therefore the solution goes to zero as t ! 1 since the eigenvalues are negative regardless of the initial
data.
3.4 - #7 Solve the following system, draw direction field and a phase portrait. Describe the behaviour of
solutions as t ! 1
! !
0 1 4 4
x = x, x(0) =
1 1 3
Solution Let’s compute the characteristic equation to find the eigenvalues.
+1 4 2
P ( ) = det(1 A) = = + 2 + 5 = 0 =) ± = 1 ± 2i
1 +1
Now that we’ve found the eigenvalues, let’s compute the eigenvectors.
! ! !
2i 4 2i 2i
ker(1 + A) = ker = span =) ~ + =
1 2i 1 1
The conjugate gives us the other vector. !
~ 2i
=
1
Thus the complex valued solution is given by
h i
x(t) = e t
C1 ~ e 2it
+ C2~ + e2it
Using Euler’s identity and assuming C2 = C 1 we obtain
" ! !#
t 2 sin(2t) 2 cos(2t)
x(t) = e A +B
cos(2t) sin(2t)
The phase portrait is given by
5
Tutorial #6 – Fall 2015 MAT 292
3.4 - #22 The electric circuit shown below is described by the system of di↵erential equations
! ! !
d I 0 1/L I
=
dt V 1/C 1/RC V
where I is the current through the inductor and V is the voltage across the capacitor.
(a) Show that the eigenvalues of the coefficient matrix are real and di↵erent if L > 4R2 C. Show that they
are complex conjugates if L < 4R2 C.
(b) Suppose that R = 1⌦, C = 1/2 F, and L = 1 H. Find the general solution of the system in this case.
(c) Find I(t) and V (t) if I(0) = 2 Amp and V (0) = 1 Volt.
Solution
(a) This follows from our previous calculation setting R1 = 0 and R2 = R, we see the eigenvalues are real and
distinct if
1 4
> 0 =) L > 4R2 C
C 2 R2 CL
and complex conjugates if the discriminate is negative, i.e. L < 4R2 C
(b) Suppose that R = 1⌦, C = 1/2 F, and L = 1 H, we see the system is
! ! !
d I 0 1 I
=
dt V 2 2 V
It’s easy to see that
1 2
P ( ) = det(A 1 )= = +2 +2
2 2
Thus the eigenvalues are
± = 1±i
The eigenvectors are found by checking the kernel.
! ! !
1 i 1 1+i 1+i
+ =) ker = span =) ~ + =
2 1 i 2 2
We may take the complex conjugate to retrieve the second
!
~ 1 i
= + =
2
Thus the complex valued solution is given by
!
I h i
= e t C1~ + eit + C2~ e it
C 1 , C2 2 C
V
To make this real valued, fix C2 = C 1 . Then using Euler’s identity we obtain
! " ! !#
I sin(t) sin(t) + cos(t)
=e t A +B , A, B 2 R
V cos(t) sin(t) 2 sin(t)
6
Tutorial #6 – Fall 2015 MAT 292
(c) Here we have I(0) = 2 and V (0) = 1, so we solve for A and B:
! " ! !#
2 0 1
= A +B =) A = 1 & B=2
1 1 0
Thus the solution is given by ! !
I t 3 sin(t) + 2 cos(t)
=e
V cos(t) 5 sin(t)
