Scribd
Upload
19 views1 page

Y y y H F F y X, X,, y .: Actual Solution

The document outlines Problem Set 4 for MTL 712, focusing on computational methods for differential equations. It includes tasks such as checking the convergence of a numerical method, implementing MATLAB scripts for various IVPs using Adam-Bashforth and Adam-Moulton methods, and comparing errors in approximations. Additionally, it involves solving specific IVPs using predictor-corrector methods and determining the order of convergence for the solutions.

Uploaded by

kinjal.iitd
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
19 views1 page

Y y y H F F y X, X,, y .: Actual Solution

The document outlines Problem Set 4 for MTL 712, focusing on computational methods for differential equations. It includes tasks such as checking the convergence of a numerical method, implementing MATLAB scripts for various IVPs using Adam-Bashforth and Adam-Moulton methods, and comparing errors in approximations. Additionally, it involves solving specific IVPs using predictor-corrector methods and determining the order of convergence for the solutions.

Uploaded by

kinjal.iitd
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

MTL 712: Computational Methods for Differential Equations

Problem Set 4

1. Check the convergence of the numerical method

yn+2 + 4yn+1 − 5yn = h(4fn+1 + 2fn )

used to solve the IVP y 0 (x) = 0, x ∈ [0, 5], y(0) = 1.

2. Write MATLAB script to implement the Adam- Bashforth and Adam-Moulton two-step methods
to solve the below given IVPs. Compare the errors in the approximations.

(i) y 0 = y/x − (y/x)2 , 1 ≤ x ≤ 2, y(1) = 1, with step size h = 0.1,


x
actual solution y(x) =
1 + ln(x)
(ii) y 0 = 1 + y/x + (y/x)2 , 1 ≤ x ≤ 3, y(1) = 0, with step size h = 0.2,
actual solution y(x) = x tan(ln x)
(iii) y 0 = −(x + 1)(x + 3), 0 ≤ x ≤ 2, y(0) = −2, with step size h = 0.2
2
actual solution y(x) = −3 +
1 + exp(−2x)
(iv) y 0 = −5y + 5x2 + 2x, 0 ≤ x ≤ 1, y(0) = 1/3, with step size h = 0.1
1
actual solution y(x) = x2 + exp(−5x)
3

3. Redo the Problem 1 using Adam- Bashforth and Adam-Moulton three and four-step methods.
Compare the numerical results of Problem 1 and Problem 2. Find the order of convergence.

4. Solve the IVP


y 0 = ey , x ∈ [0, 0.2], y(0) = 1.
using the predictor corrector method by taking mid-point method as predictor and modified
Euler’s method as corrector. What is the numerically observed order of convergence?

5. Use the predictor-corrector method with 3-step Adam-Bashforth method as predictor and 3-
step Adam-Moulton method as corrector to find the approximation of y(1.1) for the initial value
problem
y 0 = 2xy, y(1) = 1.
taking h = 0.01. Compare the approximate solution with the exact solution to check the accuracy
of the method. Also, find the order of convergence.

You might also like