Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: THEend8_
1. The BVP y ′′ ? q(x)y = r(x) with y(a) = y(b) = 0 is solved using the difference scheme
where j = 1, 2, · · · , N, rj = r(xj ), qj = q(xj ), u0 = uN+1 = 0, and h = (b ? a)/(N + 1).
(a) (20 points) Assuming continuous higher order derivatives of y, q, and r, find α?1, α0, and α1, such that the truncation error is O(h 4 ). (Note: y (4) ? [q(x)y] ′′ = r ′′(x))
(b) (10 points) Using your result from (a), give the linear system Ax = b (in component form) that must be solved to get the solution uj for j = 1, 2, · · · , N.
(c) (10 points) Prove: If q(x) ≥ Q > 0 and h is sufficiently small, then
2. Consider the following 2-step method for solving the IVP ˙y = f(y), y(0) = y0.
yn+1 + a0yn + a1yn?1 = hbfn+1
(a) (10 points) How should a0, a1, and b be chosen to ensure the method is consistent?
(b) (10 points) How should a0, a1, and b be chosen to ensure the method is stable?
3. (40 points) The IVP ˙y = f(t, y), y(0) = y0 is integrated to get y(1) = y(0) + f(t, y(t))dt. Show that approximating the integral with the Gaussian quadrature formula leads to the following Runge-Kutta method for solving the IVP.
Remark: The Runge–Kutta methods take the form.
The coefficients of the method are collected together in the Butcher tableau as follows:
Hint: apply the Gaussian quadrature formula on the following integral
where the change of variable t = tn + (tn+1 ? tn)s := tn + hs is considered.