a logo for the course

Numerical Analysis

David Gleich

Purdue University

Spring 2016

Course number CS-51400, MATH-51400

Tuesday and Thursday, noon-1:15pm

Location Lawson B155

Homework 8

Please answer the following questions in complete sentences in a typed manuscript and submit the solution on blackboard by on April 25th at noon. Submissions will be accepted and graded at any point before the final exam is given, however, any submissions turned in after April 25th cannot be guaranteed to be graded before the final exam.

Problem 0: Homework checklist

Problem 1: (Ascher and Greif, Problem 16.5)

Consider the ODE with .

  1. Apply the stretching transformation to obtain the equivalent ODE And show the relationship between and .

  2. Show that applying the forward Euler method to the ODE for in with step size is equivalent to applying the same method to the ODE for in with step size that satisfies .

Problem 2: (Ascher and Greif, Example 16.20)

Consider the following system that describes the behavior of two masses in the presence of a third mass of (small) size. The masses are and (earth and sun): The initial conditions are Let .

  1. Implement a 4th order RK method for this problem.

  2. Show what happens for , and , , steps. Discuss any interesting observations.

  3. (Not required) Use a standard integration software package such as ode45 or an equivalent to look at the result.

Problem 3: (Gautschi Machine Exercise 5.5)

Newton's equations for the motion of a particle on a planar orbit (with eccentricity , ) are

  1. Verify that the solution can be written in the form of
  2. Write a program to solve the nonlinear equation at each step using Newton's method and plot the exact solution for .
  3. Write a forward and backwards Euler solver for this problem.
  4. Evaluate the error with and time-steps using the same values of .