This preview shows page 1. Sign up to view the full content.
Unformatted text preview: AIMS Exercise Set # 7
Peter J. Olver 1. In this exercise, you are asked to find "onesided" finite difference formulas for derivatives. These are useful for approximating derivatives of functions at or near the boundary of their domain. (a) Construct a second order, onesided finite difference formula that approximates the derivative f (x) using the values of f (x) at the points x, x + h and x + 2 h. (b) Find a finite difference formula for f (x) that involves the same values of f . What is the order of your formula? (c) Test your formulas by computing approximations to the first and second derivatives of f (x) = ex at x = 1 using step sizes h = .1, .01 and .001. What is the error in your numerical approximations? Are the errors compatible with the theoretical orders of the finite difference formulae? Discuss why or why not. (d) Answer part (c) at the point x = 0. 2. (a) Design an explicit numerical method for solving the initialboundary value problem ut = uxx + s(x), u(t, 0) = u(t, 1) = 0, u(0, x) = f (x), 0 x 1,
2 for the heat equation with a source term s(x). (b) Test your scheme on the particular problem for 1 1 2 x  6  3, f (x) = 0, 1 5 2 3 x 6 , 0 x 1, 3
1 3 2 3 2 x 3, = 1 , 6 s(x) = x(1  x)(10  22 x), x 1, using space step sizes h = .1 and .05, and a suitably chosen time step k. (c) What is the long term behavior of your solution? Can you find a formula for its eventual profile? (d) Design an implicit scheme for the same problem. Does the behavior of your numerical solution change? What are the advantages of the implicit scheme? 1 c 2006 Peter J. Olver ...
View
Full
Document
This note was uploaded on 02/10/2012 for the course MATH 5485 taught by Professor Olver during the Fall '09 term at University of Central Florida.
 Fall '09
 Olver
 Derivative, Formulas

Click to edit the document details