The Computable Differential Equation Lecture ... - Bruce E. Shapiro
The Computable Differential Equation Lecture ... - Bruce E. Shapiro The Computable Differential Equation Lecture ... - Bruce E. Shapiro
40 CHAPTER 2. SUCCESSIVE APPROXIMATIONS Hence by equation 2.82, |φ(t) − φ n (t)| = ∣ ≤ ≤ ≤ M K ∞∑ i=n+1 ∞∑ i=n+1 ∞∑ i=n+1 (φ i (t) − φ i−1 (t)) ∣ (2.90) |φ i (t) − φ i−1 (t)| (2.91) K i−1 M |t − t 0 | i (2.92) i! ∞∑ i=n+1 |K(t − t 0 )| i Therefore by comparison with a Taylor series for e K(b−a) , i! (2.93) ‖φ(t) − φ n (t)‖ ≤ M K ≤ M K ∞∑ |K(b − a)| i i! i=n+1 ( ) n∑ e K(b−a) |K(b − a)| i − i! i=0 (2.94) (2.95) ≤ M K sup R n (t) (2.96) 0≤t≤KL where R n (t) is the Taylor Series Remainder for e t after n terms, sup R n (t) ≤ 0≤t≤KL sup 0≤{c,t}≤KL t n+1 e c (n + 1)! ≤ (KL)n+1 e KL (n + 1)! (2.97) for some unknown c between a and b.Hence proving the theorem. ‖φ(t) − φ n (t)‖ ≤ M K (KL) n+1 e KL (n + 1)! (2.98) The following example shows that this bounds is not very useful in practice. Example 2.7. Estimate the number of iterations required to obtain an solution to y ′ = t, y(0) = 1 on [0, 10] with a precision of no more that 10 −7 . Solution. Since f(t, y) = t we have f y = 0 and hence a Lipshitz constant is K = 1 (or any positive number), and we can use M = 10 on [0, 10]. The precision in the error is bounded by ɛ ≤ M(KL)n+1 e KL K(n + 1)! ≤ 10(10)n+1 e 10 (n + 1)! (2.99) We can determine the minimum value of n by using Mathematica. The following will print a list of values of 2.99 for n ranging from 1 to 50. Math 582B, Spring 2007 California State University Northridge c○2007, B.E.Shapiro Last revised: May 23, 2007
CHAPTER 2. SUCCESSIVE APPROXIMATIONS 41 errs = Table[{n, 10 (10) ˆ (n + 1) (E ˆ 10.)/(n + 1)!}, {n, 1, 50}] The output is a list of number pairs, which can be plotted with ListPlot or
- Page 1 and 2: The Computable Differential Equatio
- Page 3 and 4: Contents 1 Classifying The Problem
- Page 5 and 6: CONTENTS v Timeline on Computable D
- Page 7 and 8: Chapter 1 Classifying The Problem 1
- Page 9 and 10: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 11 and 12: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 13 and 14: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 15 and 16: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 17 and 18: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 19 and 20: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 21 and 22: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 23 and 24: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 25 and 26: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 27 and 28: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 29 and 30: CHAPTER 1. CLASSIFYING THE PROBLEM
- Page 31 and 32: Chapter 2 Successive Approximations
- Page 33 and 34: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 35 and 36: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 37 and 38: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 39 and 40: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 41 and 42: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 43 and 44: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 45: CHAPTER 2. SUCCESSIVE APPROXIMATION
- Page 49 and 50: Chapter 3 Approximate Solutions 3.1
- Page 51 and 52: CHAPTER 3. APPROXIMATE SOLUTIONS 45
- Page 53 and 54: CHAPTER 3. APPROXIMATE SOLUTIONS 47
- Page 55 and 56: CHAPTER 3. APPROXIMATE SOLUTIONS 49
- Page 57 and 58: CHAPTER 3. APPROXIMATE SOLUTIONS 51
- Page 59 and 60: CHAPTER 3. APPROXIMATE SOLUTIONS 53
- Page 61 and 62: CHAPTER 3. APPROXIMATE SOLUTIONS 55
- Page 63 and 64: CHAPTER 3. APPROXIMATE SOLUTIONS 57
- Page 65 and 66: CHAPTER 3. APPROXIMATE SOLUTIONS 59
- Page 67 and 68: CHAPTER 3. APPROXIMATE SOLUTIONS 61
- Page 69 and 70: CHAPTER 3. APPROXIMATE SOLUTIONS 63
- Page 71 and 72: CHAPTER 3. APPROXIMATE SOLUTIONS 65
- Page 73 and 74: CHAPTER 3. APPROXIMATE SOLUTIONS 67
- Page 75 and 76: Chapter 4 Improving on Euler’s Me
- Page 77 and 78: CHAPTER 4. IMPROVING ON EULER’S M
- Page 79 and 80: CHAPTER 4. IMPROVING ON EULER’S M
- Page 81 and 82: CHAPTER 4. IMPROVING ON EULER’S M
- Page 83 and 84: CHAPTER 4. IMPROVING ON EULER’S M
- Page 85 and 86: CHAPTER 4. IMPROVING ON EULER’S M
- Page 87 and 88: CHAPTER 4. IMPROVING ON EULER’S M
- Page 89 and 90: CHAPTER 4. IMPROVING ON EULER’S M
- Page 91 and 92: CHAPTER 4. IMPROVING ON EULER’S M
- Page 93 and 94: CHAPTER 4. IMPROVING ON EULER’S M
- Page 95 and 96: Chapter 5 Runge-Kutta Methods 5.1 T
CHAPTER 2. SUCCESSIVE APPROXIMATIONS 41<br />
errs = Table[{n, 10 (10) ˆ (n + 1) (E ˆ 10.)/(n + 1)!}, {n, 1,<br />
50}]<br />
<strong>The</strong> output is a list of number pairs, which can be plotted with ListPlot or<br />