nagle_differential_equations_ISM_Part4

nagle_differential_equations_ISM_Part4 - Use the method of...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Use the method of steps to show that the solution to the initial value problem x ( t ) =- x ( t- 1) , x ( t ) = 1 on [- 1 , 0] , is given by x ( t ) = n X k =0 (- 1) k [ t- ( k- 1)] k k ! , for n- 1 t n, where n is a nonnegative integer. (This problem can also be solved using the Laplace transform method of Chapter 7.) (c) Use the method of steps to compute the solution to the initial value problem given in (0.1) on the interval 0 t 15 for t = 3. Extrapolation When precise information about the form of the error in an approximation is known, a technique called extrapolation can be used to improve the rate of convergence. Suppose the approximation method converges with rate O ( h p ) as h 0 (cf. Section 3.6). From theoretical considerations, assume we know, more precisely, that y ( x ; h ) = ( x ) + h p a p ( x ) + O ( h p +1 ) , (0.3) where y ( x ; h ) is the approximation to ( x ) using step size h and a p ( x ) is some function that is independent of h (typically, we do not know a formula for a p ( x ), only that it exists). Our goal is to obtain approximations that converge at the faster rate O ( h p +1 ). We start by replacing h by h/ 2 in (0.3) to get y x ; h 2 = ( x ) + h p 2 p a p ( x ) + O ( h p +1 ) . If we multiply both sides by 2 p and subtract equation (0.3), we find 2 p y x ; h 2- y ( x ; h ) = (2 p- 1) ( x ) + O ( h p +1 ) . Solving for ( x ) yields ( x ) = 2 p y ( x ; h/ 2)- y ( x ; h ) 2 p- 1 + O ( h p +1 ) . Hence, y * x ; h 2 := 2 p y ( x ; h/ 2)- y ( x ; h ) 2 p- 1 has a rate of convergence of O ( h p +1 ). 10 (a) Assuming y * x ; h 2 = ( x ) + h p +1 a p +1 ( x ) + O ( h p +2 ) , show that y ** x ; h 4 := 2 p +1 y * ( x ; h/ 4)- y * ( x ; h/ 2) 2 p +1- 1 has a rate of convergence of O ( h p +2 ). (b) Assuming y ** x ; h 4 = ( x ) + h p +2 a p +2 ( x ) + O ( h p +3 ) , show that y *** x ; h 8 := 2 p +2 y ** ( x ; h/ 8)- y ** ( x ; h/ 4) 2 p +2- 1 has a rate of convergence of O ( h p +3 )....
View Full Document

This note was uploaded on 02/12/2011 for the course MA 221 taught by Professor Mazmani during the Spring '08 term at Stevens.

Page1 / 5

nagle_differential_equations_ISM_Part4 - Use the method of...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online