EDO Problemas

Download as pdf or txt
Download as pdf or txt
You are on page 1of 3

236CHAPTER 5.

NUMERICAL SOLUTION OF INITIAL-VALUE PROBLEMS

EXERCISE SET 5.3

1. Use the Midpoint method to approximate the solutions to each of the following
initial-value problems, and compare the results to the actual values.

(a) y  = te3t −2y, for 0 ≤ t ≤ 1, with y(0) = 0 and h = 0.5; actual solution
y(t) = 15 te3t − 25
1 3t
e + 251 −2t
e .

(b) y  = 1 + (t − y)2 , for 2 ≤ t ≤ 3, with y(2) = 1 and h = 0.5; actual


solution y(t) = t + 1/(1 − t).
y
(c) y  = 1 + , for 1 ≤ t ≤ 2, with y(1) = 2 and h = 0.25; actual solution
t
y(t) = t ln t + 2t.

(d) y  = cos 2t + sin 3t, for 0 ≤ t ≤ 1, with y(0) = 1 and h = 0.25; actual
solution y(t) = 12 sin 2t − 13 cos 3t + 43 .

2. Repeat Exercise 1 using Heun’s method.

3. Repeat Exercise 1 using the Modified Euler method.

4. Use the Modified Euler method to approximate the solutions to each of the
following initial-value problems, and compare the results to the actual values.

y # y $2
(a) y  = − , for 1 ≤ t ≤ 2, with y(1) = 1 and h = 0.1; actual
t t
solution y(t) = t/(1 + ln t).
y # y $2
(b) y  = 1 + + , for 1 ≤ t ≤ 3, with y(1) = 0 and h = 0.2; actual
t t
solution y(t) = t tan(ln t).

(c) y  = −(y + 1)(y + 3), for 0 ≤ t ≤ 2, with y(0) = −2 and h = 0.2;


actual solution y(t) = −3 + 2(1 + e−2t )−1 .
1
(d) y  = −5y + 5t2 + 2t, for 0 ≤ t ≤ 1, with y(0) = 3 and h = 0.1; actual
solution y(t) = t2 + 13 e−5t .

5. Use the results of Exercise 4 and linear interpolation to approximate values


of y(t), and compare the results to the actual values.

(b) y(2.1) and y(2.75)


(a) y(1.25) and y(1.93)
(c) y(1.3) and y(1.93) (d) y(0.54) and y(0.94)

6. Repeat Exercise 4 using Heun’s method.

7. Repeat Exercise 5 using the results of Exercise 6.


5.3. RUNGE-KUTTA METHODS 237

8. Repeat Exercise 4 using the Midpoint method.


9. Repeat Exercise 5 using the results of Exercise 8.
10. Repeat Exercise 1 using the Runge-Kutta method of order 4.
11. Repeat Exercise 4 using the Runge-Kutta method of order 4.
12. Use the results of Exercise 11 and Cubic Hermite interpolation to approximate
values of y(t) and compare the approximations to the actual values.

(b) y(2.1) and y(2.75)


(a) y(1.25) and y(1.93)
(c) y(1.3) and y(1.93) (d) y(0.54) and y(0.94)

13. Show that the Midpoint method, the Modified Euler method, and Heun’s
method give the same approximations to the initial-value problem
y  = −y + t + 1, 0 ≤ t ≤ 1, y(0) = 1,
for any choice of h. Why is this true?
14. Water flows from an inverted conical tank with circular orifice at the rate

dx 2
x
= −0.6πr −2g ,
dt A(x)
where r is the radius of the orifice, x is the height of the liquid level from
the vertex of the cone, and A(x) is the area of the cross section of the tank
x units above the orifice. Suppose r = 0.1 ft, g = −32.17 ft/s2 , and the tank
has an initial water level of 8 ft and initial volume of 512(π/3) ft3 .

(a) Compute the water level after 10 min with h = 20 s.


(b) Determine, to within 1 min, when the tank will be empty.

15. The irreversible chemical reaction in which two molecules of solid potassium
dichromate (K2 Cr2 O7 ), two molecules of water (H2 O), and three atoms of
solid sulfur (S) combine to yield three molecules of the gas sulfur dioxide
(SO2 ), four molecules of solid potassium hydroxide (KOH), and two molecules
of solid chromic oxide (Cr2 O3 ) can be represented symbolically by the stoi-
chiometric equation
2K2 Cr2 O7 + 2H2 O + 3S −→ 4KOH + 2Cr2 O3 + 3SO2 .
If n1 molecules of K2 Cr2 O7 , n2 molecules of H2 O, and n3 molecules of S are
originally available, the following differential equation describes the amount
x(t) of KOH after time t:
#  3
dx x $2 # x $2 3x
= k n1 − n2 − n3 − ,
dt 2 2 4
238CHAPTER 5. NUMERICAL SOLUTION OF INITIAL-VALUE PROBLEMS

where k is the velocity constant of the reaction. If k = 6.22 × 10−19 , n1 =


n2 = 2 × 103 , and n3 = 3 × 103 , how many units of potassium hydroxide will
have been formed after 0.2 s?

You might also like