E1 consider the following initial value problem for a second-order equation:

y// = 2xy/ -  ,        1 < x < ,        y(1) = 1,    y/ (1) =  .

(a)  Rewrite this problem in standard form, as a system of first-order ODEs.

(b)  The Butcher tableau (Eq.  1) defines a Runge-Kutta method.  Apply this method with one step of size h =  to the initial value problem.

Show all your work explicitly; i.e., do the entire computation on paper, not with a Java/C program or a Maple procedure.  (Of course, you can use a pocket calculator or Maple for any floating point calculations required.)

0

1

2

(1)

E2  Show that the Runge-Kutta method (Eq.  1) is of order O(h3 ).  To that end, prove an

analogue of Theorem 6.1 in the lecture notes for this particular method. Restrict yourself to the case of a single first-order ODE, that is, to the case m = 1.

E3  The C code ode_rungekutta.c is a first-order ODE solver which applies a Runge-Kutta

algorithm. It takes information on the interval [a, b] along with initial condition y(a) = α , and number of steps n to perform. It then computes an approximate solution to the ODE y/ (x) = f (x, y). Consider the ODE with (IVP)

y  + 2y/ = 2 - e-4z ,        0 < x < 1,        y(0) = 1.

which has exact solution given by

y(x) = 1 + e-4z - e-2z

(a) Alter the function f and yexact in the code ode_rungekutta.c so that it solves the

ODE. Compile and run the program and choose n = 10 for the number of steps. Write out the table as your solution.

(b) Add an additional function called RK3 which implements the Runge-Kutta method

in Eq. (1). Compile and run the program and choose n = 10 for the number of steps. Write out the table as your solution.

(c)  Advanced  By  using  the  two  Runge-Kutta  methods  together,  implement  a  basic adaptive step size (error correcting) technique.  The function included in the code (stepAdjust) implements the formula

q = ╱ ← 1/2