Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
39 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
Available Formats
Download as PDF, TXT or read online on Scribd
39 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
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