Scribd
Upload
45 views

Tutorialsheet 4

This document contains 16 numerical analysis problems involving: 1) Using the trapezoidal rule and Simpson's rule to approximate integrals with multiple subintervals. 2) Finding the degree of precision for a given quadrature formula. 3) Solving initial value problems using Euler's method, modified Euler's method, and Runge-Kutta methods to specified accuracies. 4) Approximating solutions to boundary value problems using finite difference methods.

Uploaded by

Gaurav Upadhyay
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
45 views

Tutorialsheet 4

This document contains 16 numerical analysis problems involving: 1) Using the trapezoidal rule and Simpson's rule to approximate integrals with multiple subintervals. 2) Finding the degree of precision for a given quadrature formula. 3) Solving initial value problems using Euler's method, modified Euler's method, and Runge-Kutta methods to specified accuracies. 4) Approximating solutions to boundary value problems using finite difference methods.

Uploaded by

Gaurav Upadhyay
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/ 2

Dr B R AMBEDKAR NIT JALANDHAR

NUMERICAL METHODS (MACI-203)

Tutorial Sheet-4

(1) Use the Trapezoidal rule to approximate the given integral.


Z 2
1
6
dx, n = 6.
0 1+x

(2) Use Trapezoidal rule to compute the given integral using three intervals. Compare it
with exact value
Z 2
1
dx.
1 x

(3) Find the degree of precision of the quadrature formula


Z 1 √ √
3 3
f (x)dx = f (− ) + f( ).
−1 3 3
R1
(4) The quadrature formula −1 f (x)dx = c0 f (−1) + c1 f (0) + c2 f (1) is exact for all the
polynomials of degree less than or equal to 2. Determine c0 , c1 and c2 .

(5) Calculate using Simpson’s 1/3 rule


Z π/2 √
sin xdx, n = 6.
0

R2
(6) Evaluate 0 ex dx using the Simpson’s rule with h = 1 and 1/2. Find a bound on the
error in each case. Compare with the exact solution.

R 0.5 R 0.5 sin xy


(7) Estimate 0 0
dxdy using Simpson’s rule for double integrals with both step
1 + xy
sizes = 0.25.

(8) Evalute the double integral using the Trapezoidal Rule with h = k = 0.25.
Z 1 Z 2 
2xy
2 2
dy dx.
0 1 (1 + x )(1 + y )

(9) Discretize the following initial value problems


(a). y 0 = t + y, y(1) = 0
0 2
(b). y = −y , y(1) = 1
using Euler’s method Compute y(1.2), with h = 0.1 in each case.
(10) Use Taylor’s method of order two to approximate the solutions for each of the following
initial value problems.
(a). y 0 = te3t − 2y, 0 ≤ t ≤ 1, y(0) = 0, with h = 0.5
(b). y 0 = 1 + y/t, 1 ≤ t ≤ 2, y(1) = 2, with h = 0.25
(c). y 0 = 1 + (t − y)2 , 2 ≤ t ≤ 3, y(2) = 1, with h = 0.5.

(11) Use Euler’s method to approximate the solutions for each of the following initial value
problems.
(a). y 0 = cos 2t + sin 3t, 0 ≤ t ≤ 1, y(0) = 1, with h = 0.25
0 2
(b). y = −5y + 5t + 2t, 0 ≤ t ≤ 1, y(0) = 31 , with h = 0.1.

(12) Use the Modified Euler method(Runge-Kutta method of two order) to approximate the
solution to the following initial-value problem, and compare the result to the actual
values.
(a). y 0 = 1+ y/t, 1 ≤ t ≤ 2, y(1) = 2, with h = 0.25; actual solution y(t) = t ln t +2t
1
(b).y 0 = 1 + (t − y)2 , 2 ≤ t ≤ 3, y(2) = 1, h = 0.5; actual solution y(t) = t + .
1−t
(13) Solve the initial value problem
u0 = −2tu2 , y(0) = 1
with h = 0.2 on the interval [0, 0.4]. Use the forth order classical Runge-Kutta method.
Compare with the exact solution.

(14) Use the classical Runge-Kutta method of forth order to find the numerical solution at
x = 0.8 for
dy √
= x + y, y(0.4) = 0.41.
dx
Assume the step length h = 0.2.

(15) Approximate the solutions to the following initial-value problems. Compare the results
to the actual values.
(a). y 0 = 1 + y/t + (y/t)2 , 1 ≤ t ≤ 3, y(1) = 0; actual solution y(t) = t tan(ln t).
(b). y 0 = y/t − (y/t)2 , 1 ≤ t ≤ 4, y(1) = 1; actual solution y(t) = t/(1 + ln t).
(c).y 0 = −(y + 1)(y + 3), 0 ≤ t ≤ 3, y(0) = −2, h = 0.5; actual solution y(t) =
−3 + 2(1 + e−2t )−1 .

(16) Use the linear Finite-Difference method to approximate the solution to the following
boundary-value problems .
(a). y 00 = −3y 0 + 2y + 2x + 3, 0 ≤ x ≤ 1, y(0) = 2, y(1) = 1; use h = 0.1.
(b). y 00 = −(x+1)y 0 +2y+(1−x2 )e−x , 0 ≤ x ≤ 1, y(0) = −1, y(1) = 0; use h = 0.1.

You might also like