Syllabus Ma 8452 - Statistics and Numerical Methods
Syllabus Ma 8452 - Statistics and Numerical Methods
Syllabus Ma 8452 - Statistics and Numerical Methods
UNIT - I
TESTING OF HYPOTHESIS
UNIT - II
DESIGN OF EXPERIMENTS
One way and two way classifications - Completely randomized design – Randomized block
design – Latin square design - 22 factorial design.
UNIT - III
SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS
Solution of algebraic and transcendental equations - Fixed point iteration method – Newton
Raphson method - Solution of linear system of equations - Gauss elimination method – Pivoting
- Gauss Jordan method – Iterative methods of Gauss Jacobi and Gauss Seidel - Eigenvalues of a
matrix by Power method and Jacobi’s method for symmetric matrices.
UNIT - IV
INTERPOLATION, NUMERICAL DIFFERENTIATION AND NUMERICAL
INTEGRATION
Lagrange’s and Newton’s divided difference interpolations – Newton’s forward and backward
difference interpolation – Approximation of derivates using interpolation polynomials –
Numerical single and double integrations using Trapezoidal and Simpson’s 1/3 rules.
UNIT - V
NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
Single step methods : Taylor’s series method - Euler’s method - Modified Euler’s method -
Fourth order Runge-Kutta method for solving first order equations - Multi step methods :
Milne’s and Adams - Bash forth predictor corrector methods for solving first order equations.
MA6452 STATISTICS AND NUMERICAL METHODS
QUESTION BANK
UNIT I
TESTING OF HYPOTHESIS
PART A
Solution:
On the basis of sample statistics, the study of difference between the observed sample
statistics and the hypothetical population parameter value is significant or not is defined as
test of significance.
Solution:
The statistical constants like mean µ , variance σ 2 computed from population are called
parameters where as statistical constants like mean x̄, variance s2 etc., computed from a
sample are called sample statistics.
Solution:
The probability of type - II error is denoted by β and (1 − β) is defined as power of the
statistical test.
Solution:
The probability of Type I error is called the level of significance of the test and is denoted
by α.
M A6452 QUESTION BANK
Solution:
The rejection of H0 when it is true.
Solution:
The acceptance of H0 when it is false.
Solution:
For a test statistic, the area under the probability curve, which is normal, is divided into
two regions namely the region of acceptance of H0 and the region of rejection of H0 . The
region in which H0 is rejected is called the critical region. The area of the critical region is
α, the level of significance.
Solution:
In making statistical decisions, we make assumptions or guesses about populations involved.
Such assumptions, which may be true or false are called statistical hypothesis.
Solution:
For applying the test of significance, we first set up a hypothesis which is a statement about
the population parameter. This statement is usually a hypothesis of no difference and so it
is called null hypothesis and is denoted by H0 . The purpose of null hypothesis is for possible
rejection under the assumption that it is true.
Solution:
Suppose the null hypothesis is false, then something else must be true. This is called an
alternate hypothesis and is denoted by H1 .
KCG\S&H 2
M A6452 QUESTION BANK
Solution:
One tailed test: In a test of any statistical hypothesis, if the alternate hypothesis is
one-sided then it is called one tailed test. It may be right tailed or left tailed.
Two tailed test: In a test of statistical hypothesis, if the alternate hypothesis is two sided,
then it is called a two tailed test.
Solution:
The number of elements in a sample is greater than or equal to 30, then the sample is called
a large sample and if it is less than 30, then the sample is called a small sample.
Solution:
Testing the significance of the difference between
(i) The mean of a sample and the mean of the population.
(ii) The means of two samples.
Solution:
(i) Testing the significance of the difference between the variances of two populations from
which two samples are drawn.
(ii) analysis of variance.
Solution:
KCG\S&H 3
M A6452 QUESTION BANK
16. Write conditions for the application of Chi-square test. [AU 2014]
Solution:
1. The experimental data (or sample deviations) must be independent of each other.
4. The Pconstraints
P on the cell frequencies should be linear.
eg., Oi = Ei = N ≥ 50.
Solution:
An attribute is a characteristic or a quality which may be present amongst the members of
a population.
For example, tall, short, healthy, black etc.
Solution:
Let A and B two attributes. Dividing A into A1 and A2 and B into B1 , B2 , we get the
following 2 × 2 table, called the 2 × 2 table.
A Total
B
( ) ) )
) ) )
Total ) ) N
KCG\S&H 4
M A6452 QUESTION BANK
a b
19. Give the formula for the χ2 -test of independence for [AU 2012]
c d
Solution:
A Total
B
Total
N (ad − bc)2
∴ the value of χ2 =
(a + b)(c + d)(a + c)(b + d)
20. A random sample of 200 tins of coconut oil gave an average weight of 4.95 kgs.
With a standard deviation of 0.21 kg. Do we accept the net weight is 5 kgs per tin at 5% level?
[AU 2013]
Solution:
Let H0 = µ = 5
H1 = µ 6= 5
x̄ − µ 4.95 − 5 −0.05 √
∴ Z= √ = √ = 10 2 = −3.37
s/ n 0.21/ 200 0.21
KCG\S&H 5
M A6452 QUESTION BANK
21. Mention the various steps involved in testing of hypothesis. [AU 2010]
Solution:
t − E(t)
Step 4. Compute the test statistic Z =
S.Eof (t)
Step 5. Compare the computed value of |Z| with the table value of Z and decide
the acceptance or the rejection of H0 .
Step 6. Inference.
KCG\S&H 6
M A6452 QUESTION BANK
PART B
1. A sample of 100 people during the past year showed an average life span of 71.8 years. If the
standard deviation of the population is 8.9 years, test whether the mean life span today is
greater than 70 years. [AU 2010]
2. A sample of 900 items has mean 3.4 cms and standard deviation 2.61 cms. Can the sample be
regarded as drawn from a population with mean 3.25 cms at 5% level of significance.[AU 2008, 2010]
3. A sample of heights of 6400 Englishmen has a mean of 170 cms and a standard deviation of
6.4 cms, while a sample of heights of 1600 Australians has a mean of 172 cm and standard
deviation of 6.3 cm. Do the data indicate that the Australians are on the average taller than
the Englishmen. [AU 2007]
4. A machine produces 16 imperfect articles in a sample of 500. After it was overhauled, it puts ou
3 imperfect articles in a sample of 100. Has the machine improved in its performance? [AU 2012]
5. A mathematics test was given to 50 girls and 75 boys. The girls made an average grade of 76
with a SD of 6, while boys made an average grade of 82 with a SD of 2. Test whether there is
any significant difference between the performance of boys and girls. [AU 2012]
6. Test whether there is any significant difference between the variances of the populations from
which the following samples are taken: [AU 2012]
Sample I: 20 16 26 27 23 22
.
Sample II: 27 33 42 35 32 34 38
7. A sample of 10 boys had the I.Q’s:70, 120, 110, 101, 88, 83, 95, 98, 100 and 107. Test whether
the population mean I.Q may be 100. [AU 2012]
8. Two independent sample of sizes 9 and 7 from a normal population had the following values of
the variables. [AU 2014]
Sample I: 18 13 12 15 12 14 16 14 15
.
Sample II: 16 19 13 16 18 13 15
9. Out of 8000 graduates in a town 800 are females, out of 1600 graduate employees 120 are
females. Use χ2 to determine if any distinction is made in appointment on the basis of sex.
Value of χ2 at 5% level for one degree of freedom is 3.84. [AU 2010]
10. An automobile company gives you the following information about age groups and the liking
for particular model of car which it plans to introduce. On the basis of this data can it be
concluded that the model appeal is independent of the age group. (χ20.05 (3) = 7.815) [AU 2010]
Age Group
Persons who : Below 20 20 − 39 40 − 59 60 and above
Liked the car : 140 80 40 20 .
disliked the car : 60 50 30 80
KCG\S&H 7
M A6452 QUESTION BANK
UNIT II
DESIGN OF EXPERIMENTS
PART A
Solution:
The three basic principles of experimental design are 1. Randomisation, 2. Replication and
3. Local control.
Solution:
The prime objective of design of experiments is to control the extraneous variables so that
the results could be attributed only to the experimental variables.
Solution:
It is defined as the logical construction of the experiment in which the degree of uncertainty
with which the inference is drawn, may be well defined.
Solution:
1. The samples are drawn from normal populations.
2. The samples are drawn independently from these populations.
3. All the populations have the same variance.
ANOVA should not be used if we cannot make these assumptions.
Solution:
When we have to test the differences between means of more than two samples we use
analysis of variance.
KCG\S&H 8
M A6452 QUESTION BANK
6. Write down the ANOVA table for one way classification. [2013]
Solution:
One way classification ANOVA table
Source of Sum of Degrees of Mean square Variance
Variation squares (SS) freedom (df ) (Ms) ratio
Between SSB M SB
SSB r−1 MSB= F=
samples r−1 M SW
SSW M SW
Within samples SSW N −r MSW= (or)
N −r M SB
(or error)
Total SST N −1
7. Compare one-way classification model with two-way classification model. [AU 2010]
Solution:
1. CRD analysis results in one-way classification, whereas RBD analysis results in two-way
classification.
2. Experimental errors are large in CRD compared to RBD and RBD is popular.
8. State the null and alternate hypothesis for a CRD. [AU 2009]
Solution:
H0 : µ1 = µ2 = · · · = µk = µ i.e., the treatment means are equal.
H1 : not all µ′i s are equal.
Solution:
KCG\S&H 9
M A6452 QUESTION BANK
10. Discuss the advantages and disadvantages of Randomized block design. [2010]
Solution:
Advantages:
1. It has a simple layout but it is more efficient than CRD because of reduction of
experimental error.
2. It is flexible and so any number of treatments and any number of replication may be
used.
Disadvantages:
1. If the number of treatments is large, then the size of the blocks will increase this may
cause heterogeneity within blocks.
Solution:
Suppose we want to test the effect of r fertilizers on the yield of paddy. We divide the plots
into h blocks, each block is relatively homogeneous and each block contains r plots. Within
each block the plots are selected at random and the r treatments (i.e., fertilizers) are given.
Thus in each block only one plot receives one fertiliser. This repeated for all the h blocks.
This design is called randomized block design.
12. Write any two differences between RBD and CRD. [2011]
Solution:
1. CRD analysis results in one-way classification, whereas RBD analysis results in two-way
classification.
2. Experimental errors are large in CRD compared to RBD and RBD is popular.
Solution:
Latin square design controls variation in two directions of the experimental material as rows
and columns resulting in the reduction of experimental error.
KCG\S&H 10
M A6452 QUESTION BANK
Solution:
15. Explain the situations in which randomized block design is considered an improvement over a
completely randomized design. [AU 2014]
Solution:
• RBD has a simple layout but it is more efficient than CRD because of reduction of
experimental error.
Solution:
If the variation due to heterogeneity in experimental units is so large then the sensitivity of
detecting treatment differences is reduced because of large value of s2 . A better idea would
be to "block off" variation due to these units and thus reduce the extraneous variation by
smaller homogeneous blocks.
17. State the advantages of a factorial Experiment over a simple experiment. [2010,2014]
Solution:
Factorial design is one of the fruitful advancement in the endevour to improve the logical
foundations of experimental designs. In experiments based on factorial design, the experiment
can evaluate the combined effect of two or more factors when used simultaneously.
Solution:
When there are two factors A, B and two levels ’high’ and ’low’ for each factor we have
a 22 factorial design. Inspite of its simplicity, the 22 design is a powerful tool to improve
products and processes.
KCG\S&H 11
M A6452 QUESTION BANK
Solution:
corresponding SS
Mean square=
df
Solution:
Consider a n × n Latin square design, then the degrees of freedom for SSE is
= (n2 − 1) − (n − 1) − (n − 1) − (n − 1)
= n2 − 1 − 3n + 3
= n2 − 3n + 2
= (n − 1)(n − 2)
For n = 2, d.f of SSE= 0 and hence MSE is not defined.
∴ comparisons are not possible. Hence a 2 × 2 Latin Square Design is not possible.
KCG\S&H 12
M A6452 QUESTION BANK
PART B
1. A set of data involving four "four tropical feed stuffs A, B, C, D" tried on 20 chicks is given
below. All the twenty chicks are treated alike in all respects except the feeding treatments and
each feeding treatment is given to 5 chicks. Analyze the data [AU 2010]
Weight gain baby chicks fed on different feeding materials composed of tropical feed stuffs:
TOTAL Ti
A 55 49 42 21 52 219
B 61 112 30 89 63 355
.
C 42 97 81 95 92 407
D 169 137 169 85 154 714
Total 1695
2. An experiment was planned to study the effect of sulphate of potash and super phosphate on the
yield of potatoes. All the combinations of 2 levels of super phosphate and 2 levels of sulphate
of potash were studied in a randomized block design with 4 replications for each. The yields
(per plot) obtained are given below. [AU 2010]
Analyze the data and comment on your findings. (F0.05 (3, 9) = 3.86, F0.05 (1, 9) = 5.12).
A B C D
1 44 38 47 36
2 46 40 52 43
Workers .
3 34 36 44 32
4 43 38 46 33
5 38 42 49 39
KCG\S&H 13
M A6452 QUESTION BANK
I II III
1 0.194 0.73 1.187
A B C
2 0.758 0.311 0.589 .
C A B
3 0.369 0.558 0.311
B C A
5. The following are the number of mistakes made in 5 successive days by 4 technicians working
for a photographic laboratory test at a level of significance α = 0.01. Test whether the difference
among the four sample means can be attributed to chance.[2011]
Technician
I II III IV
6 14 10 9
14 9 12 12 .
10 12 7 8
8 10 15 10
6. The following data represent the number of units of production per day turned out by different
workers using 4 different types of machines[AU 2011, 2013]
Machine type
A B C D
1 44 38 47 36
2 46 40 52 43
Workers .
3 34 36 44 32
4 43 38 46 33
5 38 42 49 39
1. Test whether the five men differ with respect to mean productivity and
2. Test whether the mean productivity is the same for the four different machine types.
7. (a) What are the basic assumptions involved in ANOVA?[2011] [4]
(b) In a Latin square experiment given below are the yields in quintals per acre on the paddy [12]
crop carried out for testing the effect of five fertilizers A, B, C, D, E. Analyze the data for
variations.[AU 2011]
8. The sales of 4 salesmen in 3 seasons are tabulated here. Carry out an analysis of variance.[AU
2012]
Salesmen
Seasons A B C D
Summer 36 36 21 35
.
Winter 28 29 31 32
Monsoon 26 28 29 29
9. A farmer wishes to test the effect of 4 fertilizers A, B, C, D on the yield of wheat. The fertilizers
are used in a LSD and the result are tabulated here. Perform an analysis of variance.[AU 2012]
A 18 C 21 D 25 B 11
D 22 B 12 A 15 C 19
.
B 15 A 20 C 23 D 24
C 22 D 21 B 10 A 17
KCG\S&H 14
M A6452 QUESTION BANK
10. Four varities A, B, C, D of a fertilizer are tested in RBD with 4 replications. The plot yields in
pounds are as follows:[AU 2012,2013,2014]
A 12 D 20 C 16 B 10
D 18 A 14 B 11 C 14
.
B 12 C 15 D 19 A 13
C 16 B 11 A 15 D 20
11. A variable trial was conducted on wheat with 4 varities in a Latin Square design. The plan of
the experiment and per plot yield are given below:[AU 2012]
C 25 B 23 A 20 D 20
A 19 D 19 C 21 B 18
.
B 19 A 14 D 17 C 20
D 17 C 20 B 21 A 15
12. The following is a Latin square of a design when 4 varieties of seeds are being tested. Set up
the analysis of variance table and state your conclusion. You may carry out suitable change of
origin and scale.[AU 2013]
13. Compare and contrast the Latin square Design with the Randomized Block Design. [AU 2013]
14. Analyse the variance in the Latin square of yields in kgs of paddy where P, Q, R, S denote the
different methods of cultivation. [AU 2014]
15. A company wants to produce cars for its own use. It has to select the make of the car out of
the four makes A, B, C and D available in the market. For this he tries four cars of each make
by assigning the cars to four drivers to run on four different routes. The efficiency of cars is
measured in terms of time in hours. The layout and time consumed is as given below.
Drivers
Routes 1 2 3 4
1 C 18 D 12 A 16 B 20
2 D 26 A 34 B 25 C 31 .
3 B 15 C 22 D 10 A 28
4 A 30 B 20 C 15 D9
Analyze the experimental data and draw conclusions. (F0.05 (3, 5) = 5.41)
KCG\S&H 15
M A6452 QUESTION BANK
UNIT III
SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS
PART A
Solution:
The order of convergence is 2.
Solution:
The condition is |f (x).f ′ (x)| < (f ′ (x))2 in a neighbourhood of the root.
Solution:
1. It can be used for finding root of both algebraic and transcendental equations.
3. It is simple and easy to deal with and it is used to improve the results obtained by
other methods.
Solution:
The augmented matrix is
1 1 : 2
[A, B] =
2 3 : 5
1 1 : 2
∼ R2 → R2 − 2R1
0 1 : 1
∴ y = 1, x + y = 2 ⇒ x = 1
∴ x = 1, y = 1
KCG\S&H 16
M A6452 QUESTION BANK
Solution:
Gauss-elimination method and Gauss-Jordan method are the direct methods.
Solution:
The process of iteration will converge if in each equation one coefficient is much larger than
the other two and the largest coefficient must be attached to different coefficients in different
equations.
Solution:
In Jacobi’s method at each iteration the values of the variables in the previous iteration
are used, where as in Gauss-Seidel method in each iteration the latest available values of
the variables are used. Hence the convergence of Gauss-Seidel method is twice faster than
Jacobi’s method.
1 3
8. Find the inverse of A = by Gauss-Jordan method.
2 7
Solution:
The augmented
matrix is
1 3 : 1 0 1 3 : 1 0
[A, B] = ∼ R2 → R2 + (−2)R1
2 7 : 0 1 0 1 : −2 1
1 0 : 7 −3
∼ R1 → R1 + (−3)R2
0 1 : −2 1
−1 7 −3
A = R1 → R1 + (−3)R2
−2 1
9. Write down all possible types of initial vectors to determine the largest eigen value and the
corresponding eigen vector of a matrix size 2 × 2
Solution:
1 1 0
, ,
1 0 1
KCG\S&H 17
M A6452 QUESTION BANK
Solution:
Sl. COMPARISON
No. Gauss Jordan Gauss Jacobi
1 Direct method Itterative method
2 This method produce exact solution This method give a sequence
after a finite number of steps of approximate solutions, which
ultimately approach the actual
solution.
3 Applicable if the coefficient matrix Applicable if the coefficient matrix
is non-singular is diagonally dominant.
Solution:
Sl. COMPARISON
No. Gauss elimination Gauss Jacobi
1 Direct method Itterative method
2 This method produce exact solution This method give a sequence
after a finite number of steps of approximate solutions, which
ultimately approach the actual
solution.
3 Applicable if the coefficient matrix Applicable if the coefficient matrix
is non-singular is diagonally dominant.
12. What is the order of convergence and also state the error term for Newton Raphson method?
[AU 2011]
Solution:
f ”(x) 2
Order of convergence is 2 and the error term is ǫk+1 ≈ ǫ
2f ′ (x) k
1 2
13. Find the dominant eigen value of the matrix by power method. [AU 2011]
3 4
Solution:
1 2 0
Given A = . Take X0 = as the initial vector (II row contains the largest
3 4 1
element and largest
sum)
1 2 0 2 0.5
Then AX0 = = =4 = 4X1
3 4 1 4 1
1 2 0.5 2.5 0.4545
AX1 = = = 5.5 = 5.5X2
3 4 1 5.5 1
KCG\S&H 18
M A6452 QUESTION BANK
1 2 0.4545 2.4545 0.4576
AX2 = = = 5.3635 = 5.3635X3
3 4 1 5.3635 1
1 2 0.4576 2.4576 0.4574
AX3 = = = 5.3728 = 5.3728X4
3 4 1 5.3728 1
1 2 0.4574 2.4574 0.4574
AX4 = = = 5.3722 = 5.3722X5
3 4 1 5.3722 1
Since X4 = X5 , we stop the iteration and the largest eigen value is 5.3722
√
14. Using Newton-Raphson method, find the iteration formula to compute N .[AU 2010]
Solution:
√
Let x = N , then x2 = N ⇒ x2 − N = 0
Let f (x) = x2 − N, then f ′ (x) = 2x
f (xn )
Newton’s iteration formula is xn+1 = xn − , n = 0, 1, 2, · · ·
f ′ (xn )
x2n − N
∴ xn+1 = xn −
2xn
2x2n − x2n + N
=
2xn
x2n + N
=
2xn
1 N
= xn +
2 xn
15. Explain the power method to determine the eigenvalue of a matrix.[AU 2010]
Solution:
Workingrule
1. Let x0 be the initial which is usually chosen as a vector with all components equal to
1. (i.e., normalised)
2. Form the product AX0 and express it in the form AX0 = λ1 X1 , where X1 is normalized
by taking out the largest component λ1 .
KCG\S&H 19
M A6452 QUESTION BANK
√
3
16. Arrive a formula to find the value of N , where N 6= 0, using Newton-Raphson method.[AU
2011]
Solution:
√
Let x = 3 N , then x3 = N ⇒ x3 − N = 0
Let f (x) = x3 − N, then f ′ (x) = 3x2
f (xn )
Newton’s iteration formula is xn+1 = xn − , n = 0, 1, 2, · · ·
f ′ (xn )
x3n − N
∴ xn+1 = xn −
3x2n
3x3n − x3n + N
=
3x2n
2x3n + N
=
3x2n
1 N
= 2xn + 2
3 xn
17. Solve the following system of equations, using Gauss-Jordan elimination method 2x + y =
3, x − 2y = −1. [AU 2012]
Solution:
Given 2x + y = 3, x − 2y = −1.
Augmented matrix
2 1 3
[A, B] =
1 −2 −1
1 1/2 3/2 R1
∼ R1 −→
1 −2 −1 2
1 1/2 3/2
∼ R2 −→ R2 − R1
0 −5/2 −5/2
1 1/2 3/2 R2
∼ R1 −→
0 1 1 −5/2
1 0 1
∼ R1 −→ R1 − (1/2)R2
0 1 1
∴ the solution is x = 1, y = 1
KCG\S&H 20
M A6452 QUESTION BANK
18. Form the divided difference table for the following data:[AU 2012]
x: 5 15 22
y: 7 36 160
Solution:
The Divided difference table
x y ∆y ∆2 y
5 7
36−7
15−5 = 2.6
17.71−2.6
15 36 22−5 = 0.88
160−36
22−15 = 17.71
22 160
19. Find the real positive root of 3x − cos x − 1 = 0 by Newton’s method correct to 6 decimal
places.[2013]
Solution:
Given eqn. is 3x − cos x − 1 = 0.
Let f (x) = 3x − cos x − 1 then f ′ (x) = 3 + sin x
f (0) = −1 < 0,
f (1) = 3 − cos 1 − 1 = 1.4597 > 0
∴ a root lies between 0 and 1.
Now f (0.5) = 3(0.5) − cos 0.5 − 1 = −0.3776, which is closer to 0.
So, the root is indeed between 0.5 and 1 and is nearer to 0.5.
f (xn )
Take x0 = 0.6 Iteration formula is xn+1 = xn − ′ , n = 0, 1, 2, · · ·
f (xn )
I approximation:
f (x0 )
x1 = x0 −
f ′ (x0 )
f (0.6)
= 0.6 − ′
f (0.6)
−0.0253356
= 0.6 −
3.564642
= 0.6 + 0.0071
= 0.6071
II approximation:
f (x1 )
x2 = x1 −
f ′ (x1 )
f (0.6071)
= 0.6071 − ′
f (0.6071)
−0.000006
= 0.6 −
3.570488
= 0.6 + 0.00000168
= 0.60710168
∵ x1 = x2 upto 4 decimal places, the root is 0.6071
KCG\S&H 21
M A6452 QUESTION BANK
Solution:
KCG\S&H 22
M A6452 QUESTION BANK
PART B
Newton-Raphson method
2. Find by Newton Raphson method a positive root of the equation 3x−cosx−1 = 0. [AU 2006, 2009]
3. Using Newton’s method, find a real root of xlog10 x = 1.2 correct to 4 decimals.[AU 2004, 2007]
4. Use Newton’s Raphson method to find a root of the equation cosx − xex = 0 correct to 3
decimal places.[AU 2000, 2006]
9. Find the iterative formula for finding the value of N1 where N is a real number, using
1
Newton-Raphson method. Hence evaluate 26 correct to 4 decimal places. [AU 2006]
1
10. Find the Newton’s iterative formula for the reciprocal of N and hence find the value of 23 ,
correct to 5 decimal places. [AU 2006, 2012]
√
11. Find an iterative
√ formula for 3 N , where N is a positive integer, using Newton’s method and
√
hence find 3 24, 3 41
1 √1 ,
12. Find the iterative formula by Newton’s formula for (i) N (ii) N
where N is a positive integer.
1 √1
Hence find 31 and 15
13. Solve for a positive root of the equation x4 −x−10 = 0 using Newton’s Raphson’s method.[AU 2010]
14. Find a root of 3x3 − 9x2 + 8 = 0 using Newton Raphson method. [AU 2003]
15. Find a root of x3 − 6x + 4 = 0 between 0 and 1 to 5 places of decimals. [AU 2000, 2008]
16. Find a root of 2x3 − 3x − 6 = 0 between 1 and 2 to five places of decimals. [AU 2008]
17. Using Newton’s method, solve x log10 x = 12.34 taking the initial value x0 as 10 . [AU 2004, 2012, 2014]
KCG\S&H 23
M A6452 QUESTION BANK
KCG\S&H 24
M A6452 QUESTION BANK
5. Solve by Gauss - Seidel method the equations starting with (0,0,0,0) as solution. Do 5
iterations only.
4x1 − x2 − x3 = 2
−x1 + 4x2 − x4 = 2
−x1 + 4x3 − x4 = 2
−x2 − x3 + 4x4 = 1 [AU 2008]
KCG\S&H 25
M A6452 QUESTION BANK
2 6 6
4. Using Gauss Jordan method, find the inverse of the matrix A= 2 8 6 [AU 2006]
2 6 8
2 2 3
5. Using Gauss Jordan method, find the inverse of the matrix A= 2 1 1 [AU 2004]
1 3 5
0 1 2
6. Using Gauss Jordan method, find the inverse of the matrix A= 1 2 3 [AU 2010]
3 1 1
0 1 1
7. Using Gauss Jordan method, find the inverse of the matrix A= 1 2 0 [AU 2010]
3 −1 −4
1 1 1
8. Using Gauss Jordan method, find the inverse of the matrix A= 4 3 1 [AU 2011]
3 5 3
KCG\S&H 26
M A6452 QUESTION BANK
UNIT IV
INTERPOLATION, NUMERICAL DIFFERENTIATION AND NUMERICAL
INTEGRATION
PART A
Solution:
(x−x1 )(x−x2 )···(x−xn )
y = f (x) = (x0 −x1 )(x0 −x2 )···(x0 −xn ) y0
+ (x(x−x
1 −x
0 )(x−x2 )···(x−xn )
0 )(x1 −x2 )···(x1 −xn )
y1
2. Use Lagrange’s formula to fit a polynomial to the data and find y at x = 1. [2013]
X −1 0 2 3
Y −8 3 1 12
Solution:
Given x0 = −1 x1 = 0 x2 = 2 x3 = 3
and y0 = −8 y1 = 3 y2 = 1 y3 = 12
(1−0)(1−2)(1−3) (1+1)(1−2)(1−3) (1+1)(1−0)(1−3)
y = f (x) = (−1−0)(−1−2)(−1−3) (−8) + (0+1)(0−2)(0−3) (3) + (2+1)(2−0)(2−3) (1)
+ (1+1)(1−0)(1−2)
(3+1)(3−0)(3−2) (12)
(1)(−1)(−2) ✟2
✟ ✟
(2)(✟
−3)(−2) (2)(1)(−2) (2)(1)(−1)
✚
= (−1)(−3)(✟✟✟ (−8) + ✟ (3)
(1)(−2)(✟ + (3)(2)(−1) (1) + ✘✘ ✘ (✚
12)
−4) −3) ✘
(4)(3)(1)
2 −4 −4 8 + 36 + 4 − 12 36
= 2+ 3+ −2= = =6
3 −2 −6 6 6
Solution:
x
Rn h
f (x)dx = 3 [(y0 + yn ) + 4(y1 + y3 + y5 + . . . + yn−1 ) + 2(y2 + y4 + . . . + yn−2 )]
x0
xn −x0
where h = n and n = even no. of sub intervals.
KCG\S&H 27
M A6452 QUESTION BANK
R1 dx
4. Using Trapezoidal rule, evaluate with h = 0.2. Hence obtain an approximate value of
0 1 + x2
π. [AU 2011]
Solution:
Given h=0.2
R1 h
∴ ydx = 2 [(y0 + y5 ) + 2(y1 + y2 + y3 + y4 )]
0
= 0.2
2 [(1 + 0.5) + 2(0.96 + 0.86 + 0.74 + 0.61 + 0.50)]
= 0.1 [(1.5) + 2(3.67)]
= 0.1 [(1.5) + 7.34]
= 0.884
5. State the formula to find the second order derivative using the forward differences. [AU 2011]
Solution:
h i
1 2u−1 3u2 −6u+2 4u3 −18u2 +22u−6
∆2 y ∆3 y ∆4 y
1. y ′ (x)∆y0 +
= h 0+ 2 0+ 6 24 0 + ...
h 2 i
′′ 1 2 3 6u −18u+11 4
2. y (x) = h2 ∆ y0 + (u − 1) ∆ y0 + 12 ∆ y0 + . . .
R1 dx
6. Evaluate x by Trapezoidal rule, dividing the range into equal 4 parts.[AU 2012]
0.5
Solution:
1 − 0.5
Here n=4 ∴ h = = 0.125
4
x 0.5 0.625 0.75 0.875 1
1
y= 2 1.6 1.33 1.14 1
x
y0 y1 y2 y3 y4
Trapezoidal rule is
R1
f (x)dx = h2 [(y0 + y4 ) + 2(y1 + y2 + y3 )]
0.5
0.125
= 2 [(2 + 1) + 2(1.6 + 1.33 + 1.14)]
0.125
= 2 [(3) + 2(4.08)]
= 0.697
KCG\S&H 28
M A6452 QUESTION BANK
7. Find the area under the curve passing through the points (0,0),(1,2),(2,2.5),(3,2.3),(4,2),(5,1.7)
and (6,1.5).
Solution:
Given
x 0 1 2 3 4 5 6
y 0 2 2.5 2.3 2 1.7 1.5
y0 y1 y2 y3 y4 y5 y6
6−0
Here n=6 ∴ h = =1
6
Trapezoidal rule is
R6
ydx = h2 [(y0 + y6 ) + 2(y1 + y2 + y3 + y4 + y5 )]
0
1
= 2 [(0 + 1.5) + 2(2 + 2.5 + 2.3 + 2 + 1.7 + 1.5)]
1
= 2 [(1.5) + 2(12)]
= 12.75
Solution:
1. The divided differences are symmetrical in all their arguments, that is, the value of
any difference is independent of the order of the arguments.
2. The divided difference of the product of a constant and a functions is equal to the
product of the constant and the divided difference of the function.
R2 dx
9. Evaluate 1+x2
by Trapezoidal rule, taking h = 0.5.[AU 2013]
1
Solution:
Given h = 0.5
x 1 1.5 2
1
y= 0.5 0.31 0.2
1 + x2
y0 y1 y2
R2 h 0.5 0.5
∴ ydx = 2 [(y0 + y2 ) + 2(y1 )] = 2 [(0.5 + 0.2) + 2(0.31)] = 2 [(0.7) + 0.62] = 0.33
1
KCG\S&H 29
M A6452 QUESTION BANK
Solution:
It is the process of finding a value of x for the corresponding value of y and we use Lagrange’s
interpolation formula by taking the independent variable as y and the dependent variable
as x. It is the inverse process of direct interpolation in which we find the values of y
corresponding to a value of x, not present in the table.
Solution:
The process of computing intermediate values of a function from a given set of tabular
values of the function.
Newton’s forward and backward interpolation formula can be used only when the values
of the independent variable x are equally spaced and also be when the differences of the
dependent variable y become smaller finally. But Lagrange’s interpolation formula can be
used whether the values of x, the independent variable, are equally spaced or not,whether
the differences of y becomes smaller or not.
12. Find the parabola of the form y = ax2 + bx + c passing through the points (0, 0), (1, 1), (2, 20).
[AU 2011]
Solution:
Given
x0 x1 x2
x 0 1 2
y 0 1 20
y0 y1 y2
(x − x1 )(x − x2 )
∴ y = f (x) = y0
(x0 − x1 )(x0 − x2 )
(x − x0 )(x − x2 )
+ y1
(x1 − x0 )(x1 − x2 )
(x − x0 )(x − x1 )
+ y2
(x2 − x0 )(x2 − x1 )
x2 −2x 2
x −x ✚10
=0+ + (2)(1)
✚
20 = −x2 + 2x + 10x2 − 10x = 9x2 − 7x
✚
(1)(−1) ✚
R1
13. Write down trapezoidal rule to evaluate f (x)dx with h = 0.5, function f (x) is unknown.[AU
6
2005]
KCG\S&H 30
M A6452 QUESTION BANK
Solution:
6−1
Here h=0.5 ∴ n = = 10 x0 = 1
0.5
x 1 1.15 2 2.5 3 3.5 4 4.5 5 5.5 6
f (x) f (1) f (1.5) f (2) f (2.5) f (3) f (3.5) f (4) f (4.5) f (5) f (5.5) f (6)
y0 y1 y2 y3 y4 y5 y6 y7 y8 y9 y10
Trapezoidal rule is
R6
f (x)dx = h2 [(y0 + y10 ) + 2(y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9 )]
1
0.5
= 2 [(f (1) + f (6)) + 2(f (1.5) + f (2) + f (2.5) + f (3) + f (3.5)
+f (4) + f (4.5) + f (5) + f (5.5))]
14. What are the errors in trapezoidal and simpson’s rule of numerical integration.[AU 2005]
Solution:
Rb (b − a)2 2
When evaluating f (x)dx, the error in the trapezoidal rule is < h M, where
a 12
b−a
h= , n is the number of subintervals of (a, b), and
n ′′ |}, y ′′ = f ′′ (x )
M = max{|y0′′ |, |y1′′ |, · · · , |yn−1 r r
(b − a) 4 b−a
The error in Simpson’s rule is < h M, where h = , 2n is the number of
180 2n
subintervals of (a, b), and
(1) (4) (4) (4)
M = max{|y0 |, |y1 |, · · · , |yn−1 |}, yr = f (4) (xr )
1
15. What is the order of error in Simpson’s rule?
3
Solution:
Total error in Simpson’s forumula is of order h4 where h is the length of the subintervals.
Solution:
−h2 ′′
The local error term in the interval (x0 , x1 ) is y
12 0
KCG\S&H 31
M A6452 QUESTION BANK
Solution:
Rb
If the integrand f(x) in f (x)dx is a second degree polynomial, then Simpson’s formula will
a
give exact result.
Solution:
y = f (xn + hv)
Solution:
y = f (x0 + hu)
Rb
20. In order to evaluate f (x)dx by Trapezoidal and simpson’s rule. What is the restriction on
a
the number of intervals?
Solution:
In trapezoidal rule number of intervals can be any positive integer, where as in Simpson’s
rule the number of intervals is even.
KCG\S&H 32
M A6452 QUESTION BANK
PART B
1. Using Newton’s forward interpolation formula find the cubic polynomial which takes the
following values:
x 0 1 2 3
[AU 2000, 2009]
y 1 2 1 10
2. A third degree polynomial passes through the points (0, −1), (1, 1), (2, 1) and (3, −2). Using
Newton’s forward formula, find the polynomial. Hence find the value at 15. [AU 2000].
x 4 6 8 10
y 1 3 8 16
4. Find a polynomial of degree two for the data by Newton’s forward difference method.
x 0 1 2 3 4 5 6 7
[AU 2007].
y 1 2 4 7 11 16 22 29
5. The population of a city in a Census taken once in 10 years is given below. Estimate the
population in the year 1955.
7. From the data given below find the number of students whose weight is between 60 to 70.
9. From the following table of half yearly premium for policies maturing at different ages estimate
the premium for policies maturing at age 46 and 63. [AU 2006, 2008, 2011].
KCG\S&H 33
M A6452 QUESTION BANK
Age x 45 50 55 60 65
Premium y 114.84 96.16 83.32 74.48 68.48
1. Find f (x) as a polynomial in x from the given data and find f(8).
x 3 7 9 10
[AU 2008].
f(x) 168 120 72 63
2. Find the polynomial f (x) by using Lagrange’s formula and hence f (3) for the following values
of x and y.
x 0 1 2 5
y 2 3 12 147
3. Using Lagrange’s interpolation formula find y(10) from the following table:
x 5 6 9 11
[AU 2008, 2011].
y 12 13 14 16
4. Using Lagrange’s interpolation formula to fit a polynomial to the following data hence find
y(1).
x −1 0 2 3
[AU 2010].
y -8 3 1 12
x 5 7 11 13 17
y 150 392 1452 2366 5202
6. Use Lagrange’s method to find log10 656, given that log10 654 = 2.8156, log10 658 = 2.8182,
log10 659 = 2.8189 and log10 661 = 2.8202 [AU 2012].
x 4 5 7 10 11 13
f(x) 48 100 294 900 1210 2028
Also find f (6), f (9), f (15) [AU 2005, 2006, 2011, 2013].
2. Use Newton’s divided difference formula to calculate f (8), f ′ (3) f ′′ (3) from the following table
[AU 2010].
x 0 1 2 4 5 6
f(x) 1 14 15 5 6 19
KCG\S&H 34
M A6452 QUESTION BANK
3. If f (0) = 0, f (1) = 0, f (2) = −12, f (4) = 0, f (5) = 600, f (7) = 7308, find a polynomial that
satisfies this data using Newton’s divided difference formula. Hence find f (6). [AU 2007].
4. From the following table f (x) and hence f (6) using Newton’s divided difference formula.
x 1 2 7 8
[AU 2007].
y 1 5 5 4
x −4 −1 0 2 5
[AU 2004].
y 1245 33 5 9 1335
6. Using Newton’s divided difference formula, find f (x) from the following table and hence f (4).
x 0 1 2 5
[AU 2012].
f(x) 2 3 12 147
θ◦ 31 32 33 34
[AU 2004, 2011].
tanθ 0.6008 0.6249 0.6494 0.6745
2. Find the value of cos 1.74. Using the values given in the table below:
x 0 1 2 3 4
[AU 2012].
y 1 2.718 7.381 20.086 54.598
4. Find the first and second derivatives of y w.r.to x at x=10 from the data given below.
x 3 5 7 9 11
[AU 2008].
y 31 43 57 41 27
5. A jet fighter’s position on an aircraft carrier’s runway was timed during landing.
dy
where y is the distance from the end of the carrier. Estimate velocity dydt and acceleration dt2
at (i) t=1.1, (ii) t=1.6 using numerical differentiation. [AU 2007, 2008, 2010].
KCG\S&H 35
M A6452 QUESTION BANK
x 3 5 11 27 34
[AU 2011].
f(x) −13 23 899 17315 35606
2. Using the following data, f ′ (1), f ′ (5) and f ′ (6) [AU 2007, 2009, 2010, 2012].
x 0 2 3 4 7 9
f(x) 4 26 58 112 466 922
Rπ
1. Dividing the range into 10 equal parts, find the approximate value of sinxdx
0
by (i) Trapezoidal rule, (ii) Simpson’s rule. [AU 2006, 2010, 2012]
R2 dx
2. Evaluate x2 +4
using trapezoidal rule taking h = 0.25 [AU 2007] .
0
R2 dx
3. Evaluate x2 +x+1
to three decimal places, dividing the range of integration into 8 equal parts
0
using Simpson’s rule [AU 2012] .
4. A rocket is launched from the ground. Its acceleration is registered during the first 80 seconds
and it is in the table below. Using trapezoidal rule and Simpson’s 31 rule, find the velocity of
the rocket at t = 80secs. [AU 2010] .
t secs 0 10 20 30 40 50 60 70 80
f: cm/sec 30 31.63 33.34 35.47 37.75 40.33 43.25 46.69 40.67
5. Evaluate the length of the curve 3y = x3 from (0, 0) to (0, 13 ) rule and using 8 sub intervals
[AU 2010] .
R6 1
6. Evaluate 1+x2
dx using Trapezoidal rule. Verify the answer with direct integration.[AU 2010] .
0
7. The velocity v of a particle at a distance s from a point on its path is given by the table
s ft: 0 10 20 30 40 50 60
v: ft/sec: 47 58 64 65 61 52 38
1
Estimate the time taken to travel 60ft by using Simpson’s 3 rule. Compare the result with
Simpson’s 38 rule.[AU 2003] .
t sec : 0 2 4 6 8 10 12
v : m/sec : 4 6 16 34 60 94 136
KCG\S&H 36
M A6452 QUESTION BANK
Find the distance moved by the particle in 12 sec and also the acceleration at t = 2 sec.[AU 2003, 2011] .
R2 R2
1. Evaluate f (x, y)dxdy by Trapezoidal rule for the following data:
0 0
❍❍ x
❍❍ 0 0.5 1 1.5 2
y ❍
❍
0 2 3 4 5 5 [AU 2005] .
1 3 4 6 9 11
2 4 6 8 11 14
R2 R2 dxdy
2. Use trapezoidal rule to evaluate x+y taking 4 subintervals [AU 2004, 2007] .
1 1
R1 R1 1
3. Evaluate 1+x+y dxdy by Trapezoidal rule [AU 2009, 2011] .
0 0
R 1.4
1.2
1
R
4. Use trapezoidal rule to evaluate x+y dxdy by taking h = k = 0.1 [AU 2010] .
1 1
R1 R2 2xy
5. Evaluate (1+x2 )(1+y 2 )
dxdy by Trapezoidal rule with h = k = 0.25 [AU 2004, 2009] .
0 1
R2 R2 dxdy
6. Evaluate x+y by using Trapezoidal rule with h=k=0.5 and h=k=0.25. Improve the
1 1
estimate by Romberg’s formula[AU 2007] .
1
R1 R1 1
7. Using Simpson’s 3 rule, evaluate 1+x+y dxdy by taking h = k = 0.5. [AU 2004, 2009, 2013] .
0 0
R 4.4
2.4 R
8. Evaluate xydxdy using Simpson rule, given h = k = 0.1 [AU 2008, 2010] .
2 4
KCG\S&H 37
M A6452 QUESTION BANK
UNIT V
NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
PART A
Solution:
Given y ′ = f (x, y) = −y
x0 = 0, y0 = 1
By Euler’s method,
y1 = y0 + hf (x0 , y0 )
= y0 + h(−y0 )
= 1 − h(1)
=1−h
2. Write down the Milne’s predictor-corrector formula for solving initial value problem in first
order differential equation. [AU 2010,2011]
Solution:
4h
′ ′
yn+1 = yn−3 + 3 2yn−2 − yn−1 + 2yn′
h
′
yn+1 = yn−1 + 3 yn−1 + 4yn′ + yn+1
′
3. Given y ′ = −y and y(0) = 1, determine the value of y(0.01) by Euler’s method.[AU 2011]
Solution:
Given y ′ = f (x, y) = −y
x0 = 0, y0 = 1 h = 0.01
By Euler’s method,
y1 = y0 + hf (x0 , y0 )
= 1 + 0.01(−1)
= 1 − 0.01
= 0.99
∴ y(0.01) = 0.99
KCG\S&H 38
M A6452 QUESTION BANK
Solution:
k1 = hf (xn , yn )
k1
k2 = hf (xn + h2 , yn + 2 )
k2
k3 = hf (xn + h2 , yn + 2 )
k4 = hf (xn + h, yn + k3 )
[k1 +2k2 +2k3 +k4 ]
∆y = 6
yn+1 = yn + ∆y n = 0, 1, 2, 3, . . .
dy
5. Using Runge-Kutta method of fourth order to find the value of y when x=1, given that =
dx
y−x
with y(0) = 1 [AU 2014]
y+x
Solution:
dy y−x
Given = f (x, y) = , y(0) = 1. Take h = 1
dx y+x
Fourth order Runge-Kutta formula are
k1 = hf (xn , yn )
k1
k2 = hf (xn + h2 , yn + 2 )
k2
k3 = hf (xn + h2 , yn + 2 )
k4 = hf (xn + h, yn + k3 )
[k1 +2k2 +2k3 +k4 ]
∆y = 6
yn+1 = yn + ∆y n = 0, 1, 2, 3, . . .
Put n = 0 to find y1= y(1)
1−0
k1 = hf (x0 , y0 ) = 1 =1
1+0
!
3 1
k1 2 − 2
k2 = hf (x0 + h2 , y0 + 2 ) = 1f (0 + 12 , 1 + 21 ) = f ( 12 , 32 ) = 3 1 = 12
2 + 2
!
5 1
k2 4 − 2
k3 = hf (x0 + h2 , y0 + 2 ) = 1f (0 + 12 , 1 + 41 ) = f ( 12 , 54 ) = 5 1 = 37
4 + 2
!
10 3
7 −1
k4 = hf (x0 + h, y0 + k3 ) = 1f (0 + 1, 1 + 73 ) = f (1, 10
7 )= 10 = 17 7
= 3
17
7 +1 7
KCG\S&H 39
M A6452 QUESTION BANK
y1 = y0 + ∆y = 1 + 0.43 = 1.4342
∴ y(1) = 1.4342
6. Given that y ′ = y + x2 , y(0) = 1 find y(0.02) using modified Euler’s method. [AU 2014]
Solution:
Given h = 0.02 x0 = 0 y0 = 1
We have, MODIFIED EULER’S METHOD:
yn+1 = yn + hf xn + h2 , yn + h2 f (xn , yn ) ,
n = 0, 1, 2, 3, . . .
= 1.0202
7. State the merits of RK-method over Taylor series method. [AU 2012]
Solution:
Runge-kutta method is a single step method. This method does not require prior calculation
of higher derivatives likes Taylor’s series method. To find the value of y at x = xr+1 we
need the value of y at xr only.
dy d2 y
8. Write the central difference approximations for , .[AU 2012]
dx dx2
Solution:
yi+1 −yi
• yi′ = h
yi+1 −2yi +yi−1
• yi′′ = h2
KCG\S&H 40
M A6452 QUESTION BANK
9. Bring out the merits and demerits of Taylor series method. [AU 2012]
Solution:
• Taylor series method is a powerful single step method if we are able to get the
successive derivatives easily. This method is useful to give some initial values for
powerful numerical methods like Runge-kutta method.
• The disadvantage is that it became tedious if the higher derivatives are complicated.
dy
10. Find y(0.1) by Euler’s method, if = x2 + y 2 , y(0) = 1. [AU 2012]
dx
Solution:
dy
11. Using Taylor’s series method, find y at x = 0.1 given = x2 −y, y(0) = 1. [AU 2000,2006,2013]
dx
Solution:
Given y ′ = x2 − y, y(0) = 1
i.e., x0 = 0 and y0 = 1
Taylor’s series is
h ′ h2 ′′ h3 ′′′
yr+1 = yr + 1! yr + 2! yr + 3! yr + ... r = 0, 1, 2, 3, . . . − − − − − − − (1)
At (x0 , y0 ) = (0, 1)
∵ y ′ = x2 − y y ′ = 02 − 1 = 0 − 1 = −1
∴ y ′′ = 2x − y ′ y ′′ = 2(0) − (−1) = 1
y ′′′ = 2 − y ′′ y ′′′ = 2 − 1 = 1
Putting r = 0 in (1), we get
h ′ 2 3
y1 = y0 + 1! y0 + h2! y0′′ + h3! y0′′′ + . . .
∴ y(0.1) = 0.905
KCG\S&H 41
M A6452 QUESTION BANK
12. Compute y at x = 0.25 by Modified Euler method given y ′ = 2xy, y(0) = 1. [AU 2013]
Solution:
Given h = 0.25 x0 = 0 y0 = 1
yn+1 = yn + hf xn + h2 , yn + h2 f (xn , yn ) ,
n = 0, 1, 2, 3, . . .
dy
13. Using Taylor’s method, find y at x=1.1 given = x3 + y, Y (1) = 1. [AU 2013]
dx
Solution:
KCG\S&H 42
M A6452 QUESTION BANK
d2 y
14. Obtain the finite difference scheme for differential equation 2 2 + y = 5[AU 2013]
dx
Solution:
1
Given 2y ′′ + y = 5 ⇒ y ′′ = (5 − y)
2
1
⇒ yi+1 −2y
h2
i +yi−1
= (5 − yi )
2
h2
⇒ yi+1 − 2yi + yi−1 = (5 − yi )
2 2
h2
h
⇒ yi+1 + − 2 yi + yi−1 = 5
2 2
15. Convert the differential equation y ′′ (x) + y ′ (x) + y = 0 into finite difference equivalent form.
Solution:
16. How many pairs of prior values are required to predict the next value in Milne’s method?
Solution:
Four pairs of prior values are required to apply Milne’s Predictor formula.
17. State the Taylor’s series formula to find y(x1 ) for solving y ′ = f (x, y), y(x0 ) = y0 [AU 2007]
Solution:
h ′ h2 ′′ h3 ′′′
y1 = y(x1 ) = y0 + 1! y0 + 2! y0 + 3! y0 + ...
18. Write down the modified Euler’s method for ordinary differential equations. [AU 2007]
Solution:
yn+1 = yn + hf xn + h2 , yn + h2 f (xn , yn ) ,
n = 0, 1, 2, 3, . . .
19. State the difference between single step and multi step methods in solving ODE numerically?
KCG\S&H 43
M A6452 QUESTION BANK
Solution:
Euler and RK methods are single step methods, because they use only information from
the last step computed. But multistep methods use more than one past values to compute
the next value of y i.e., two or more step values are required for computation of the next.
dy
20. Write down the Euler algorithm to solve differential equation = f (x, y)
dx
Solution:
KCG\S&H 44
M A6452 QUESTION BANK
PART B
dy
1. Given = 1 + y 2 , y(0) = 0 by Taylor series method. Find the values of y at x = 0.2, x = 0.4
dx
and x = 0.6 [AU 2010]
2. Solve y ′ = x + y, y(0) = 1 by Taylor series method. Find the values of y at x = 0.1 and
x = 0.2 [AU 2005, 2006, 2010]
dy
3. Solve dx = x + y, y(0) = 1 by Taylor series method. Find y(1.1) and y(1.2) by Taylor series
method [AU 2010]
dy
4. Using Taylor series method, find y at x = 0.1 given dx = x2 − y, y(0) = 1 correct to 4 decimal
places.[AU 2000, 2006]
dy
5. Using Taylor method, find y(0.2) and y(0.4) given dx = 1 − 2xy, and y(0) = 1 by taking
h = 0.2 [AU 2011]
dy
6. Using Taylor series method solve dx = x2 − y, y(0) = 1 at x = 0.1, 0.2, 0.3 Also compare the
values with exact solution. correct to 4 decimal places.[AU 2012]
dy
7. Using Taylor series method find y at x = 0, if dx = x2 y − 1, y(0) = 1. [AU 2006, 2008]
8. Using Taylor series method find y(1.1) and y(1.2) correct to four places given y ′ = xy 1/3
[AU 2006]
9. Given y ′ = xy + y 2 , y(0) = 1, use Taylor series method to get the values of y(0.1), y(0.2) and
y(0.3). [AU 2009]
dy
10. Use Taylor series method to find y(0.1), y(0.2), given that dx = 3ex + 2y, y(0) = 0 correct to
4 decimal accuracy. [AU 2008, 2010]
dy
1. Using Euler’s method, find y(0.2) and y(0.4) from dx = x + y, y(0) = 1 with h = 0.2
[AU 2001, 2010, 2012]
2. Apply modified Euler’s method to find y(0.2) and y(0.4) given y ′ = x2 + y 2 , y(0) = 1 by
taking h = 0.2 [AU 2007, 2010]
dy
3. Using modified Euler method find y(0.1) and y(0.2) given dx = x2 + y 2 , y(0) = 1 by taking
h = 0.2 [AU 2010, 2011]
dy
4. Solve dx = log10 (x + y), y(0) = 2 by Euler’s modified method and find the values of y(0.2),
y(0.4), and y(0.6) by taking h = 0.2 [AU 2007]
5. Using modified Euler’s method solve given that y ′ = y − x2 + 1, y(0) = 0.5 find y(0.2)
[AU 2003, 2006, 2012]
dy √
6. Given that dx = 2 + xy and y = 1 when x = 1. Find approximate value of y at x = 2 in
steps of 0.2, using Euler’s modified method. [AU 2004]
KCG\S&H 45
M A6452 QUESTION BANK
2x
7. Using Eulers modified method compute y(0.1) with h = 0.1 from y ′ = y − y , y(0) = 1.
[AU 2011]
8. Using modified Euler method solve, given that y ′ = 1 − y, y(0) = 0 find y(0.1), y(0.2) and
y(0.3). [AU 2005]
dy
9. Using Modified Euler’s method, find y(4.1) and y(4.2) if 5x dx +y 2 −2 = 0; y(4) = 1. [AU 2012]
R-K method
1. Given y ′′ + xy ′ + y = 0, y(0) = 1, y ′ (0) = 0 find the value of y(0.1) by R.K method of fourth
order.[AU 2007, 2013]
2. Consider the second order initial value problem y ′′ − 2y ′ + 2y = e2t sin t with y(0) = −0.4 and
y ′ (0) = −0.6. Using fourth order Runge-Kutta method, find y(0.2) [AU 2003]
3. Compute y(0.1) and y(0.2) by Runge-Kutta method of 4th order for the differential equation
dy 2
dx = xy + y , y(0) = 1 [AU 2006, 2011]
dy dz
4. Solve for y(0.1) and z(0.1) from the simultaneous differential equations dx = 2y+z; dx = y−3z;
y(0) = 0, z(0) = 0.5 using Runge-Kutta method of fourth order. [AU 2012]
5. Solve y ′′ = xy ′ − y given y(0) = 3, y ′ (0) = 0 to find the value of y(0.1), using RK-method of
order 4. [AU 2012]
dy y 2 −x2
6. Using Runge-Kutta method of fourth order solve dx = y 2 +x2
with y(0) = 1 at x = 0.2, 0.4
[AU 2004, 2005, 2006, 2007, 2010, 2011, 2012, 2013]
7. Apply Runge-Kutta method to find approximate value of y for x = 0.2 in steps of 0.1 if
dy 2
dx = x + y given that y = 1 when x = 0. [AU 2008, 2010]
8. Find y(0.8) given that y ′ = y − x2 , y(0.6) = 1.739 by using Runge-Kutta method of fourth
order. Take h = 0.1 [AU 2000, 2012]
dy
9. Given that dx = x3 + y, y(0) = 2. Compute y(0.2), y(0.4), y(0.6) by Runge-Kutta method of
fourth order. [AU 2004]
dy y 2 −2x
10. Using Runge-Kutta method of fourth order solve dx = y 2 +x
with y(0) = 1 find y at x = 0.2
taking h = 0.2 [AU 2009]
dy
11. Apply 4th order Runge-Kutta method to determine y(0.2) with h = 0.1 from dx = x2 + y 2 ,
y(0) = 1. [AU 2000]
1. Solve y ′ = 21 (1 + x2 )y 2 , y(0) = 1, y(0.1) = 1.06, y(0.2) = 1.12, y(0.3) = 1.21. Compute y(0.4),
using Milne’s Predictor Corrector formula. [AU 2006, 2010]
2. Given 5xy ′ + y 2 = 2, y(4) = 1, y(4.1) = 1.0049, y(4.2) = 1.0097, y(4.3) = 1.0143 Compute
y(4.4) using Miline’s method. [AU 2004, 2007, 2010, 2011]
KCG\S&H 46
M A6452 QUESTION BANK
4. Given y ′ = xy + y 2 , y(0) = 1, find y(0.1) by Taylor’s method, y(0.2) by Euler’s method, y(0.3)
by Runge-Kutta method and y(0.4) by Milne’s method. [AU 2003, 2006, 2008]
dy
5. Given dx = x3 + y, y(0) = 2. y(0.2) = 2.073, y(0.4) = 2.452, y(0.6) = 3.023. Compute y(0.8)
by Milne’s predictor-corrector method taking h = 0.2 [AU 2004, 2006]
dy
6. Given dx = x2 (1 + y) and y(1) = 1, y(1.1) = 1.233, y(1.2) = 1.548, y(1.3) = 1.979, evaluate
y(1.4) by Milne’s predictor corrector method. [AU 2009]
7. Using Milne’s predictor-corrector method, find y(0.2) given y ′ = xy(1 + y) given y(−0.2) =
1.0412, y(0) = 1.0108 and y(−0.1) = 1.0108 [AU 2009]
dy
8. Given the initial value problem dx = 21 (x + y), y(0) = 2, y(0.5) = 2.636, y(1) = 3.595,
y(1.5) = 4.968 find y(2) by Milne method. [AU 2000, 2009]
9. Given y ′ = x(x2 + y 2 )e−x , y(0) = 1, find y at x = 0.1, 0.2 and 0.3 by Taylors series method
and compute y(0.4) by Milne’s method. [AU 2007]
Finite differences
∗ ∗ ∗ ∗ ∗ ∗ ∗∗
KCG\S&H 47