Numerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)

Size: px

Start display at page:

Download "Numerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)"

Spencer Douglas
5 years ago
Views:

1 Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation method)

2 We have seen in Section 1.1 that by drawing a direction field we can visualize qualitatively the behavior of solutions of a first order DE. When the grid is fine then we obtain a direction field that can convey a sense of the behavior of solutions for various initial conditions. For example, for DE we have direction field dy = y 2-2t dt For a first order DE we can sketch a direction field and visualize the behavior of solutions. This has the advantage of being a relatively simple process, even for complicated equations. However, direction fields do not lend themselves to quantitative computations or comparisons.

3 By choosing different initial conditions we can see approximately how some solutions behave. In this case we are considering initial value problems. Initial condition x = -2, y = -1 Initial condition x = -1, y = -1 To sketch these solutions of the IVP a numerical method was used.

4 For quantitative approximations of an IVP we need to consider numerical approximation schemes. Such techniques generate approximate values for the solution of the IVP at a set of discrete values of independent variable t; t 0, t 1, t 2, etc. So the output from such numerical algorithms consists of a set of ordered pairs (t k, y k ) where y k is an approximation to the true value of the solution at t k. (y k y(t k )) One algorithm that will allow you to estimate values of the solution y is given by Euler s Method sometimes called the tangent line method. (Euler (pronounced Oiler) lived in the 18th century is considered to be one of the greatest mathematicians who ever lived.) Observe that the direction field has many tangent line segments at successive values of t that almost touch each other. It takes only a bit of imagination to consider starting at a point specified by an initial condition and linking line segments for successive values of t in the grid, thereby producing a piecewise linear graph. Such a graph would provide an approximation to the solution of the IVP.

5 To see how Euler s method works, recall from calculus a section called Linear Approximation. In that section you did the following: Construct the tangent line to a differentiable function f(x) at a point (a, f(a)) and approximate values of the function at nearby values of x using the equation of the tangent line. Suppose we have the point (a, f(a)) and the slope of f(x) at this point is the value m. Then the tangent line is obtained using the point-slope method for the equation of a line. Tangent line is: y f(a) = m(x a) y = f(a) + m(x a). We can estimate the value of f(x) at x = a + Δx by computing f(a) + m(a + Δx a) = f(a) +m Δx. Euler s method uses this idea applied to an initial value problem as we illustrate next.

We start Euler s approximation to the IVP at the initial point (t 0, y 0 ) and use the fact that the slope of the exact solution at the initial point is f(t 0, y 0 ) so the equation of the tangent

6 We start Euler s approximation to the IVP at the initial point (t 0, y 0 ) and use the fact that the slope of the exact solution at the initial point is f(t 0, y 0 ) so the equation of the tangent line at the initial point is y y 0 = f(t 0, y 0 ) (t t 0 ). Rewrite this as y = y 0 + f(t 0, y 0 ) (t t 0 ). The tangent line is a good approximation to the actual solution curve on a short interval over which the slope of the solution does not change appreciably from its value at the initial point. Exact solution Thus, if t 1 is close enough to t 0, we can approximate the exact solution at t = t 1 by the value y 1 determined by substituting t 1 for t in the equation of the tangent line. Doing this we get y 1 = y 0 + f(t 0, y 0 ) (t 1 t 0 ). So far we have points (t 0, y 0 ) and (t 1, y 1 ). If we connect these points with a straight line we have a linear approximation to the exact solution φ(t) over interval [t 0, t 1 ]. linear approximation

7 To proceed further, we can try to repeat the process. Unfortunately, we do not know the value of the exact solution value at t 1. The best we can do is to use the approximate value y 1 instead. Thus we construct the line through (t 1, y 1 ) with the slope f(t 1, y 1 ). For a point t 2 close to t 1, we approximate the exact solution at t 2 using the line passing through (t 1, y 1 ) with slope f (t 1, y 1 ). We then compute y 2 = y 1 + f(t 1, y 1 ) (t 2 t 1 ) which approximates the exact solution at t 2. So now we have three points (t 0, y 0 ),(t 1, y 1 ), and (t 2, y 2 ). Using this process at the second step we are really approximating a new IVP; As along as the integral curve for the original IVP is not much different that than the integral curve for the new IVP, we can obtain an estimate of the solution of the original IVP. This is where some error creeps in. Thus we create a discrete sequence y n of approximations to the exact solution at the discrete points t 1, t 2, t 3, : where f k = f (t k, y k ) is the slope of the tangent line at t = t k for k =0, 1,, n.

8 The sequence of calculations generates the sequence of points To graph an Euler approximation, we plot the points (t 0, y 0 ), (t 1, y 1 ),, (t n, y n ), and then connect these points with line segments to provide a piecewise linear approximation to the solution of the original IVP. If the discrete values t 1, t 2, t 3, are equispaced we say the step size is a fixed value often denoted by h. For a fixed step size h = t k t k-1, for each value of k. Thus Euler s formula becomes

9 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' i = f(t i,y i ) = f i Slope computation y i+1 =y i + h f(t i,y i ) Estimate of soln at t i+1 Compute the slope at t 0 = 0 and use the fact that y 0 = y(0)= 1 in the Euler formula.

10 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' i = f(t i,y i ) = f i Slope computation y i+1 =y i + h f(t i,y i ) Estimate of soln at t i y 1 =y 0 + h f(t 0,y 0 ) = y1 y0 f0 h 1 ( ) ( 0. 2) ( 0. 2) 1. 5 Substituted t = 0 and y = 1 into formula for f(t, y)

11 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i Slope computation y i+1 =y i + h f(t i,y i ) Estimate of soln at t i y 1 =y 0 + h f(t 0,y 0 ) = y1 y0 f0 h 1 ( )( 0. 2) ( 0. 2) 1. 5 Compute the slope at t 1 = 0.2 and use the fact that we have the approximation y 1 = 1.5 to use in the Euler formula.

12 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y1 y0 f0 h 1 ( ) ( 0. 2) ( 0. 2) y2 y1 f1 h ( 0. 2) Substituted t = 0.2 and y = 1.5 into formula for f(t, y)

13 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y1 y0 f0 h 1 ( ) ( 0. 2) ( 0. 2) y2 y1 f1 h ( 0. 2) Compute the slope at t 2 = 0.4 and use the fact that we have the approximation y 2 = 1.87 to use in the Euler formula.

14 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y 3 =y 2 + h f(t 2,y 2 ) = y1 y0 f0 h 1 ( ) ( 0. 2) ( 0. 2) y2 y1 f1 h ( 0. 2) y3 y2 f2 h ( 0. 2)

15 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y 3 =y 2 + h f(t 2,y 2 ) = y y f h 1 ( )( 0. 2) ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) Continue in a similar fashion.

16 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y 3 =y 2 + h f(t 2,y 2 ) = y 4 =y 3 + h f(t 3,y 3 ) = y y f h 1 ( )( 0. 2) ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2)

17 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y 3 =y 2 + h f(t 2,y 2 ) = y 4 =y 3 + h f(t 3,y 3 ) = y y f h 1 ( )( 0. 2) ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2)

18 Example: For the initial value problem use Euler s method with h = 0.2 to approximate the solution at t = 0.2, 0.4, 0.6, 0.8, and 1.0. Index i t i y i y' I = f(t i,y i ) = f i y i+1 =y i + h f(t i,y i ) y 1 =y 0 + h f(t 0,y 0 ) = y 2 =y 1 + h f(t 1,y 1 ) = y 3 =y 2 + h f(t 2,y 2 ) = y 4 =y 3 + h f(t 3,y 3 ) = y 5 =y 4 + h f(t 4,y 4 ) = y y f h 1 ( )( 0. 2) ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2) y y f h ( 0. 2)

19 We can find the exact solution to our IVP is We approximated the solution of our IVP at t = 0.2, 0.4, 0.6, 0.8, 1.0 so now let s compare our results to the exact solution at those values of t. Exact y in red Approximate y in blue 2.5 y t

20 Euler Approx. 3 Direction Field for y'=3-2*t-.5*y Exact Note here we are looking over interval 0 t 1.

21 Accuracy of Euler s Method: Euler s method should not be used for real world approximations of IVPs. There are much more accurate methods available. But these other methods are used in much the same way as Euler s method. Euler s method is a recursive procedure used to obtain discrete approximations to IVPs. It generates a piecewise linear approximation by pasting together tangent line approximations. Its accuracy depends on how closely the t-values t 1, t 2, t 3, are spaced. For a very small spacing Euler s method requires a lot of steps and hence can be considered expensive. Even with small spacing certain IVPs can not be well approximated with Euler s method. Next we illustrate how Euler s method behaves for a sequence of decreasing values of the stepsize h. The IVP we use is 2 dy = y - 2t, y(-1.5) = -2.5 dt

22 h = 0.25

23 h = 0.25 h = 0.125

24 h = 0.25 h = h =

25 h = 0.25 h = h = h = 0.01

26 h = 0.25 h = h = h = 0.01 h = 0.001

27 Here is the direction field where we have used a much more accurate method called Runge-Kutta Fehlberg to generate the approximation to the solution of the IVP.

Numerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)

Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation