21.04.2015 Views

The Computable Differential Equation Lecture ... - Bruce E. Shapiro

The Computable Differential Equation Lecture ... - Bruce E. Shapiro

The Computable Differential Equation Lecture ... - Bruce E. Shapiro

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!