The Computable Differential Equation Lecture ... - Bruce E. Shapiro
The Computable Differential Equation Lecture ... - Bruce E. Shapiro
The Computable Differential Equation Lecture ... - Bruce E. Shapiro
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 2. SUCCESSIVE APPROXIMATIONS 39<br />
Where K is any number larger than sup (a,b) f y . If we choose the endpoints a and<br />
b such that |b − a| < 1/K we have K|b − a| < 1. Thus T is a contraction. By the<br />
contraction mapping theorem it has a fixed point; call this point φ. <strong>Equation</strong> 2.73<br />
follows immediately.<br />
<strong>The</strong>orem 2.12 (Error Bounds on Picard Iterations). Under the same conditions<br />
as before, let φ n be the n th Picard iterate, and let φ be the solution of the IVP. <strong>The</strong>n<br />
|φ(t) − φ n (t)| ≤ M |K(t − t 0)| n+1<br />
e K|t−t 0|<br />
K(n + 1)!<br />
(2.80)<br />
where M = sup D |f(t, y)| and K is a Lipshitz constant. Furthermore, if L = |b − a|<br />
then<br />
‖φ(t) − φ n (t)‖ ≤ M [KL]n+1 e KL<br />
(2.81)<br />
K(n + 1)!<br />
where ‖ · ‖ denotes the sup-norm.<br />
Proof. We begin by proving the conjecture<br />
For n = 1, equation 2.82 says that<br />
|φ n − φ n−1 | ≤ Kn−1 M<br />
|t − t 0 | n (2.82)<br />
n!<br />
|φ 1 − y 0 | ≤ M|t − t 0 | (2.83)<br />
which follows immediately from equation 2.73. Next, make the inductive hypothesis<br />
2.82 and calculate<br />
∫ t<br />
|φ n+1 − φn| =<br />
∣ [f(s, φ n (s)) − f(s, φ n−1 (s))] ds<br />
∣ (2.84)<br />
t 0<br />
≤ K<br />
∫ t<br />
t 0<br />
|φ n (s) − φ n−1 (s)| ds (2.85)<br />
by the definition of φ n and the Lipshitz condition. Applying the inductive hypothesis<br />
and then integrating,<br />
|φ n+1 − φn| ≤ Kn M<br />
n!<br />
which proves conjecture 2.82. Now let<br />
≤<br />
φ n (t) = φ 0 (t) +<br />
∫ t<br />
t 0<br />
|s − t 0 | n ds (2.86)<br />
Kn M<br />
(n + 1)! |t − t 0| n+1 (2.87)<br />
n∑<br />
[φ i (t) − φ i−1 (t)] (2.88)<br />
<strong>The</strong>n since the sequence of Picard iterates converges to the solution,<br />
i=1<br />
φ(t) = lim<br />
n→∞ φ n(t) = φ 0 (t) +<br />
∞∑<br />
[φ i (t) − φ i−1 (t)] (2.89)<br />
i=1<br />
c○2007, B.E.<strong>Shapiro</strong><br />
Last revised: May 23, 2007<br />
Math 582B, Spring 2007<br />
California State University Northridge