1 / 92

Solutions to ODEs

Solutions to ODEs. A general form for a first order ODE is Or alternatively. A general form for a first order ODE is Or alternatively We desire a solution y(x) which satisfies (1) and one specified boundary condition. .

Jimmy
Download Presentation

Solutions to ODEs

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Solutions to ODEs

  2. A general form for a first order ODE is • Or alternatively

  3. A general form for a first order ODE is • Or alternatively • We desire a solution y(x) which satisfies (1) and one specified boundary condition.

  4. To do this we divide the interval in the independent variable x for the interval [a,b] into subintervals or steps.

  5. To do this we divide the interval in the independent variable x for the interval [a,b] into subintervals or steps. • Then the value of the true solution is approximated at n+1 evenly spaced values of x. • Such that, • where,

  6. We denote the approximation at the base pts by so that

  7. We denote the approximation at the base pts by so that • The true derivation dy/dx at the base points is approximated by or just • where

  8. We denote the approximation at the base pts by so that • The true derivation dy/dx at the base points is approximated by or just • where • Assuming no roundoff error the difference in the calculated and true value is the truncation error,

  9. The only errors which will be examined are those inherent to the numerical methods.

  10. Numerical algorithms for solving 1st order odes for an initial condition are based on one of two approaches: • Direct or indirect use of the Taylor series expansion for the solution function • Use of open or closed integration formulas.

  11. The various procedures are classified into: • One-step: calculation of given the differential equation and • Multi-step: in addition the previous information they require values of x and y outside of the interval under consideration

  12. The first method considered is the Taylor series method. • It forms the basis for some of the other methods.

  13. Taylor Expansion

  14. For this method we express the solution about some starting point using a Taylor expansion.

  15. where

  16. Where • And • etc.

  17. Consider the example, where • And initial condition

  18. Consider the example, where • And initial condition • Differentiating and then applying initial conditions,

  19. Consider the example, where • And initial condition • Differentiating and then applying initial conditions,

  20. Consider the example, where • And initial condition • Differentiating and then applying initial conditions,

  21. Consider the example, where • And initial condition • Differentiating and then applying initial conditions,

  22. Continuing,

  23. Continuing,

  24. Continuing, . . . . Here the third derivative and higher are zero.

  25. Substituting these values into Taylor expansion gives,

  26. Substituting these values into Taylor expansion gives, • To determine the error lets consider compare this to the analytic solution.

  27. Substituting these values into Taylor expansion gives, • To determine the error lets consider compare this to the analytic solution. • To do this we will integrate the differential equation.

  28. Integrating,

  29. Integrating,

  30. Integrating,

  31. Integrating, • Simplifying,

  32. In this case the two expressions are identical, therefore there is no truncation error.

  33. Homework (Due Thursday) • Consider the function • Initial condition • Write the Taylor expansion and show that it can be written in the form, • Show if terms up to are retained, the error is,

  34. Stepping from to follows from the Taylor expansion of about .

  35. Stepping from to follows from the Taylor expansion of about . • Equivalently stepping from to can be accomplished from the Taylor expansion of about .

  36. Stepping from to follows from the Taylor expansion of about . • Equivalently stepping from to can be accomplished from the Taylor expansion of about .

  37. Algorithms for which the last term in expansion is dropped are of order hn. • The error is or order hn+1. • The local truncation error is bounded as follows • where

  38. However differentiation of can be complicated.

  39. However differentiation of can be complicated. • Direct Taylor expansion is not used other than the simplest case, • Big Oh is the order of the algorithm.

  40. Usually only the value of is the only value of that is known, therefore must be replaced by . • Thus, • In general, which is known as Euler’s method.

  41. Euler’s Method

  42. The problem with the simple Euler method is the inherent inaccuracy in the formula.

  43. y y1 y(x) y(x0) y(x1 h x0 x1 x • The geometric interpretation is shown the diagram.

  44. y y1 y(x) y(x0) y(x1 h x0 x1 x • The solution across the interval[x0,x1] is assumed to follow the line tangent to y(x) at x0.

  45. When the method is applied repeatedly across several intervals in sequence, the numerical solution traces out a polygon segment with sides of slope fi, i=0,1,2,…,(n-1).

  46. The simple Euler method is a linear approximation, and only works if the function is linear (or at least linear in the interval).

  47. The simple Euler method is a linear approximation, and only works if the function is linear (or at least linear in the interval). • This is inherently inaccurate. • Because of this inaccuracy small step sizes are required when using the algorithm.

  48. From the graph, as h -> 0, the approximation becomes better because the curve in the region becomes approx. linear. y y1 y(x) y(x0) y(x1 h x0 x1 x

  49. The modified Euler method improves on the simple Euler method.

  50. Modified Euler

More Related