mth603 Solved Subjective Question For Mid Term
Uploaded by
kamran473473mth603 Solved Subjective Question For Mid Term
Uploaded by
kamran473473lOMoARcPSD|42838227
MTH603 Solved Subjective Question for mid term
Numerical analysis (Virtual University of Pakistan)
Scan to open on Studocu
Studocu is not sponsored or endorsed by any college or university
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.1
File Version v1.5.23 Update: (Dated: 23-May-2011)
This version of file contains:
Content of the Course (Done)
FAQ updated version.(These must be read once because some
very basic definition and question are being answered)
Glossary updated version.(These must be read once because
some very basic terms are being explained which you even
might not found in the book) (Done)
Solved Past Assignment Selected for MID Term. (Done)
Solved Question From Mid-Term Papers (Done)
4 Current Mid-Papers (Spring-2011) Included
MCQs GIGA File (with references) (Done)
Note: Use Table Of Content to view the Topics,
In PDF(Portable Document Format) format ,
you can check Bookmarks menu.
Disclaimer: There might be some human errors, if you find please let me know at
pak.nchd@gmail.com , duplication of data may be possible but at least possible level.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.2
----:Table of Content:----
Table of Content
FILE VERSION UPDATE: (DATED: 23-MAY-2011) ...................................................................................................................... 1
TABLE OF CONTENT ..................................................................................................................................................................... 2
INTRODUCTION TO NUMERICAL ANALYSIS .................................................................................................................................. 7
COURSE CONTENT:................................................................................................................................................................ 7
EACH TOPIC ONE SOLVED QUESTION WITH STEPS:COMPILED FROM PAST ASSIGNMENT............................................................ 7
PRELIMINARY MATERIAL: ................................................................................................................................................................ 7
Binary Number System (Base 2)=0,1 ..................................................................................................................................... 7
Octal Number System (Base 8)=0,1,2,3,4,5,6,7...................................................................................................................... 7
ERROR ................................................................................................................................................................................. 7
Decimal Number System (Base 10)=0,1,2,3,4,5,6,7,8,9 .......................................................................................................... 7
Hexa Number System (Base 16)=0,1,2,3,4,5,6,78,9,A,B,C,D,E,F ............................................................................................. 7
NON-LINEAR EQUATIONS ............................................................................................................................................................... 7
Method of Iterations .......................................................................................................................................................... 14
Bisection Method ................................................................................................................................................................. 7
Regula Falsi Method .......................................................................................................................................................... 11
Newton-Raphson Method .................................................................................................................................................. 16
Secant Method................................................................................................................................................................... 17
Muller’s Method ................................................................................................................................................................ 20
Graeffe’s Root Square Method ........................................................................................................................................... 21
LINEAR EQUATIONS ..................................................................................................................................................................... 21
Gaussian Elimination Method............................................................................................................................................. 21
Guass-Jordan Elimination Method ...................................................................................................................................... 24
Jacobi’s Iterative Method ................................................................................................................................................... 25
Gauss-Seidel Iteration Method ........................................................................................................................................... 26
Relaxation Method............................................................................................................................................................. 29
Matrix Inversion ................................................................................................................................................................. 29
EIGEN VALUE PROBLEMS .............................................................................................................................................................. 29
Power Method ................................................................................................................................................................... 29
Jacobi’s Method ................................................................................................................................................................. 32
INTERPOLATION .......................................................................................................................................................................... 36
Forward Differences ........................................................................................................................................................... 36
Backward Differences......................................................................................................................................................... 39
Divided Differences ............................................................................................................................................................ 41
By Newton’s Forward Difference Formula ........................................................................................................................... 45
b) By Lagrange’s Formula .................................................................................................................................................. 46
Langrange’s Interpolation .................................................................................................................................................. 47
Differentiation Using Difference Operators ......................................................................................................................... 47
NUMERICAL INTEGRATION ............................................................................................................................................................ 48
Trapezoidal Rule ................................................................................................................................................................ 48
Simpson’s 1/3 and 3/8 rules ............................................................................................................................................... 48
DIFFERENTIAL EQUATIONS ............................................................................................................................................................ 48
Taylor Series Method ......................................................................................................................................................... 48
Euler Method ..................................................................................................................................................................... 48
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.3
Runge-Kutta Method.......................................................................................................................................................... 48
Milne’s Predictor Corrector Method.................................................................................................................................... 48
Adam Moultan’s Predictor Corrector Method ..................................................................................................................... 48
FAQ UPDATED VERSION. ................................................................................................................................................ 48
QUESTION: WHAT IS BRACKETING METHOD? ............................................................................................................ 48
QUESTION: WHAT IS AN OPEN METHOD? ................................................................................................................... 48
QUESTION: EXPLAIN MULLER ’S METHOD BRIEFLY. ................................................................................................... 48
QUESTION: EXPLAIN THE DIFFERENCE BETWEEN THE LINEAR AND NON-LINEAR EQUATIONS. ............................... 48
QUESTION: EXPLAIN WHICH VALUE IS TO BE CHOOSED AS X0 IN N-R METHOD . .................................................... 49
QUESTION: DEFINE ITERATIVE METHOD OF SOLVING LINEAR EQUATIONS WITH TWO EXAMPLES. ......................... 49
QUESTION: DEFINE PIVOTING...................................................................................................................................... 49
QUESTION: WRITE THE TWO STEPS OF SOLVING THE LINEAR EQUATIONS USING GAUSSIAN ELIMINATION
METHOD. 49
QUESTION: DESCRIBE GAUSS-JORDAN ELIMINATION METHOD BRIEFLY.................................................................. 49
QUESTION: DESCRIBE BRIEFLY CROUT’S REDUCTION METHOD . .............................................................................. 49
QUESTION: DESCRIBE BRIEFLY THE JACOBI’S METHOD OF SOLVING LINEAR EQUATIONS. .................................... 49
QUESTION: WHAT IS THE DIFFERENCE BETWEEN JACOBI’S METHOD AND GAUSS SEIDAL METHOD?................... 49
QUESTION: WHAT IS THE BASIC IDEA OF RELAXATION METHOD? ........................................................................... 49
QUESTION: HOW THE FAST CONVERGENCE IN THE RELAXATION METHOD IS ACHIEVED? ...................................... 50
QUESTION: WHICH MATRIX WILL HAVE AN INVERSE? ............................................................................................... 50
QUESTION: WHAT ARE THE POPULAR METHODS AVAILABLE FOR FINDING THE INVERSE OF A MATRIX? .............. 50
QUESTION: EXPLAIN GAUSSIAN ELIMINATION METHOD FOR FINDING THE INVERSE OF A MATRIX. ...................... 50
QUESTION: WHAT ARE THE STEPS FOR FINDING THE LARGEST EIGEN VALUE BY POWER METHOD. ...................... 50
QUESTION: WHAT IS THE METHOD FOR FINDING THE EIGEN VALUE OF THE LEAST MAGNITUDE OF THE MATRIX
[A]? 50
QUESTION: WHAT IS INTERPOLATION? ....................................................................................................................... 50
QUESTION: WHAT IS EXTRAPOLATION? ...................................................................................................................... 50
QUESTION: WHAT HAPPENS WHEN SHIFT OPERATOR E OPERATES ON THE FUNCTION. .......................................... 50
QUESTION: WHAT IS THE BASIC CONDITION FOR THE DATA TO APPLY NEWTON’S INTERPOLATION METHODS? .. 50
QUESTION: WHEN IS THE NEWTON’S FORWARD DIFFERENCE INTERPOLATION FORMULA USED? .......................... 51
QUESTION: WHEN IS THE NEWTON’S FORWARD DIFFERENCE INTERPOLATION FORMULA USED? .......................... 51
QUESTION: WHEN THE NEWTON’S BACKWARD DIFFERENCE INTERPOLATION FORMULA IS USED?........................ 51
QUESTION: IF THE VALUES OF THE INDEPENDENT VARIABLE ARE NOT EQUALLY SPACED THEN WHICH FORMULA
SHOULD BE USED FOR INTERPOLATION? ........................................................................................................................... 51
QUESTION: TO USE NEWTON’S DIVIDED DIFFERENCE INTERPOLATION FORMULA, WHAT SHOULD THE VALUES OF
INDEPENDENT VARIABLES BE? ........................................................................................................................................... 51
QUESTION: WHICH DIFFERENCE FORMULA IS SYMMETRIC FUNCTION OF ITS ARGUMENTS? .................................. 51
QUESTION: IS THE INTERPOLATING POLYNOMIAL FOUND BY LAGRANGE’S AND NEWTON’S DIVIDED DIFFERENCE
FORMULAE IS SAME? ........................................................................................................................................................... 51
QUESTION: WHICH FORMULA INVOLVES LESS NUMBER OF ARITHMETIC OPERATIONS? NEWTON OR
LAGRANGE’S? 51
QUESTION: WHEN DO WE NEED NUMERICAL METHODS FOR DIFFERENTIATION AND INTEGRATION? ................... 51
QUESTION: IF THE VALUE OF THE INDEPENDENT VARIABLE AT WHICH THE DERIVATIVE IS TO BE FOUND APPEARS
AT THE BEGINNING OF THE TABLE OF VALUES, THEN WHICH FORMULA SHOULD BE USED? ......................................... 52
QUESTION: WHY WE NEED TO USE RICHARDSON’S EXTRAPOLATION METHOD? ..................................... 52
QUESTION: TO APPLY SIMPSON’S 1/3 RULE, WHAT SHOULD THE NUMBER OF INTERVALS BE? ............................ 52
QUESTION: TO APPLY SIMPSON’S 3/8 RULE, WHAT SHOULD THE NUMBER OF INTERVALS BE? ............................ 52
QUESTION: WHAT IS THE ORDER OF GLOBAL ERROR IN SIMPSON’S 1/3 RULE? .................................................... 52
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.4
QUESTION: WHAT IS THE ORDER OF GLOBAL ERROR IN TRAPEZOIDAL RULE? ....................................................... 52
QUESTION: WHAT IS THE FORMULA FOR FINDING THE WIDTH OF THE INTERVAL? ................................................. 52
QUESTION: WHAT TYPE OF REGION DOES THE DOUBLE INTEGRATION GIVE?.......................................................... 52
QUESTION: COMPARE THE ACCURACY OF ROMBERG’S INTEGRATION METHOD TO TRAPEZOIDAL AND SIMPSON’S
RULE. 52
QUESTION: WHAT IS THE ORDER OF GLOBAL ERROR IN SIMPSON’S 3/8 RULE? .................................................... 52
QUESTION: WHICH EQUATION MODELS THE RATE OF CHANGE OF ANY QUANTITY WITH RESPECT TO ANOTHER? 52
QUESTION: BY EMPLOYING WHICH FORMULA, ADAM-MOULTON P-C METHOD IS DERIVED? ................................. 52
QUESTION: WHAT ARE THE COMMONLY USED NUMBER SYSTEMS IN COMPUTERS? ................................................ 53
QUESTION: IF A SYSTEM HAS THE BASE M, THEN HOW MANY DIFFERENT SYMBOLS ARE NEEDED TO REPRESENT
AN ARBITRARY NUMBER? ALSO NAME THOSE SYMBOLS. ................................................................................................. 53
QUESTION: WHAT IS INHERENT ERROR AND GIVE ITS EXAMPLE. ............................................................................. 53
QUESTION: WHAT IS LOCAL ROUND-OFF ERROR? ...................................................................................................... 53
QUESTION: WHAT IS MEANT BY LOCAL TRUNCATION ERROR? .................................................................................. 53
QUESTION: WHAT IS TRANSCENDENTAL EQUATION AND GIVE TWO EXAMPLES. ...................................................... 53
QUESTION: WHAT IS MEANT BY INTERMEDIATE VALUE PROPERTY? ........................................................................ 53
QUESTION: WHAT IS DIRECT METHODS OF SOLVING EQUATIONS? ........................................................................... 53
QUESTION: WHAT IS ITERATIVE METHOD OF SOLVING EQUATIONS ? ....................................................................... 53
QUESTION: IF AN EQUATION IS A TRANSCENDENTAL, THEN IN WHICH MODE THE CALCULATIONS SHOULD BE
DONE? 53
QUESTION: WHAT IS THE CONVERGENCE CRITERION IN METHOD OF ITERATION? ................................................. 54
QUESTION: WHEN WE STOP DOING ITERATIONS WHEN TOLL IS GIVEN? .................................................................. 54
QUESTION: WHAT IS A TRANCE DENTAL EQUATION? ........................................................ ERROR! BOOKMARK NOT DEFINED.
QUESTION: HOW THE VALUE OF H IS CALCULATED IN INTERPOLATION? ................................................................. 54
QUESTION: WHAT IS AN ALGEBRAIC EQUATION? ....................................................................................................... 54
QUESTION: WHAT IS DESCARTES RULE OF SIGNS?.................................................................................................... 54
QUESTION: WHAT ARE DIRECT METHODS? ................................................................................................................ 54
QUESTION: WHAT IS MEANT BY ITERATIVE METHODS? ............................................................................................. 54
QUESTION: WHAT IS GRAPHICALLY MEANT BY THE ROOT OF THE EQUATION? ....................................................... 54
QUESTION: Q. WHAT IS THE DIFFERENCE BETWEEN OPEN AND BRACKETING METHOD? ....................................... 54
QUESTION: CONDITION FOR THE EXISTENCE OF SOLUTION OF THE SYSTEM OF EQUATIONS. ................................ 54
QUESTION: SHOULD THE SYSTEM BE DIAGONALLY DOMINANT FOR GAUSS ELIMINATION METHOD? ..................... 54
QUESTION: WHAT IS MEANT BY DIAGONALLY DOMINANT SYSTEM? .......................................................................... 55
QUESTION: STATE THE SUFFICIENT CONDITION FOR THE CONVERGENCE OF THE SYSTEM OF EQUATION BY
ITERATIVE METHODS........................................................................................................................................................... 55
QUESTION: THE CALCULATION FOR NUMERICAL ANALYSIS SHOULD BE DONE IN DEGREE OR RADIANS. ............... 55
QUESTION: HOW WE CAN IDENTIFY THAT NEWTON FORWARD OR BACKWARDS INTERPOLATION FORMULA IS TO
BE USED. 55
QUESTION: WHAT IS MEANT PRECISION AND ACCURACY? ......................................................................................... 55
QUESTION: WHAT IS THE CONDITION THAT A ROOT WILL LIE IN AN INTERVAL........................................................ 55
QUESTION: HOW THE DIVIDED DIFFERENCE TABLE IS CONSTRUCTED? ................................................................... 55
QUESTION: WHAT IS GAUSS-SEIDEL METHOD ............................................................................................................ 55
QUESTION: WHAT IS PARTIAL AND FULL PIVOTING? .................................................................................................. 56
QUESTION: HOW THE INITIAL VECTOR IS CHOOSE IN POWER METHOD? .................................................................. 56
QUESTION: WHAT IS THE RELATION SHIP BETWEEN P=0 AND NON ZERO P IN INTERPOLATION. ............................. 56
QUESTION: WHAT IS CHOPPING AND ROUNDING OFF? .............................................................................................. 56
QUESTION: WHEN THE FORWARD AND BACKWARD INTERPOLATION FORMULAE ARE USED? ................................. 56
QUESTION: WHAT IS FORWARD AND BACKWARD DIFFERENCE OPERATOR AND THE CONSTRUCTION OF THEIR
TABLE. 57
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.5
QUESTION: WHAT IS JACOBI’S METHOD?................................................................................................................... 57
QUESTION: WHAT IS SIMPSON'S 3/8TH RULE. ........................................................................................................... 57
QUESTION: WHAT IS CLASSIC RUNGE-KUTTA METHOD ............................................................................................... 57
QUESTION: WHAT IS MEANT BY TOL? ........................................................................................................................ 57
QUESTION: WHAT IS MEANT BY UNIQUENESS OF LU METHOD. .................................................................................. 58
QUESTION: HOW THE VALU OF H IS CALCULATED FROM EQUALLY SPACED DATA. ................................................... 58
GLOSSARY (UPDATED VERSION) .................................................................................................................................. 58
Absolute Error : .................................................................................................................................................................. 58
Accuracy : .......................................................................................................................................................................... 58
Algebraic equation : ........................................................................................................................................................... 58
Bisection method : ............................................................................................................................................................. 58
Bracketing Method : .......................................................................................................................................................... 58
Crout’s Method : ................................................................................................................................................................ 58
Direct methods :................................................................................................................................................................. 59
Gauss Elimination Method : ............................................................................................................................................... 59
Gauss Seidel iterative method : .......................................................................................................................................... 59
Graeffee’s root squaring method : ...................................................................................................................................... 59
Intermediate value property : ............................................................................................................................................. 59
Inverse of a matrix : ........................................................................................................................................................... 59
Iterative methods : ............................................................................................................................................................. 59
Jacobie’s iterative method: : ............................................................................................................................................... 59
Muller’s Method : .............................................................................................................................................................. 59
Newton Raphson Method. : ................................................................................................................................................ 59
Non singular matrix :.......................................................................................................................................................... 59
Open methods :.................................................................................................................................................................. 60
Pivoting : ........................................................................................................................................................................... 60
Precision : .......................................................................................................................................................................... 60
Regula –Falsi method : ....................................................................................................................................................... 60
Relaxation method : ........................................................................................................................................................... 60
Secant Method :................................................................................................................................................................. 60
Significant digits :............................................................................................................................................................... 60
Singular matrix :................................................................................................................................................................. 60
Transcendental equation : .................................................................................................................................................. 61
Truncation Error : ............................................................................................................................................................... 61
IMPORTANT FORMULA FOR MTH603 ........................................................................................................................................ 61
Bisection Method ............................................................................................................................................................... 61
Muller Method ................................................................................................................................................................... 61
Regula Falsi Method (Method of False position) ................................................................................................................. 61
Newton Rophson method ................................................................................................................................................... 61
Secant Method................................................................................................................................................................... 61
Newton’s Formula .............................................................................................................................................................. 62
Graffee root squaring method ............................................................................................................................................ 62
SPRING 2011 LATEST PAPERS (CURRENT)................................................................................................................................... 62
PAPERS NUMBER 01 ................................................................................................................................................................... 62
Question: What are the steps for finding the largest eigen value by power method? ..................................................... 62
Jacobi method only ist iteration calculate kerni thi 2 marks................................................................................................. 62
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.6
Jacobi method 2 iteration calculation 5 marks .................................................................................................................... 63
Newton Backward difference method 5 marks.................................................................................................................... 63
Backward difference table construct 2 marks ..................................................................................................................... 63
Backward difference table construct 3 marks ..................................................................................................................... 63
PAPERS NUMBER 02 ................................................................................................................................................................... 63
1 question was to find the value of p from the given table using Interpolation Newton's forward Difference Formula. ......... 63
2) we have to calculate the values of x, y and z for the given system of equations using Gauss Siede iterative method taking
initial solution vector as (0, 0, 0)t? ...................................................................................................................................... 63
Both of the questions were to find the value of theta for the given matrices using Jacobi's method ? 2x3 marks ................ 63
1) Construct backward difference table for the given values ….5 marks ............................................................................... 63
2) Find the solution of following system of Equations using Jacobi's iterative method upto three decimal places ? 5 marks . 63
PAPERS NUMBER 03 ................................................................................................................................................................... 64
Regula Falsi method ka aik 5 marks ka sewal...................................................................................................................... 64
Orthogonal metrix ki defination 3 mark ki .......................................................................................................................... 64
Crute's method ko briefly define kerna tha 3 marks ............................................................................................................ 64
SOLVED_MID-TERM PAST PAPERS ............................................................................................................................................. 65
SHORT QUESTION (SET-1) ............................................................................................................................................................ 65
SHORT QUESTION (SET-2) ............................................................................................................................................................ 65
SHORT QUESTION (SET-4) ............................................................................................................................................................ 66
SHORT QUESTION (SET-5) ............................................................................................................................................................ 66
SOLVED_MCQZ .......................................................................................................................................................................... 66
MCQZ (SET-1) .......................................................................................................................................................................... 66
MCQZ (SET-2) .......................................................................................................................................................................... 67
MCQZ (SET-3) ................................................................................................................................. ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-4) ................................................................................................................................. ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-5) ................................................................................................................................. ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-6) .......................................................................................................................................................................... 69
MCQZ (SET-7) ................................................................................................................................. ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-8) .......................................................................................................................................................................... 71
MCQZ (SET-9) ................................................................................................................................. ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-10) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-11) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-12) ........................................................................................................................................................................ 73
MCQZ (SET-13) ........................................................................................................................................................................ 73
MCQZ (SET-14) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-15) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-16) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
MCQZ (SET-17) ............................................................................................................................... ERROR! BOOKMARK NOT DEFINED.
=========================================================>
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.7
Introduction To Numerical Analysis
Course Content:
Number systems, Errors in computation, Methods of solving non-linear equations,
Solution of linear system of equations and matrix inversion, Eigen value problems, power
method, Jacobi’s method, Different techniques of interpolation, Numerical differentiation
and integration, Numerical integration formulas, different methods of solving ordinary
differential equations.
Each Topic One Solved Question with Steps:Compiled from Past assignment
Preliminary material:
Representation of numbers
Binary Number System (Base 2)=0,1
Used in Computer Internal Operations and Calculation
Octal Number System (Base 8)=0,1,2,3,4,5,6,7
ERROR
Types of error are,
Inherent errors,
Local round-off errors
Local truncation errors
Decimal Number System (Base 10)=0,1,2,3,4,5,6,7,8,9
World Wide used and the most common system
Hexa Number System (Base 16)=0,1,2,3,4,5,6,78,9,A,B,C,D,E,F
Used in Processor Register Addressing
Errors in computations
The initial approximation
It may be found by two methods either by graphical method or
Analytical method
Graphical method
Non-Linear Equations
Bisection Method
Procedure in Detail:
Step,Take two Initial Approximation such that f(x1).F(x2)<0.Means both must have
opposite signs.Take their mean by x3=(x1+x2)/2
Next Take two element from which 1 will be x3 and another from x1 or x2 such that both
x3 and the other element should have opposite sign.
Repeat the above process to the required numbers of iterations.
Question 1
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.8
Find the root of the equation given below by bisection method.
x3 x 2 x 7 0
(Note: accuracy up to three decimal places is required.) Marks: 10
SOLUTION
Let
f ( x) x3 x 2 x 7
Now f (1) (1)3 (1) 2 1 7 6 0
f (2) (2)3 (2) 2 2 7 1 0
and f (3) (3)3 (3) 2 3 7 14 0
As f (2) f (3) 14 0
Therefore a real root lies between 2 and 3
Iteration 1
Let
x0 2 and x1 3
then
x0 x1
x2 ( )
2
23
( ) 2.5
2
Now f (2.5) (2.5)3 (2.5) 2 2.5 7 4.875 0
As f (2) f (2.5) 1(4.875) 4.875 0
Therefore a real root lies between 2 and 2.5
Iteration 2
x0 2 and x2 2.5
then
x0 x2
x3 ( )
2
2 2.5
( ) 2.25
2
Now f (2.25) (2.25)3 (2.25) 2 2.25 7 1.578125 0
As f (2) f (2.25) 1(1.578125) 1.578125 0
Therefore a real root lies between 2 and 2.25
Iteration 3
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.9
x0 2 and x3 2.25
then
x0 x3
x4 ( )
2
2 2.25
( ) 2.125
2
Now f (2.125) (2.125)3 (2.125) 2 2.125 7 0.205078 0
As f (2) f (2.125) 1(0.205078) 0.205078 0
Therefore a real root lies between 2 and 2.125
Iteration 4
x0 2 and x4 2.125
then
x0 x4
x5 ( )
2
2 2.125
( ) 2.0625
2
Now f (2.0625) (2.0625)3 (2.0625) 2 2.0625 7 0.417725 0
As f (2.0625) f (2.125) 0.417725(0.205078) 0.085666 0
Therefore a real root lies between 2.0625 and 2.125
Iteration 5
x5 2.0625 and x4 2.125
then
x5 x4
x6 ( )
2
2.0625 2.125
( ) 2.09375
2
Now f (2.09375) (2.09375)3 (2.09375)2 2.09375 7 0.111481 0
As f (2.09375) f (2.125) 0.111481(0.205078) 0.022862 0
Therefore a real root lies between 2.09375 and 2.125
Iteration 6
x6 2.09375 and x4 2.125
then
x6 x4
x7 ( )
2
2.09375 2.125
( ) 2.109375
2
Now f (2.109375) (2.109375)3 (2.109375) 2 2.109375 7 0.0455 0
As f (2.09375) f (2.109375) 0.111481(0.0455) 0.00507 0
Therefore a real root lies between 2.09375 and 2.109375
Iteration 7
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.10
x6 2.09375 and x7 2.109375
then
x6 x7
x8 ( )
2
2.09375 2.109375
( ) 2.1015625
2
Now f (2.1015625) (2.1015625)3 (2.1015625) 2 2.1015625 7 0.0333 0
As f (2.1015625) f (2.109375) 0.0333(0.0455) 0.01515 0
Therefore a real root lies between 2.1015625 and 2.109375
Iteration 8
x8 2.1015625 and x7 2.109375
then
x8 x7
x9 ( )
2
2.1015625 2.109375
( ) 2.10546875
2
Now f (2.10546875) (2.10546875)3 (2.105468752 2.10546875 7 0.006 0
As f (2.1015625) f (2.10546875) 0.0333(0.006) 0.0002 0
Therefore a real root lies between 2.1015625 and 2.10546875
Iteration 9
x8 2.1015625 and x9 2.10546875
then
x8 x9
x10 ( )
2
2.1015625 2.10546875
( ) 2.103515625
2
Now f (2.103515625) (2.103515625)3 (2.103515625)2 2.103515625 7 0.01367 0
As f (2.103515625) f (2.10546875) 0.01367(0.006) 0.00008 0
Therefore a real root lies between 2.103515625 and 2.10546875
Iteration 10
x10 2.103515625 and x9 2.10546875
then
x10 x9
x11 ( )
2
2.103515625 2.10546875
( ) 2.104492188
2
Now f (2.104492188) (2.104492188) 3 (2.104492188) 2 2.104492188 7 0.0038 0
As f (2.104492188) f (2.10546875) 0.0038(0.006) 0.0000228 0
Therefore a real root lies between 2.104492188 and 2.10546875
Iteration 11
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.11
x11 2.104492188 and x9 2.10546875
then
x11 x9
x12 ( )
2
2.104492188 2.10546875
( ) 2.104980469
2
Now f (2.104980469) (2.104980469) 3 (2.104980469) 2 2.104980469 7 0.0011 0
As f (2.104492188) f (2.104980469) 0.0038(0.0011) 0.00000418 0
Therefore a real root lies between 2.104492188 and 2.104980469
As the next root lies between 2.104492188 and 2.104980469 and these roots are equal up to
three decimal places.
So, the required root up to three decimal places is 2.104
Example From Handout at page # 4
Solve x3 −9x+1=0 for the root between x=2 and x=4 by bisection method
Solution:
Here we are given the interval (2,4) so we need not to carry out intermediate value
property to locate initial approximation.
Here f(x) =x 3 9 x 1 0
Now f(2)=-9 and f(4)=29.
Here f(2).f(4)<0 So the root lies between 2 and 4.
So, xo 2, x1 4
xo x 1
and x 2 =3 (This formula predicts next iteration)
2
Now f(3)=1,Here f(2).f(3)<0 So the root lies between 2 and 3.
Repeat above process to required number of iteration.
Step,Take two Initial Approximation such that f(x1).F(x2)<0.Means both must have
opposite signs.Take their mean by x3=(x1+x2)/2
Next Take two element from which 1 will be x3 and another from x1 or x2 such that both
x3 and the other element should have opposite sign.
Repeat the above process to the required numbers of iterations.
Regula Falsi Method
x 2 x1
Formula for the Regula Falsi Method is x3 x 2 f ( x 2) or
f ( x 2) f ( x1)
xn xn 1
x( n1) xn f xn
f xn f xn1
Steps: If the interval is given , check whether the root lies in between the interval or not
by the principle that states if both number have opposite sign , then the root lies in
between interval.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.12
Find the next approximation with the help of the formula.
Question#2 Marks 10
Use Regula-Fasli method to compute the root of the equation f ( x) cos x xe In the
x
interval [0, 1] after third iteration.
Solution:
As f ( x) cos x xe x
f (0) cos 0 0e0 1 0
f (1) cos1 1e1 0.5403 2.718 2.1779 0
So root of the eq. lie between 0 and 1.Let x0 0 and x1 1
Thus f ( x0 ) 1 and f ( x1 ) 2.1779
The formula for finding the root of the function f(x) by Regula-Fasli is given by.
xn xn 1
x( n1) xn f xn
f xn f xn1
For n=1 we have,
x1 xo 1 0 2.17797
x2 x1 f x1 1 -2.17797 1 1 0.68533 0.31467
f x1 f xo -2.17797 1 3.17797
f x2 cos(0.31467) 0.31467 e0.31467 0.950899 0.31467 1.369806 =0.519863
Since f (x1) and f (x2) are of opposite signs, the root lies between x1 and x2.
Therefore, for n=2 we have,
x2 x1 0.31467 1
x3 x2 f x2 0.31467 0.519863
f x2 f x1 0.519863 -2.17797
0.68533
0.31467 0.519863 0.44673
2.697833
f x3 cos(0.44673) 0.44673 e0.44673 0.901864 0.44673 1.563191 =0.20353
Since f (x1) and f (x3) are of opposite signs, the root lies between x1 and x3.
Therefore, for n=3 we have,
x3 x1 0.44673 1
x4 x3 f x2 0.44673 0.20353
f x3 f x1 0.20353 -2.17797
0.54327
0.44673 0.20353 0.493159
2.3815
The required root after 3rd iteration using Regula-Falsi method is 0.493159
Question 2
Use the Regula Falsi (method of false position) to solve the equation x3 4 x 9 0
(Note: accuracy up to four decimal places is required) Marks: 10
SOLUTION
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.13
Let
f ( x) x 3 4 x 9
Now
f (0) (0)3 4(0) 9 0 0 9 9 0
f (1) (1)3 4(1) 9 1 4 9 12 0
f (2) (2)3 4(2) 9 8 8 9 9 0
f (3) (3)3 4(3) 9 27 12 9 6 0
Now, Since f (2) and f (3) are of opposite signs,
therefore, the real root lies between 2 and 3
Now,
let x1 2 and x2 3
First iteration
x2 x1
x3 x2 f ( x2 )
f ( x2 ) f ( x1 )
3 2
3 (6)
6 (9)
1
3 (6)
15
2
3
5
2.6
f ( x3 ) f (2.6) (2.6)3 4(2.6) 9
1.824
Second iteration
x3 x2
x4 x3 f ( x3 )
f ( x3 ) f ( x2 )
2.6 3
2.6 (1.824)
1.824 6
2.69325
f ( x4 ) f (2.69325) (2.69325)3 4(2.69325) 9
0.23725
Third iteration
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.14
x4 x3
x5 x4 f ( x4 )
f ( x4 ) f ( x3 )
2.69325 2.6
2.69325 0.2061
0.23725 (1.824)
2.70536
f ( x5 ) f (2.70536) (2.70536)3 4(2.70536) 9
0.02098
Fourth iteration
x5 x4
x6 x5 f ( x5 )
f ( x5 ) f ( x4 )
2.70536 2.69325
2.70536 (0.02098)
0.02098 (0.23725)
2.70653
f ( x6 ) f (2.70653) (2.70653)3 4(2.70653) 9
0.0000368
Fifth iteration
x6 x5
x7 x6 f ( x6 )
f ( x6 ) f ( x5 )
2.70653 2.70536
2.70653 (0.0000368)
0.0000368 (0.02098)
2.70653
Hence,
The required root of the given equation correct to 4 decimal places is 2.7065
Method of Iterations
Question#1 Marks 10
x
Use the method of iterations to determine the real root of the equation e 10 x in the
interval [0, 1], correct to four decimal places after four Iterations.
Solution:
Since e x 10 x.Let f ( x ) e x 10 x.we see that f (0) 1 and f (1) 9.6321
As e x 10 x ,we can easily find the value of x, thus
e x e x
x , Let x
10 10
By taking its derivative we get
e x
x , we see that x 1 for all values in [0,1].
10
Therefore we can apply method of iterations to the given function.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.15
Take any value within [0,1].Let
x0 0.5
e 0.5
x1 0.5 0.0607 , f ( x1 ) 0.09516
10
e 0.0607
x2 0.0607 0.0914 , f ( x2 ) 0.00869
10
e 0.0914
x3 0.0914 0.0941 , f ( x3 ) 0.000790477 7.90477 x10 4
10
e 0.0941
x4 0.0941 0.0913 , f ( x4 ) 0.0000275784 2.75784 x10 5
10
x4 0.0913
Example From the Handout at page #
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.16
Newton-Raphson Method
Procedural Detail:
Find the limit if not provided by starting from x=0 to ,…
Find two consecutive numbers for f(x) should have opposite sign.
Use the Newton Raphson Formula to find next approximation .
Question#3 Marks 10
Find the real root of the equation x3 3x 5 0 using Newton-Raphson method in the
interval [2,3] after third iteration.
Solution:
As f ( x) x 3 3x 5
f (2) 23 3.2 5 3 0
f (3) 3 3.3 5 13
3
0
So root of the eq. will lie in 2,3
now f / ( x) 3 x 2 3 and f / / ( x) 6 x
f / (2) 3.22 3 9 and f / (3) 3.32 3 24
f / / (2) 6.2 12 and f / / (3) 6.3 18
Since f (3) and f / / (3)are of same sign.So we choose x0 3
so by Newton ' s method we have
f (x ) 13
x1 x0 / 0 3 2.4583
f ( x0 ) 24
f ( x1 ) (2.4583)3 3(2.4583) 5 14.8561 7.3749 5 2.4812
f / ( x1 ) f / (2.4583) 3(2.4583) 2 3 18.1297 3 15.1297
f ( x1 ) 2.4812
x2 x1 /
2.4583 2.2943
f ( x1 ) 15.1297
f ( x2 ) (2.2943)3 3(2.2943) 5 12.0767 6.8829 5 0.1938
f / ( x2 ) f / (2.2943) 3(2.2943) 2 3 15.7914 3 12.7914
f ( x2 ) 0.1938
x3 x2 /
2.2943 2.2791
f ( x2 ) 12.7914
Example From the Handout at page #
Value of e=2.7182
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.17
Secant Method
Question 1
Do Four iterations of Secant method, with an accuracy of 3 decimal places to find the
.
root of
f ( x) x 3 3x 1 0, x0 1, x1 0.5 Marks: 10
Solution:
FORMULA OF SECANT METHOD
xn1 xn1 ( fxn ) xn f ( xn 1)
( fxn ) f ( xn 1)
ITERATION 1
n 1
f ( x0 ) f 1 1
f ( x1 ) f 0.5 0.375
x2 xo ( fx1 ) x1 f ( xo )
( fx1 ) f ( xo )
x2 (1)( 0.375) (0.5)( 1)
( 0.375) ( 1)
x2 (1)( 0.375) (0.5)(1)
( 0.375) (1)
x2 ( 0.375) (0.5)
0.625
x2 0.125
0.625
x2 0.2
ITERATION 2
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.18
n 2
f ( x2 ) f 0.2 0.408
x3 x1 ( fx2 ) x2 f ( x1 )
( fx2 ) f ( x1 )
x3 (0.5)(0.408) (0.2)( 0.375)
(0.408) ( 0.375)
x3 (0.204) (0.075)
(0.408) (0.375)
x3 0.279
0.783
x3 0.3563
ITERATION 3
n 3
f ( x3 ) f 0.3563 0.02367
x4 x2 f ( x3 ) x3 f ( x2 )
f ( x3 ) f ( x2 )
x4 (0.2) f (0.3563) (0.3563) f (0.408)
f (0.3563) f (0.408)
x4 (0.2)( 0.02367) (0.3563)(0.408)
( 0.02367) (0.408)
x4 0.004737 0.1453
0.43167
x4 0.150034
0.43167
x4 0.3477
ITERATION 4
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.19
n 4
f ( x3 ) f 0.3477 0.02367
x4 0.3477 f ( x4 ) 0.00107
x5 x3 f ( x4 ) x4 f ( x3 )
f ( x4 ) f ( x3 )
x5 (0.3563)(0.00107) (0.3563)(0.02367)
( 0.00107) ( 0.02367)
x5 0.000381 0.0823
0.00107 0.02367
x5 0.7849
0.0226
x5 0.3473
x5 0.3473 f ( x5 ) 0.0000096
Hence, the root after four iterations is 0.347
================================================================
Question 2
Use the secant method to solve the equation e x 3x 2 for 0 x 1 . (Perform only 3 iterations.)
Solution
f ( x) e x 3 x 2 , x0 0 , x1 1
both aretheinitialapproximations
f (0) e0 3(0) 1
f (1) e1 3(1) 0.281
nowwecalculatethesec ond approximation
x f ( x1 ) x1 f ( x0 ) (0)(0.281) (1)1
x2 0 0.7806
f ( x0 ) f ( x1 ) 0.281 1
f (0.7806) e0.7806 3(0.7806) 2 2.1827 1.82800 0.3546
now x1 1 and x2 0.786
f (1) 0.281 f (0.786) 0.3546
x f ( x2 ) x2 f ( x1 ) 1(0.3546) (0.786)( 0.281)
x3 1 0.9052
f ( x2 ) f ( x1 ) 0.3546 0.281
now
x3 0.9052f(0.9052)
x2 0.786 f (0.786) 0.3546
x2 f ( x3 ) x3 f ( x2 ) 0.7806(0.0074) (0.9052)(0.3546)
x4 0.9076
f ( x3 ) f ( x2 ) 0.0074 0.3546
sorootis0.9076
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.20
Muller’s Method
Question 1
Solve the equation x 7 x 14 x 6 by using Muller’s method only perform three
3 2
Iterations.( x0 0.5, x1 1, x2 0 )
Solution
1st iteration
x0 0.5 , x1 1 and x2 0
f ( x0 ) f 0 f (0.5) (0.5)3 7(0.5) 2 14(0.5) 6 0.625
f ( x1 ) f1 f (1) 13 7(1) 2 14(1) 6 2
f ( x2 ) f 2 f (0) (0)3 7(0) 2 14(0) 6 6
c f 0 0.625
h1 x1 x0 1 0.5 0.5
h2 x0 x2 0.5 0 0.5
h2 f1 (h1 h2 ) f 0 h1 f 2
a
h1h2 (h1 h2 )
(0.5)(2) (0.5 0.5)(0.625) (0.5)(6) 1 (1)( 0.625) (0.5)(6)
a = 5.5
(0.5)(0.5)(0.5 0.5) (0.5)(0.5)(1)
f1 f 0 ah12
b
h1
2 (0.625) (5.5)(0.5) 2 2 0.625) (5.5)(0.25) 2
b 8
0.5 0.5
2c
x x0
b b 2 4ac
2(0.625)
x 0.5
8 (8) (4)(5.5)(0.625)
2
x 1.8721
2nd iteration
x0 0.5 , x1 1.8721 and x2 1
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.21
f ( x0 ) f 0 f (0.5) (0.5)3 7(0.5) 2 14(0.5) 6 0.625
f ( x1 ) f1 f (1) 1.87213 7(1.8721) 2 14(1.8721) 6 2.228
f ( x2 ) f 2 f (0) (0)3 7(0) 2 14(0) 6 6
c f 0 0.625
h1 x1 x0 1.8721 0.5 1.3721
h2 x0 x2 0.5 1 0.5
h2 f1 (h1 h2 ) f 0 h1 f 2
a
h1h2 (h1 h2 )
(0.5)(2.228) (1.3721 0.5)(0.625) (1.3721)( 6)
a 14.70
(1.3721)(0.5)(1.3721 0.50.5)
f1 f 0 ah12
b
h1
2.228 (0.625) (14.70)(1.3721) 2
b 17.72
0.5
2c
x x0
b b 2 4ac
1.25
x 0.5 0.29
12.62 (17.72)2 (4)(14.71)(0.625)
x 1.8721
For 3rd iteration
x0 0.29 , x1 1.8721 , x2 1
For third iteration we will proceed in the same manner.
Graeffe’s Root Square Method
Linear Equations
Gaussian Elimination Method
Question 2
Using Gaussian Elimination Method, solve the following system of equations
x1 x2 2 x3 9
x1 3x2 2 x3 13
3x1 x2 3x3 14
Solution:
Marks: 10
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.22
The Augmented Matrix of the given system of equations is
1 1 2 9
Ab 1 3 2 13
3 1 3 14
1 1 2 9
R 0 2 0 4 by R2 R1 and R3 3R1
0 2 3 13
1 1 2 9
R 0 1 0 2 by
1
R2
2
0 2 3 13
1 1 2 9
R 0 1 0 2 by R 3 (2) R2
0 0 3 9
1 1 2 9
R 1
0 1 0 2 by R3
3
0 0 1 3
Which shows that from the third and second rows
Z=3, y = 2
And from the first row
X+y+2z=9
Using the values of y and z, we get x = 1
Hence the solution of the given system is
X = 1, y = 2, z = 3
Question 1
Using Gaussian Elimination Method, solve the following system of equations
2x y 2z 2
x 10 y 3z 5
x yz 3
Marks: 10
Solution:
The Augmented Matrix of the given system of equations is
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.23
2 1 2 2
Ab 1 10 3 5
1 1 1 3
1 10 3 5
R 2 1 2 2 by R12
1 1 1 3
1 10 3 5
R 0 21 8 8 by R2 2 R1 , R3 R1
0 11 2 2
1 10 3 5
8 8 1
R
0 1 by R2
21 21 21
0 11 2 2
1 10 3 5
0 8 8
R
1 by R3 11R2
21 21
46 46
0 0
21 21
Using Gaussian Elimination method, by backward substitution, we get as follows
From the third row, we get
46 46
z
21 21
z 1
From the second row, we get
8 8
y z
21 21
8 8
y z
21 21
Putting the value of z , we get
8 8
y (1)
21 21
8 8
21 21
0
And finally from the first row, we get
x 10 y 3z 5
Putting the values of y and z, we get
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.24
x 10(0) 3( 1) 5
x35
x2
So, the solution is
x 2, y 0, z 1
Guass-Jordan Elimination Method
Question#1 Marks 10
By using Gauss –Jordan elimination method, solve the following system of equations,
x yz 7
3x 3 y 4 z 24
2 x y 3z 16
Solution:
The given system in matrix form is
1 1 1 x 7
3 3 4 y 24
2 1 3 z 16
A X B
1 1 1 : 7
A | B 3 3 4 : 24
2 1 3 : 16
1 1 1 : 7
2 1 3 : 16 R23
3 3 4 : 24
1 1 1 : 7
0 1 1 : 2 R2 R2 2 R1 , R3 R3 3R1
0 0 1 : 3
1 0 2 : 9
0 1 1 : 2 R1 R1 R2
0 0 1 : 3
1 0 0 : 3
0 1 0 : 1 R1 R1 2 R3 , R2 R2 R3
0 0 1 : 3
1 0 0 : 3
0 1 0 : 1 R2
0 0 1 : 3
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.25
So,
x 3, y 1, z 3
Jacobi’s Iterative Method
Solve the following system of equations
20 x y 2 z 17
3x 20 y z 18
2 x 3 y 20 z 25
By Jacobi’s iterative method taking the initial starting of solution vector as (0, 0, 0)T
and perform the first three iterations.
Solution:
20 x y 2 z 17
3x 20 y z 18
2 x 3 y 20 z 25
17 y 2 z
x
20
18 3 x z
y
20
25 2 x 3 y
z
20
startingwith(0,0,0)
Iteration#01
17 0 2 0 17
x 0.85
20 20
18 3 0 0 18
y 0.9
20 20
25 2 0 3 0 25
z 1.25
20 20
Iteration#02
x 0.85, y 0.9, z 1.25
17 0.9 2 1.25 20.4
x 1.02
20 20
18 3 0.85 1.25 19.3
y 0.965
20 20
25 2 0.85 3 0.9 20.6
z 1.03
20 20
Iteration#03
x 1.02, y 0.965, z 1.03
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.26
17 0.965 2 1.03 20.025
x 1.00125
20 20
18 3 1.02 1.03 20.03
y 1.0015
20 20
25 2 1.02 3 0.965 20.065
z 1.00325
20 20
Gauss-Seidel Iteration Method
Question#3 Marks 10
Solve Question No. #2 by Gauss-Seidel iterative method and perform first three
iterations. What you see the difference after solving the same question by two different
iterative methods? Give your comments.
Solution:
The above system of linear equations is diagonally dominant; therefore, Gauss-Seidel
iterative method could be applied to find out real roots
20 x y 2 z 17
3x 20 y z 18
2 x 3 y 20 z 25
The above system of equations could be written in the form
17 y 2 z
x
20
18 3 x z
y
20
25 2 x 3 y
z
20
startingwith(0,0,0)
Iteration#01
17 0 2 0 17
x 0.85
20 20
18 3 0.85 0 20.55
y 1.0275
20 20
25 2 0.85 3 1.0275 20.2175
z 1.010875
20 20
Iteration#02
x 0.85, y 1.0275, z 1.010875
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.27
17 1.0275 2 1.010875 20.04925
x 1.0024625
20 20
18 3 1.0024625 1.010875 19.9965125
y 0.999825625
20 20
25 2 1.0024625 3 0.999825625 19.995598125
z 0.99977990625
20 20
Iteration#03
x 1.0024625, y 0.999825625, z 0.99977990625
17 0.999825625 2 0.99977990625 19.9904160625
x 0.9999708
20 20
18 3 0.9999708 0.99977990625 20.00013249375
y 1.0000066247
20 20
25 2 0.9999708 3 1.0000066247 20.000038526
z 1.0000019263
20 20
In Gauss Siedal Method, the newly computed values in each iteration are directly
involved to find the other value of the system of equation and save memory for
computation and hence results are more accurate.
==============================================================
Question 2
Do five iterations to solve the following system of equations by Gauss-Seidal iterative
method
10 x 2 y 3z 305
2 x 10 y 2 z 154
2 x y 10 z 120
Marks: 10
Solution:
Since, the given system is diagonally dominant; hence we can apply here the Gauss-
Seidal method.
From the given system of equations
1
x r 1 305 2 y r 3 z r
10
1
y r 1 154 2 x r 1 2 z r
10
1
z r 1 120 2 x r 1 y r 1
10
ITERATION 1
For r =0
Taking y=z=0 on right hand side of first equation. In second equation we take z=0 and
current value of x. In third equation we take current value of both x and y.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.28
1 305
x1 305 2(0) 3(0) 30.5
10 10
1 1 215
y1 154 2(30.5) 2(0) 154 61 21.5
10 10 10
1 1 202.5
z1 120 2(30.5) (21.5) 120 61 21.5 20.25
10 10 10
ITERATION 2
Similar procedure as used in Iteration 1 will be used for Iterations 2, 3, 4 and 5.
1 408.75
x2 305 2(21.5) 3(20.25) 40.875
10 10
1 276.25
y 2 154 2(40.875) 2(20.25) 27.625
10 10
1 229.375
z 2 120 2(40.875) (27.625) 22.938
10 10
ITERATION 3
1 429.064
x 3 305 2(27.625) 3(22.938) 42.906
10 10
1 285.688
y 3 154 2(42.906) 2(22.938) 28.569
10 10
1 234.381
z 3 120 2(42.906) (28.569) 23.438
10 10
ITERATION 4
1 432.452
x 4 305 2(28.569) 3(23.438) 43.245
10 10
1 287.366
y 4 154 2(43.245) 2(23.438) 28.737
10 10
1 235.227
z 4 120 2(43.245) (28.737) 23.523
10 10
ITERATION 5
1 433.043
x 5 305 2(28.737) 3(23.523) 43.304
10 10
1 287.654
y 5 154 2(43.304) 2(23.523) 28.765
10 10
1 235.373
z 5 120 2(43.304) (28.765) 23.537
10 10
Above Results are summarized in tabular form as
Variables
Iterations x y z
1 30.5 21.5 20.25
2 40.875 27.625 22.938
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.29
3 42.906 28.569 23.438
4 43.245 28.737 23.523
5 43.304 28.765 23.537
Hence the solution of the given system of equations, after five iterations, is
x 43.304
y 28.765
z 23.537
Relaxation Method
Matrix Inversion
Eigen Value Problems
Power Method
Question 1
Find the largest eigen value of the matrix
4 1 0
1 20 1 MARKS 10
0 0 4
And the corresponding eigenvector, by Power Method after fourth iteration starting with
the initial vector v(0) =(0,0,1)T
SOLUTION
4 1 0
Let A = 1 20 1
0 0 4
Choosing an initial vector as
v(0) =(0,0,1)T then
ITERATION 1
4 1 0 0 0
u (1)
[ A]v (0)
1 20 1 0 1
0 0 4 1 4
Now we normalize the resultant vector to get
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.30
0
1
u (1) 4 q1v (1)
4
1
Continuing this procedure for subsequent Iterations , we have
ITERATION 2
1
0
4 1 0 4
(1) 1
u (2) [ A]v 1 20 1 6
0 0 4 4 4
1
1
24
u (2) 6 1 q2 v (2)
2
3
ITERATION 3
1
4 1 0 24 1.167
u (3) [ A]v (2) 1 20 1 1 20.708
0 0 4 2 2.667
3
0.056
u (3) 20.708 1 q3v
(3)
0.129
ITERATION 4
4 1 0 0.056 1.224
u (4) [ A]v (3) 1 20 1 1 20.185
0 0 4 0.129 0.516
0.061
u 20.185 1
(4)
q4 v
(4)
0.026
Therefore, the largest eigen value and the corresponding eigen vector accurate to three
decimals places are
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.31
20.185
0.061
( X ) 1
0.026
=======================================================================
Question 1 marks 10
Find the largest eigen value of the matrix
6 5 2
3 1 4
1 10 3
And the corresponding eigen vector, by Power Method after fourth iteration starting with
the initial vector v(0) =(0,1,0)T
SOLUTION
6 5 2
Let A= 3 1 4
1 10 3
Choosing an initial vector
v(0) =(0,1,0)T then
ITERATION 1
6 5 2 0 5
u [ A]v 3 1 4 1 1
(1) (0)
1 10 3 0 10
Now we normalize the resultant vector to get
0.5
u (1) 10 0.1 q1v (1)
1
Continuing this procedure for subsequent Iterations, we have
ITERATION 2
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.32
6 5 2 0.5 5.5
u (2) [ A]v (1) 3 1 4 0.1 5.6
1 10 3 1 4.5
5.5
5.6
0.982
5.6
u (2) 5.6 1 q2 v
(2)
5.6
0.804
4.5
5.6
ITERATION 3
6 5 2 0.982 12.5
u [ A]v 3 1 4 1
(3) (2)
7.162
1 10 3 0.804 13.394
0.933
u 13.394 0.535 q3v (3)
(3)
1
ITERATION 4
6 5 2 0.933 10.273
u [ A]v 3 1 4 0.535 7.334
(4) (3)
1 10 3 1 9.283
1
u 10.273 0.714 q4 v (4)
(4)
0.904
Therefore, the largest eigen value and the corresponding eigen vector accurate to four
decimals places are
10.273
and
1
( X ) 0.714
0.904
Jacobi’s Method
Question 1
Using Jacobi’s method, find all the eigenvalues and the corresponding eigenvectors of the
following matrix,
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.33
1 2 2
2 1 2 Marks: 10
2 2 1
Note: Give results at the end of third rotation.
Solution.
1 2 2
Let A 2 1 2
2 2 1
1 2 2
Then, A 2 1 2
t
2 2 1
Hence,
Matrix A is real and Symmetric and Jacobi’s method can be applied.
Rotation 1
In Matrix A, all the off-diagonal elements are found be 2.
So, the largest off-diagonal element is found to be
a12 a13 a21 a23 a31 a32 2 .
Therefore, we can choose any one of them as the largest element.
Suppose, we choose a 12 as the largest element
Then, we compute the rotation angle as,
2a12 2 2
tan 2
a11 a22 1 1
Therefore
tan 2
2
2
π
θ=
4
Therefore, we construct an orthogonal matrix S1 such that,
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.34
cos 4 sin 4 0
S1 sin 4 cos 4 0
0 0 1
12 12 0
12 1
2
0
0 0 1
Now the first rotation gives,
D1 S11 AS1
12 1
2
0
1
Here, S1 = 12 1
2
0
0 0 1
So,
D1 S11 AS1
i.e.
12 1
2
0
1 2 2 2
1
1
2
0
2 1 2 1
D1 12 1
0 1
0
2 2 2
0 0 1
2 2 1 0 0 1
3 0 2.828
0 1 0
2.828 0 1
To check that we are right in our calculations, we can see that the sum of the diagonal
elements is, 3 1 1 3, which is same as the sum of the diagonal elements of the
original matrix A.
Rotation 2
For the second rotation we choose the largest off-diagonal element d13 d 31 2.828 , then
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.35
2d13 2 2.828 5.657
tan 2 2.828
d11 d33 3 1 2
So,
2 tan 1 2.828 70.526
35.263
Therefore, we construct an orthogonal matrix S 2 such that,
cos 35.263 0 sin 35.263
S2 0 1 0
sin 35.263 0 cos 35.263
0.817 0 0.577
0 1 0
0.577 0 0.817
0.817 0 0.577
S2 1 = 0 1 0
0.577 0 0.817
Now the rotation 2 gives,
D2 S21D1S2
0.817 0 0.577 3 0 2.828 0.817 0 0.577
D2 0 1 0 0 1 0 0 1 0
0.577 0 0.817 2.828 0 1 0.577 0 0.817
5 0 0
0 1 0
0 0 1
Again to check that we are right in our calculations, we can see that the sum of the
diagonal elements is 5 1 1 3 which is same as the sum of the diagonal
elements of the original matrix A.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.36
Rotation 3
We can see that in above iteration that D2 is a diagonal matrix, so we stop here and take
the Eigen Values as
λ1 = 5 , λ2 = -1 , λ3 = -1
Now the Eigenvectors are the columns vectors of the matrix S S1 S 2 , which are,
12 12 0 0.817 0 0.577
S 12 1
0 0 1 0
2
0
0 1 0.577 0 0.817
0.577 0.707 0.408
0.577 0.707 0.408
0.577 0 0.817
Therefore, the corresponding EigenVectors are
0.577 0.707 0.408
X 1 0.577 , X 2 0.707 , X 3 0.408
0.577 0 0.817
Interpolation
For a given table of values ( , ), 0,1,2,..., k k x y k= nwith equally spaced abscissas of a
function y= f(x),we define the forward difference operator Δ as follows,
These differences are called first differences of the function y and are denoted by the
symbol Δyi Here, Δ is called the first difference operator.
Similarly, rth Difference operator would be
Leading term =yo
Leading Diffenrence=∆y
Backward Difference Operators:
Central Difference is given by,
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.37
Shift operator, E
Let y = f (x) be a function of x, and let x takes the consecutive values x, x + h, x + 2h, etc.
We then define an operator having the property
E f( x) = f(x+h)
Thus, when E operates on f (x), the result is the next value of the function. Here, E is
called the shift operator. If we apply the operator E twice on f (x), we get
Thus, in general, if we apply the operator ‘E’ n times on f (x), we get
Enf(x)= f (x+nh)
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.38
Newton Forward Difference Interpolation.
Forward Differences
Question 2
Construct a forward difference table from the following values of x and y.
x 1.0 1.2 1.4 1.6 1.8 2.0 2.2
y=f(x) 2.0 3.5 -1.7 2.3 4.2 6.5 5.7
Marks:
10
Solution.
Forward-difference table
x y ∆y ∆2y ∆ 3y ∆ 4y ∆ 5y ∆ 6y
1.0 2
1.5
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.39
1.2 3.5 -6.7
-5.2 15.9
1.4 -1.7 9.2 -27.2
4 -11.3 41
1.6 2.3 -2.1 13.8 -60.8
1.9 2.5 -19.8
1.8 4.2 0.4 -6
2.3 -3.5
2.0 6.5 -3.1
-0.8
2.2 5.7
Backward Differences
Question 2
Construct a backward difference table from the following values of x and y
.
x 1.0 1.2 1.4 1.6 1.8 2.0 2.2
y=f(x) 2.0 3.5 -1.7 2.3 4.2 6.5 5.7
Marks:
10
Solution.
Backward-difference table
x y ∆y ∆2y ∆ 3y ∆ 4y ∆ 5y ∆ 6y
1.0 2
1.5
1.2 3.5 -6.7
-5.2 15.9
1.4 -1.7 9.2 -27.2
4 -11.3 41
1.6 2.3 -2.1 13.8 -60.8
1.9 2.5 -19.8
1.8 4.2 0.4 -6
2.3 -3.5
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.40
2.0 6.5 -3.1
-0.8
2.2 5.7
=====================================================================
Question 2 marks
10
Form a table of forward and backward differences of the function
f ( x) x 3 3x 2 5 x 7 For x 1, 0,1, 2,3, 4,5
SOLUTION
For the given function, the values of y for the given values of x are calculated as
for x 1
f (1) (1)3 3(1) 2 5(1) 7 6
for x 0
f (0) (0)3 3(0) 2 5(0) 7 7
for x 1
f (1) (1)3 3(1) 2 5(1) 7 14
for x 2
f (2) (2)3 3(2) 2 5(2) 7 21
for x 3
f (3) (3)3 3(3) 2 5(3) 7 22
for x 4
f (4) (4)3 3(4) 2 5(4) 7 11
for x 5
f (5) (5)3 3(5) 2 5(5) 7 18
So, the table of values of x and y is
x -1 0 1 2 3 4 5
y f ( x) -6 -7 -14 -21 -22 -11 18
Forward difference table
Forward difference table for the table of values of x and y f ( x) is shown below
x y ∆y ∆2y ∆ 3y ∆ 4y ∆ 5y ∆ 6y
-1 -6
-1
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.41
0 -7 -6
-7 6
1 -14 0 0
-7 6 0
2 -21 6 0 0
-1 6 0
3 -22 12 0
11 6
4 -11 18
29
5 18
Backward difference table
Backward difference table for the table of values of x and y f ( x) is shown below
x y y 2 y 3 y 4 y 5 y 6 y
-1 -6
-1
0 -7 -6
-7 6
1 -14 0 0
-7 6 0
2 -21 6 0 0
-1 6 0
3 -22 12 0
11 6
4 -11 18
29
5 18
Divided Differences
Question 1 marks
10
For the following table of values, estimate f (2.5) using Newton’s forward difference
interpolation formula.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.42
x 1 2 3 4 5 6 7 8
y f ( x) 1 8 27 64 125 216 343 512
SOLUTION
Forward difference table
Forward difference table for the given values of x and y is shown below
x y ∆y ∆2y ∆ 3y ∆ 4y ∆ 5y ∆ 6y
1 1
7
2 8 12
19 6
3 27 18 0
37 6 0
4 64 24 0 0
61 6 0
5 125 30 0
91 6
6 216 36
127
7 343 42
169
8 512
Newton’s forward difference interpolation formula is given by
p( p 1) 2 p( p 1)( p 2) 3
yx y0 py0 y0 y0
2! 3!
p( p 1)( p 2)( p 3) 4 p( p 1)( p 2)........( p n 1) n
y0 y0
4! n!
Here
x x0 2.5 1 1.5
p 1.5
h 1 1
And
y0 1, y0 7, 2 y0 12, 3 y0 6, 4 y0 5 y0 6 y0 0
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.43
So, by putting the above values in Newton’s forward difference interpolation formula,
We have
p( p 1) 2 p( p 1)( p 2) 3
yx y0 py0 y0 y0
2! 3!
1.5(1.5 1) 1.5(1.5 1)(1.5 2)
1 1.5(7) (12) (6)
2! 3!
1 10.5 4.5 0.375
15.625
i.e.
y2.5 15.625
===========================================================
Question 2 marks
10
Compute f (1.5) for the following data by using Newton’s divided difference interpolation
formula.
x 1 2 4 5 8
f ( x) 5 14 28 45 92
SOLUTION
The divided difference table for the given data is given by
x y 1st D.D. 2nd D.D. 3rd D.D. 4th D.D.
1 5
9
2 14 -2/3
7 1
4 28 10/3 -29/126
17 -11/18
5 45 -1/3
47/3
8 92
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.44
Newton ' s Divided Difference formula is
f ( x) y0 ( x x0 ) y[ x0 , x1 ] ( x x0 )( x x1 ) y[ x0 , x1 , x2 ] ( x x0 )( x x1 )( x x2 ) y[ x0 , x1, x2 , x3 ]
( x x0 )( x x1 )( x x2 )( x x3 ) y[ x0 , x1 , x2 , x3 , x4 ]
2
f (1.5) 5 + (1.5 - 1) 9 + (1.5-1)(1.5-2)(- ) + (1.5-1) (1.5-2)(1.5-4)1
3
29
+(1.5-1) (1.5-2)(1.5-4)(1.5-5)(- )
126
=5+4.5+0.1667+0.625+0.5035
=10.7952
Question 3 marks
10
Find y / (6) and y // (6) from the following table of values.
x 1 2 3 4 5 6
y 3 9 17 27 40 55
SOLUTION
Backward difference table
Backward difference table for the given values of x and y is shown below
x y y 2 y 3 y 4 y 5 y
1 3
6
2 9 2
8 0
3 17 2 1
10 1 -3
4 27 3 -2
13 -1
5 40 2
15
6 55
By backward difference formula for first derivative, we have
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.45
1 2 yn 3 yn 4 yn
yn yn
h 2 3 4
1 2 yn 3 yn 4 yn 5 yn
yn yn
h 2 3 4 5
1 2 1 2 3
y(6) 15
1 2 3 4 5
15 1 0.333 0.5 0.6
14.567
By backward difference formula for second derivative, we have
5 5
yn
1
h2
2
y n 3
y n
11 4
12
y n
6
yn
Therefore
1 11 5
y // (6) 2 [2 (1) (2) (3)]
1 12 6
2 1 1.833 2.5
3.333
Question 1
Marks 10
Find the interpolating polynomial for the following data by
a) Newton’s Forward Difference Formula
b)Lagrange’s Formula
Hence show that both the methods give raise to the same polynomial.
x 0 1 2 3
y 1 3 7 13
Solution:
By Newton’s Forward Difference Formula
The Newton’s Forward difference table is given as
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.46
x y y 2 y 3 y
0 1
2
1 3 2
4 0
2 7 2
6
3 13
We have Newton’s forward difference interpolation formula as
p( p 1) 2 p( p 1)( p 2) 3
y y0 py0 y0 y0 ... (1)
2! 3!
Since, the third and higher order differences are zero; Therefore, Newton’s forward difference
Interpolation formula reduces to
p( p 1) 2
y y0 py0 y0 (2)
2!
Here,
x0 =0, y0 =1, y0 =2, 2 y0 =2, h =1
And p is given by
x x0 x 0
p x
h 1
Substituting these values in Eq.(2), we have
x( x 1)
y 1 x (2) (2)
2
1 2 x x( x 1)
1 2x x2 x
x2 x 1
Hence, the required polynomial is
y x2 x 1.
Yahan Tak Midterm K Papers Khatm Ho Gai
b) By Lagrange’s Formula
Lagrange’s Interpolation formula is given by
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.47
( x x1 )( x x2 )( x x3 ) ( x x0 )( x x2 )( x x3 )
y y0 y1
( x0 x1 )( x0 x2 )( x0 x3 ) ( x1 x0 )( x1 x2 )( x1 x3 )
( x x0 )( x x1 )( x x3 ) ( x x0 )( x x1 )( x x2 )
y2 y3
( x2 x0 )( x2 x1 )( x2 x3 ) ( x3 x0 )( x3 x1 )( x3 x2 )
By putting the values in the formula, we get
( x 1)( x 2)( x 3) ( x 0)( x 2)( x 3)
y (1) (3)
(0 1)(0 2)(0 3) (1 0)(1 2)(1 3)
( x 0)( x 1)( x 3) ( x 0)( x 1)( x 2)
(7) (13)
(2 0)(2 1)(2 3) (3 0)(3 1)(3 2)
x3 6 x 2 11x 6 3x 3 15 x 2 18 x x 3 6 x 2 11x 6 3 x 3 15 x 2 18 x
y y
(1)(2)(3) (1)(1)(2) 6 2
7( x3 4 x 2 3 x) 13( x 3 3 x 2 2 x) 7 x 28 x 21x 13 x 39 x 2 26 x
3 2 3
+ +
(2)(1)(1) (3)(2)(1) 2 6
x 3 6 x 2 11x 6 3x 3 15 x 2 18 x 7 x 3 28 x 2 21x 13 x 3 39 x 2 26 x
y
6 2 2 6
x 6 x 11x 6 9 x 45 x 54 x 21x 84 x 63x 13x 39 x 2 26 x
3 2 3 2 3 2 3
6
x3 9 x3 21x 3 13x 3 6 x 2 45 x 2 84 x 2 39 x 2 11x 54 x 63x 26 x 6
y
6
6x 6x 6
2
6
6( x x 1)
2
6
x x 1
2
Thus we have the equation
y x2 x 1
This is the same polynomial as obtained in Newton’s Forward difference interpolating
formula. Hence it is proved that both the methods give rise to the same polynomial
Langrange’s Interpolation
Differentiation Using Difference Operators
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.48
Numerical Integration
Trapezoidal Rule
Simpson’s 1/3 and 3/8 rules
Differential Equations
Taylor Series Method
Euler Method
Runge-Kutta Method
Milne’s Predictor Corrector Method
Adam Moultan’s Predictor Corrector Method
FAQ updated version.
Question: What is Bracketing method?
Answer: Methods such as bisection method and the false position method of finding
roots of a nonlinear equation f(x) = 0 require bracketing of the root by two guesses. Such
methods are called bracketing methods. These methods are always convergent since they
are based on reducing the interval between the two guesses to zero in on the root.
Question: What is an Open method?
Answer: In the Newton-Raphson method, the root is not bracketed. Only one initial
guess of the root is needed to get the iterative process started to find the root of an
equation. Hence, the method falls in the category of open methods.
Question: Explain Muller’s method briefly.
Answer: In Muller’s method, f (x) = 0 is approximated by a second degree polynomial;
that is by a quadratic equation that fits through three points in the vicinity of a root. The
roots of this quadratic equation are then approximated to the roots of the equation f (x)
0.This method is iterative in nature and does not require the evaluation of derivatives as
in Newton-Raphson method. This method can also be used to determine both real and
complex roots of f (x) = 0.
Question: Explain the difference between the linear and non-linear equations.
Answer: Linear Equation An algebraic equation is said to be linear in which each
term is either a constant or the product of a constant and the first power of a single
variable. One or more variables can be involved in the linear equations. e.g. x+3y+z=0 2x-
y+4z=7 etc. Non-Linear Equation An equation is said to Non-Linear equation if it is not
linear. Equations involving the power of the variable 2 or higher, transcendental,
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.49
logarithmic and trigonometric equations etc lie in the category of Non-Linear equations.
e.g. x2+5x+3=0 sinx+3y+9=0 xlogx-7x+4y=2 etc.
Question: Explain which value is to be choosed as X0 in N-R method.
Answer: If, for a given function, f(a)*f(b)<0, then any value between a and b inclusive
can be chosen as x0.
Question: Define iterative method of solving linear equations with two examples.
Answer: Under iterative methods, the initial approximate solution is assumed to be
known and is improved towards the exact solution in an iterative way. We consider
Jacobi, Gauss-Seidel and relaxation methods under iterative methods.
Question: Define Pivoting.
Answer: The Gaussian elimination method fails if any one of the pivot elements
becomes zero. In such a situation, we rewrite the equations in a different order to avoid
zero pivots. Changing the order of equations is called pivoting.
Question: Write the two steps of solving the linear equations using Gaussian
Elimination method.
Answer: In this method, the solution to the system of equations is obtained in two
stages.
i) the given system of equations is reduced to an equivalent upper triangular form
using elementary transformations
ii) the upper triangular system is solved using back substitution procedure
Question: Describe Gauss-Jordan elimination method briefly
Answer: This method is a variation of Gaussian elimination method. In this method,
the elements above and below the diagonal are simultaneously made zero. That is a given
system is reduced to an equivalent diagonal form using elementary transformations.
Then the solution of the resulting diagonal system is obtained. Sometimes, we normalize
the pivot row with respect to the pivot element, before elimination. Partial pivoting is also
used whenever the pivot element becomes zero.
Question: Describe briefly Crout’s reduction method.
Answer: Here the coefficient matrix [A] of the system of equations is decomposed into
the product of two matrices [L] and [U], where [L] is a lower-triangular matrix and [U] is
an upper-triangular matrix with 1’s on its main diagonal.
Question: Describe briefly the Jacobi’s method of solving linear equations.
Answer: This is an iterative method, where initial approximate solution to a given
system of equations is assumed and is improved towards the exact solution in an
iterative way.
Question: What is the difference between Jacobi’s method and Gauss Seidal
method?
Answer: The difference between jacobi’s method and gauss Seidel method is that in
jacobi’s method the approximation calculated are used in the next iteration for next
approximation but in Gauss-seidel method the new approximation calculated is instantly
replaced by the previous one.
Question: What is the basic idea of Relaxation method?
Answer: We can improve the solution vector successively by reducing the largest
residual to zero at that iteration. This is the basic idea of relaxation method.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.50
Question: How the fast convergence in the relaxation method is achieved?
Answer: To achieve the fast convergence of the procedure, we take all terms to one
side and then reorder the equations so that the largest negative coefficients in the
equations appear on the diagonal.
Question: Which matrix will have an inverse?
Answer: Every square non-singular matrix will have an inverse.
Question: What are the popular methods available for finding the inverse of a
matrix?
Answer: Gauss elimination and Gauss-Jordan methods are popular among many
methods available for finding the inverse of a matrix.
Question: Explain Gaussian Elimination Method for finding the inverse of a
matrix.
Answer: In this method, if A is a given matrix, for which we have to find the inverse;
at first, we place an identity matrix, whose order is same as that of A, adjacent to A
which we call an augmented matrix. Then the inverse of A is computed in two stages. In
the first stage, A is converted into an upper triangular form, using Gaussian elimination
method In the second stage, the above upper triangular matrix is reduced to an identity
matrix by row transformations. All these operations are also performed on the adjacently
placed identity matrix. Finally, when A is transformed into an identity matrix, the
adjacent matrix gives the inverse of A. In order to increase the accuracy of the result, it is
essential to employ partial pivoting
Question: What are the steps for finding the largest eigen value by power method.
Answer:
Procedure Step 1: Choose the initial vector such that the largest element is unity.
Step 2: The normalized vector is pre-multiplied by the matrix [A].
Step 3: The resultant vector is again normalized.
Question: What is the method for finding the eigen value of the least magnitude
of the matrix [A]?
Answer: For finding the eigen value of the least magnitude of the matrix [A], we have
to apply power method to the inverse of [A].
Question: What is interpolation?
Answer: The process of estimating the value of y, for any intermediate value of x,
lying inside the table of values of x is called interpolation.
Question: What is extrapolation?
Answer: The method of computing the value of y, for a given value of x, lying outside
the table of values of x is known as extrapolation.
Question: What happens when shift operator E operates on the function.
Answer: When shift operator E operates on the function it results in the next value of
the function.Ef(x)=f(x+h)
Question: What is the basic condition for the data to apply Newton’s
interpolation methods?
Answer: To apply Newton’s interpolation methods, data should be equally spaced.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.51
Question: When is the Newton’s forward difference interpolation formula used?
Answer: Newton’s forward difference interpolation formula is mainly used for
interpolating the values of y near the beginning of a set of tabular values and for
extrapolating values of y, a short distance backward from y0.
Question: When is the Newton’s forward difference interpolation formula used?
Answer: Newton’s forward difference interpolation formula is mainly used for
interpolating the values of y near the beginning of a set of tabular values and for
extrapolating values of y, a short distance backward from y0.
Question: When the Newton’s backward difference interpolation formula is used?
Answer: For interpolating the value of the function y = f (x) near the end of table of
values, and to extrapolate value of the function a short distance forward from yn,
Newton’s backward Interpolation formula is used.
Question: What is the formula for finding the value of p in Newton’s forward
difference interpolation formula?
Answer: P=(x-x0)/h
Question: What is the formula for finding the value of p in Newton’s backward
difference interpolation formula?
Answer: P=(x-xn)/h
Question: If the values of the independent variable are not equally spaced then
which formula should be used for interpolation?
Answer: If the values of the independent variable are not given at equidistant
intervals, then the Lagrange’s interpolation formula should be used.
Question: To use Newton’s divided difference interpolation formula, what should
the values of independent variables be?
Answer: To use Newton’s divided difference interpolation formula, the values of
independent variables should not be equally spaced.
Question: Which difference formula is symmetric function of its arguments?
Answer: Divided difference formula is symmetric function of its arguments.
Question: Is the interpolating polynomial found by Lagrange’s and Newton’s
divided difference formulae is same?
Answer: Yes. The interpolating polynomial found by Lagrange’s and Newton’s divided
difference formulae is one and the same.
Question: Which formula involves less number of arithmetic operations? Newton
or Lagrange’s?
Answer: Newton’s formula involves less number of arithmetic operations than that of
Lagrange’s.
Question: When do we need Numerical Methods for differentiation and
integration?
Answer: If the function is known and simple, we can easily obtain its derivative (s) or
can evaluate its definite integral However, if we do not know the function as such or the
function is complicated and is given in a tabular form at a set of points x0,x1,…,xn, we
use only numerical methods for differentiation or integration of the given function.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.52
Question: If the value of the independent variable at which the derivative is to be
found appears at the beginning of the table of values, then which formula should be
used?
Answer: If the value of the independent variable at which the derivative is to be
found appears at and near the beginning of the table, it is appropriate to use formulae
based on forward differences to find the derivatives.
Question: If the value of the independent variable at which the derivative is to be
found occurs at the end of the table of values, then which formula should be used?
Answer: If the value of the independent variable at which the derivative is to be
found occurs at the end of the table of values, it is appropriate to use formulae based on
backward differences to find the derivatives.
Question: Why we need to use RICHARDSON’S EXTRAPOLATION METHOD?
Answer: To improve the accuracy of the derivative of a function, which is computed
by starting with an arbitrarily selected value of h, Richardson’s extrapolation method is
often employed in practice.
Question: To apply Simpson’s 1/3 rule, what should the number of intervals be?
Answer: To apply Simpson’s 1/3 rule, the number of intervals must be even.
Question: To apply Simpson’s 3/8 rule, what should the number of intervals be?
Answer: To apply Simpson’s 3/8 rule, the number of intervals must be multiple of 3.
Question: What is the order of global error in Simpson’s 1/3 rule?
Answer: The global error in Simpson’s 1/3 rule is of the order of 0(h4).
Question: Is the order of global error in Simpson’s 1/3 rule equal to the order of
global error in Simpson’s 3/8 rule?
Answer: Yes. The order of global error in Simpson’s 1/3 rule equal to the order of
global error in Simpson’s 3/8 rule.
Question: What is the order of global error in Trapezoidal rule?
Answer: The global error in Trapezoidal rule is of the order of 0(h2).
Question: What is the formula for finding the width of the interval?
Answer: Width of the interval, h, is found by the formula h=(b-a)/n
Question: What type of region does the double integration give?
Answer: Double integration gives the area of the rectangular region.
Question: Compare the accuracy of Romberg’s integration method to trapezoidal
and Simpson’s rule.
Answer: Romberg’s integration method is more accurate to trapezoidal and
Simpson’s rule
Question: What is the order of global error in Simpson’s 3/8 rule?
Answer: The global error in Simpson’s 3/8 rule is of the order of 0(h4).
Question: Which equation models the rate of change of any quantity with respect
to another?
Answer: An ordinary differential equation models the rate of change of any quantity
with respect to another.
Question: By employing which formula, Adam-Moulton P-C method is derived?
Answer: Adam-Moulton P-C method is derived by employing Newton’s backward
difference interpolation formula.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.53
Question: What are the commonly used number systems in computers?
Answer: Binary Octal Decimal Hexadecimal
Question: If a system has the base M, then how many different symbols are
needed to represent an arbitrary number? Also name those symbols.
Answer: M 0, 1, 2, 3, … , M-1
Question: What is inherent error and give its example.
Answer: It is that quantity of error which is present in the statement of the problem
itself, before finding its solution. It arises due to the simplified assumptions made in the
mathematical modeling of a problem. It can also arise when the data is obtained from
certain physical measurements of the parameters of the problem. e.g. if writing 8375
instead of 8379 in a statement lies in the category of inherent error.
Question: What is local round-off error?
Answer: At the end of computation of a particular problem, the final results in the
computer, which is obviously in binary form, should be converted into decimal form-a
form understandable to the user-before their print out. Therefore, an additional error is
committed at this stage too. This error is called local round-off error.
Question: What is meant by local truncation error?
Answer: Retaining the first few terms of the series up to some fixed terms and
truncating the remaining terms arise to local truncation error.
Question: What is transcendental equation and give two examples.
Answer: An equation is said to be transcendental equation if it has logarithmic,
trigonometric and exponential function or combination of all these three. For example it
is a transcendental equation as it has an exponential function These all are the examples
of transcendental equation.
Question: What is meant by intermediate value property?
Answer: If f(x) is a real valued continuous function in the closed interval if f(a) and
f(b) have opposite signs once; that is f(x)=0 has at least one root such that Simply If f(x)=0
is a polynomial equation and if f(a) and f(b) are of different signs ,then f(x)=0 must have
at least one real root between a and b.
Question: What is direct methods of solving equations?
Answer: Those methods which do not require any information about the initial
approximation of root to start the solution are known as direct methods. The examples of
direct methods are Graefee root squaring method, Gauss elimination method and Gauss
Jordan method. All these methods do not require any type of initial approximation.
Question: What is Iterative method of solving equations?
Answer: These methods require an initial approximation to start. Bisection method,
Newton raphson method, secant method, jacobi method are all examples of iterative
methods
Question: If an equation is a transcendental, then in which mode the calculations
should be done?
Answer: All the calculations in the transcendental equations should be done in the
radian mode.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.54
Question: What is the convergence criterion in method of iteration?
Answer: If be a root of f(x) =0 which is equivalent to I be any interval containing the
point x= and will converge to the root provided that the initial approximation is chosen in
I
Question: When we stop doing iterations when Toll is given?
Answer: Here if TOL is given then we can simply find the value of TOL by subtracting
both the consecutive roots and write it in the exponential notation if the required TOL is
obtained then we stop.
Question: How the value of h is calculated in interpolation?
Answer: There are two types of data in the interpolation one is equally spaced and
other is unequally spaced. In equally spaced data we need to calculate the value of h that
is calculated by subtracting any two consecutive values and taking their absolute value.
Question: What is an algebraic equation?
Answer: An algebraic equation is an equation which is purely polynomial in any
variable. Supposed x2+3x+2=0, x4+3x2=0, y3+6y2=0 all are algebraic equations as these
are purely polynomial in x and y variable.
Question: What is Descartes rule of signs?
Answer: The number of positive roots of an algebraic equation f(x) =0 can not exceed
the no of changes in signs. Similarly the no of negative roots of negative roots of and
algebraic equation can not exceed the no of changes in sign of equation f (-x) =0.
Question: What are direct methods?
Answer: The numerical methods which need no information about the initial
approximation are known as direct methods like Graffee’s root squaring method.
Question: What is meant by iterative methods?
Answer: The methods which need one or more iterations are known as iterative
methods like bisection method, Newton raphson method, and many other methods.
Question: What is graphically meant by the root of the equation?
Answer: If the graph of a function f(x) =0 cuts the x-axis at a point a then a is known
as the root of the equation.
Question: Q. What is the difference between open and bracketing method?
Answer: In open methods we need only one initial approximation of the root that may
be any where lying and if it is not very close then we have to perform more iteration and
the example of open method is Newton raphson method. In bracketing method we
bracket the root and find that interval in which root lies means we need two initial
approximations for the root finding. Bisection method is an example of the bracketing
method.
Question: Condition for the existence of solution of the system of equations.
Answer: If the |A| is not equal to zero then the system will have a unique solution if
|A|=0 then the system will have no solution
Question: Should the system be diagonally dominant for gauss elimination
method?
Answer: The system of equation need not to be diagonally dominant for Gauss
elimination method and gauss Jordan method for both the direct method it is not
necessary for the system to be diagonally dominant .it should be diagonally dominant for
iterative methods like jacobie and gauss seidel method.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.55
Question: What is meant by diagonally dominant system?
Answer: A system a1x+b1y+c1z=d1, a2x+b2y+c2z=d2, a3x+b3y+c3z=d3 is said to be
diagonally dominant if the following condition holds. |a1| => (greater or equal)
|b1|+|c1| |b2| => (greater or equal) |a2|+|c2| |c3| => (greater or equal) |a3|+|b3|
Question: State the sufficient condition for the convergence of the system of
equation by iterative methods.
Answer: A sufficient condition for convergence of iterative solution to exact solution
is |a1| => (greater or equal) |b1|+|c1| |b2| => (greater or equal) |a2|+|c2| |c3| =>
(greater or equal) |a3|+|b3| For the system a1x+b1y+c1z=d1, a2x+b2y+c2z=d2,
a3x+b3y+c3z=d3 Similarly for the system with more variables we can also construct the
same condition
Question: The calculation for numerical analysis should be done in degree or
radians.
Answer: All the calculation for numerical analysis should be done in radians not in
degrees set your calculator in radians mode and suppose the value of pi=3.14.
Question: How we can identify that Newton forward or backwards interpolation
formula is to be used.
Answer: If the value at which we have to interpolate is in the start of the table then
we will use Newton’s forward interpolation formula if it is at the end of the table then we
will use the Newton’s backward interpolation formula.
Question: What is meant precision and accuracy?
Answer: Precision and accuracy are two terms which are used in numerical
calculations by precision we mean that how the values in different iterations agree to
each other or how close are the different values in successive iterations. For example you
have performed 3 different iterations and result of all the iteration are
1.32514,1.32516,31.32518 these three values are very precise as these values agree with
each other . Accuracy means the closeness to the actual value. Suppose that you have
calculated an answer after some iteration and the answer is 2.718245 and the actual
answer is 2.718254 and the answer calculated is very accurate but if this answer is
2.72125 then it is not accurate.
Question: What is the condition that a root will lie in an interval.
Answer: Suppose that you have a function f(x) and an interval [a,b]and you calculate
both f(a)and f(b)if f(a)f(b)<0 then there must exist a root for this function in this
interval.In simple words f(a) and f(b)must have opposite signs.
Question: How the divided difference table is constructed?
Answer: x y First difference Second difference Third difference 1 0.2 y0
y1-y0/x1-x0=sy0 sy1- sy0/x2-x0=s2y0 s2y0-s2y1/x3-x0=s3y0 s3y1-
s3y0/x4x0 =s4y0 2 0.69 y1 y2-y1/x2-x1=sy1 sy2- sy1/x3-x1=s2y1 s2y2-
s2y1/x4-x1=s3y1 3 0.985 y2 y3-y2/x3-x2=sy2 sy3- sy2/x4-x1=s2y2 4
1.2365 y3 y4-y3/x4-x3=sy3 5 2.3651 y4 This is the complete table
showing the central differences . Here s stands for difference
Question: What is Gauss-seidel Method
Answer: It is also an iterative method and in this we also check either the system is
diagonally dominant and, if so then we proceed on the same lines as in the Jacobi’s
method. The difference in both is when we calculate a value of any variable that is
instantly replaced by the previous value. If we have initial values like x1=0,x2=0 and
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.56
x3=0 in the case of three variables in the first iteration we will use all the three values to
calculate x1 but to calculate the value of x2 we will use the recent value of x1 which is
calculated in the first iteration.
Question: What is partial and full pivoting?
Answer: Partial and full pivoting, In gauss elimination method when you have any of
diagonal element aii zero it means the solution does not exist to avoid this we change the
equation so that a non zero pivot is achieved .Now you have an argument matrix in
which he elements in the first column are 1,3,4 respectively here in case of partial
pivoting we will replace first element with the last element it is done by replacing the first
equation with last equation it is known as partial pivoting. In full pivoting we change
rows and columns but that is not implemented manually it is used in computers.
Question: How the initial vector is choose in power method?
Answer: Choose the initial vector such that the largest element is unity. But no
problem with this you will precede with this vector as it is and will not change it. In the
given question the initial vector is provided so you have to proceed with the given vector
and will normalize this vector and normalizing this vector you will keep in mind that the
greatest value should be unity. An eigenvector Vis said to be normalized if the coordinate
of largest magnitude is equal to unity Actually this condition is for the normalization of
the vectors when you have to normalize the vector you keep in mind that the larget
entery in the vector must be 1 and so you take the largest element and divide the
remaining values by the greatest value.
Question: what is the relation ship between p=0 and non zero p in interpolation.
Answer: fp=fo+pDfo+1/2(p2-p)D2fo+1/6(p3-3p2+2p)D3fo+1/24(p4-6p3+11p2-
6p)D4fo+… This is the interpolation formula For the derivative you will have to take the
derivative of this formula w.r.t p and you will get f’p=1/h{Dfo+1/2(2p-1)D2fo+1/6(3p2-
6p+2)D3fo+…} Now put here p=0 f’p=1/h{Dfo+1/2(0-1)D2fo+1/6(0-0+2)D3fo+…} So the
formula becomes f’p=1/h{Dfo+D2fo/2+D3fo/3+…} so this is the relation ship between
non zero and zero p no when you have to calculate p the you use formula p=x-xo/h so
this is the impact of x.
Question: What is chopping and Rounding off?
Answer: Chopping and rounding are two different techniques used to truncate the
terms needed according to your accuracy needs. In chopping you simply use the
mentioned number of digits after the decimal and discard all the remaining terms.
Explanation (1/3 - 3/11) + 3/20=(0.333333…-0.27272727..)+0.15=(0.333-0.272)+0.15
This is the three digit chopping.
Question: When the forward and backward interpolation formulae are used?
Answer: In interpolation if we have at the start then we use the forward difference
formula and the formula to calculate p is x-x0/h. If the value of x lies at the end then we
use Newton’s backward formula and formula to calculate the value of p is x-xn/h. Now I
come to your question as in this case the value lies at the end so 6 will be used as the xn.
This procedure has been followed by the teacher in the lectures. But some authors also
use another technique that is if you calculate the value of p and that is negative then the
origin is shifted to that value for which the value of p becomes positive. And then
according to that origin the values of differences are used and you need not follow this
procedure.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.57
Question: What is forward and backward difference operator and the construction
of their table.
Answer: For forward Dfr =fr+1 –fr Df0 = f1-f0 In terms of y Dyr+1=yr+1-yr D stands
for the forward difference operator For backward Dfr =fr –fr-1 Df1 = f1-f0 In terms of y
Dy1=y1-y0 Here D stands for backwards operator Now the construction of the difference
table is based on X Y 1st forward 2nd forward 3rd forward x1 Y1 Y2-
Y1=Dy0 x2 Y2 Y3-Y2=Dy1 x3 Y3 Y4-Y3=Dy2 x4 Y4 Now consider
the construction of table for the backward table X Y 1st forward 2nd forward 3rd
forward x1 Y1 Y2-Y1=Dy1 x2 Y2 Y3-Y2=Dy2 x3 Y3 Y4-Y3=Dy3 x4 Y4
Dear student this is the main difference in the construction of the forward and
backwards difference table when you proceed for forward difference table you get in first
difference the value Dy0 but in the construction of backwards difference table in the first
difference you get Dy1 and in the second difference in the forward difference table you get
D2 y0 and in the backward difference table the first value in the second difference is D2
y1. I think so you have made it clear.
Question: What is Jacobi’s method?
Answer: Jacobi’s Method It is an iterative method and in this method we first of all
check either the system is diagonally dominant and, if the system is diagonally dominant
then we will calculate the value of first variable from first equation in the form of other
variables and from the second equation the value of second variable in the form of other
variables and so on. We are provided with the initial approximations and these
approximations are used in the first iteration to calculate first approximation of all the
variables. The approximations calculated in the first iteration are used in the second
iteration to calculate the second approximations and so on.
Question: what is Simpson's 3/8th rule.
Answer: The general formula for Simpson’s 3/8th rule is
3h/8[f0+3(f1+f2)+2f3+(3f4+f5)+2f6+…+3(fn-2+fn-1)+fn] Now if we have to calculate the
integral by using this rule then we can simply proceed just write first and last vale and
distribute all the remaining values with prefix 3 and 2 Like you have f0,f1,f2,f3,f4 Then
the integral can be calculated as 3h/8[f0+f4+3(f1+f2)+2f4] If we have values like
f0,f1,f2,f3,f4,f5 Then integral can be calculated as 3h/8[f0+3(f1+f2)+2f3+3f4+f5] Similarly
proceeding in this fashion we can calculate the integral in this fashion
Question: what is classic runge-kutta method
Answer: The fourth order Runge-Kutta method is known as the classic formula of
classic Runge-Kutta method. yn+1=yn+1/6(k1+2k2+2k3+k4) Where k1=hf(xn,yn)
k2=hf(xn+h/2+yn+k1/2) k3=hf(xn+h/2,yn+k’2/2) k4=hf(xn+h,yn+k3)
Question: What is meant by TOL?
Answer: The TOL means the extent of accuracy which is needed for the solution. If
you need the accuracy to two places of decimal then the TOL will be 10-2 .similarly the
10-3 means that the accuracy needed to three places of decimal. Suppose you have the
root from last iteration 0.8324516 and 0.8324525 if we subtract both and consider
absolute value of the difference 0.0000009 now it can be written as 0.09*10-5 so the TOL
in this case is 10-5.similarly if we have been provided that you have to for the TOL 10-2
you will check in the same way. In the given equation you will solve the equation by any
method and will consider some specific TOL and try to go to that TOL. Some time no TOL
is provided and you are asked to perform to some specific no of iterations.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.58
Question: what is meant by uniqueness of LU method.
Answer: An invertible (whose inverse exists) matrix can have LU factorization if and
only if all its principal minors(the determinant of a smaller matrix in a matrix) are non
zero .The factorization is unique if we require that the diagonal of L or U must have
1’s.the matrix has a unique LDU factorization under these condition . If the matrix is
singular (inverse does not exists) then an LU factorization may still exist, a square matrix
of rank (the rank of a matrix in a field is the maximal no of rows or column) k has an LU
factorization if the first k principal minors are non zero. These are the conditions for the
uniqueness of the LU decomposition.
Question: how the valu of h is calculated from equally spaced data.
Answer: Consider the following data x y 1 1.6543 2 1.6984 3 2.4546 4 2.9732 5
3.2564 6 3.8765 Here for h=2-1=3-2=1 x y 0.1 1.6543 0.2 1.6984 0.3 2.4546 0.4 2.9732
0.5 3.2564 0.6 3.8765 Here for the calculation of h =0.2-0.1=0.3-0.2=0.1 I think so that
you may be able to understand .
Glossary (Updated Version)
Absolute Error :
The absolute error is used to denote the actual value of a quantity less it’s rounded
value if x and x* are respectively the rounded and actual values of a quantity, then
absolute error is defined by AE=|x-x*|
Accuracy :
The Extent of the closeness between the actual value and estimated value is
known as accuracy. Suppose you have taken readings like 2.1234,2.1236and 2.1238 and
2.45,2.52,2.63 Now if the actual root is 2.65 then the second values are more accurate
and not precise but first set of values is precise but not accurate as it differs from the
actual value.
Algebraic equation :
The equation f(x)=0 is known as algebraic equation if it is purely a polynomial in x.
Like f(x)=3 ,f(x)=x3-3x-4 are some examples of algebraic equations and the if variable is
changed to y x or any one then under the same conditions it will be algebraic.
Bisection method :
It is a bracketing method which is used to locate the root of an equation and it
make use of intermediate value property to locate the interval in which the root lies,
formula which is used is (a+b)/2 where a and b are points such that f(a)f(b)<0.
Bracketing Method :
The iterative methods which require two initial approximation for it’s first iteration
are known as bracketing methods. Bisection method is the example of the bracketing
method as it requires a interval for the approximation of the root.
Crout’s Method :
This method is used to solve the system of the equations by decomposing the
system into two matrices L and U where L is a lower triangular matrix and U is an upper
triangular matrix. In U all the elements in the main diagonal are 1.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.59
Descartes rule of signs : The no of positive roots of an algebraic equation f(x)=0 with real
coefficients cannot exceed the no of changes in the in sign of the coefficients in the
polynomial f(x)=0,similarly the no of negative roots of f(x)=0 can not exceed the no of
changes in the sign of the coefficients of f(-x)=0
Direct methods :
These are the methods which do not need the knowledge of the initial
approximation are known as direct methods.
Gauss Elimination Method :
It is a direct method which is used to solve a system of equations
Gauss Seidel iterative method :
It is an iterative method in which an initial approximation is given. First of all
system should be checked either it is diagonally dominant or not if it is not then it is
made diagonally dominant. Secondly the first variable is calculated in terms of other
variables from first equation and Second variable from second equation and so on .the
previous value of the variable is replaced by new value instantly as it is obtained. This is
the difference in the gauss seidel and jacobie iterative method.
Graeffee’s root squaring method :
This method is used to find all the roots of the polynomial equation.
Intermediate value property :
If for an equation f(x)=0 for two values a and b we have such that f(a)f(b)<0 then
there must exist a root between a and b in the interval [a,b].
Inverse of a matrix :
The matrix B is said to be inverse of a matrix A if the product of and B is I identity
matrix.
Iterative methods :
Iterative methods are those type of methods which always require an initial
approximation to start an iteration
Jacobie’s iterative method: :
It is an iterative method in which an initial approximation is given. First of all
system should be checked either it is diagonally dominant or not if it is not then it is
made diagonally dominant. Secondly the first variable is calculated in terms of other
variables from first equation and Second variable from second equation and so on after
first iteration in the second iteration all the variables are replaced by the previous one.
Muller’s Method :
In Muller’s method f(x)=0 is approximated by a second degree polynomial, that is
by quadratic equation that fits through three points in the vicinity of a root .Then roots of
this quadratic equation are then approximated to the roots of equation f(x)=0
Newton Raphson Method. :
It is an open method which is used to locate a root of the equation, it needs only
one initial approximation for it’s first iteration. xn+1=xn-(fxn/f’xn)
Non singular matrix :
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.60
A matrix is said to be non singular if the determinant of the matrix is non zero.
|A| is not equal to zero and the inverse of non singular exists.
Open methods :
The methods which require only one initial approximation to start the first
iteration, for example the Newton’s Raphsom method is known as the open method as it
requires only one initial approximation.
Pivoting :
If any of the pivot elements in gauss elimination become zero then this method
fails so to Avoid such type of situation equation are rearranged to get rid of zero pivot
element, this procedure is known as pivoting.
Polynomial:
An expression of the form where n is a
positive integer and are real constants, such type of expression is called
an nth degree polynomial in x if ao ≠ 0
Precision :
The extent of closeness of different measurement taken to estimate a certain
value. Suppose you have done different iterations to measure the root of an equation and
take different values as 2.1234, 2.1236, and 2.1238 these all the values are very close to
each other so these are very precise.
Regula –Falsi method :
It is also an iterative method and is a bracketing method and use intermediate
value property to get it’s initial guess. The formula used for this is xn+1=xn-(xn-xn-
1/f(xn)-f(xn-1))
Relative Error : It is the ratio of the absolute error to the actual value of the quantity.
Thus RE=AE=AE\|x*|
Relaxation method :
It is also an iterative method and in this method you solve the system of equation
by making the greatest residual to zero.
Root : If you have an equation f(x)=0 then the no a is said to be the root of the equation if
f(a)=0 Suppose you have an equation f(x)=x2-4 the 2 is a root of the equation as f(2)=4-
4=0
Secant Method :
The secant method is also an open method and it takes two initial values for it’s
first approximation, the formula used for this is known as xn+1={xn-1f(xn)-xnf(xn-
1)}/{f(xn)-f(xn-1)}
Significant digits :
A significant digit in an approximate no is a digit, which gives reliable information
about the size of the number. In other words a significant digit is used to express
accuracy, that Is how many digits in the no have meaning.
Singular matrix :
A matrix s said to be a singular matrix if the determinant of the matrix is zero
|A|=0,the inverse of the singular matrix do not exist.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.61
Transcendental equation :
An equation f(x) =0 is said to be transcendental equation if it contains
trigonometric, logarithmic and exponential functions
Truncation Error :
It is defined as the replacement of one series by another series with fewer terms.
The error arising by this approximation is known as truncation error.
Important Formula For MTH603
Bisection Method
x0 x1
x2
2
Muller Method
x0 2c
x
b b 2 4ac
in this formula x0,x1, x2 will be and U will just put the values according to the given
formulas.
h2 f 1 (h1 h 2) f 0 h1 f 2
a
h1h 2(h1 h 2)
f 1 f 0 ah12
b
h1
c f ( x0)
h1 x1 x0
h2 x0 x 2
Regula Falsi Method (Method of False position)
x 2 x1
x3 x 2 f ( x 2)
f ( x 2) f ( x1)
Newton Rophson method
f ( x0)
x1 x0
f ' ( x0)
Secant Method
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.62
x0 f ( x1) x1 f ( x0)
x2
f ( x1) f ( x0)
Newton’s Formula
1
x1 ( x0 n )
2 x0
in this formula
x0= 2 perfect square near to 12 such like 9 and 16
Graffee root squaring method
f ( x) f ( x) a32 x 6 (a2 2 2a1a3 ) x 4 (a12 2a0 a2 ) x 2 x02
the truncation error (TE) is given by
The TE is independent of the computer used.
===================================================================
Spring 2011 Latest Papers (Current)
Papers Number 01
Total 26 questions
20 mcqs
Question: What are the steps for finding the largest eigen value by power method?
Answer: Procedure
Step 1: Choose the initial vector such that the largest element is unity.
Step 2: The normalized vector is pre-multiplied by the matrix [A].
Step 3: The resultant vector is again normalized.
Step 4: This process of iteration is continued and the new normalized vector is
repeatedly pre-multiplied by the matrix [A] until the required accuracy is obtained.
Jacobi method only ist iteration calculate kerni thi 2 marks
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.63
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Jacobi method 2 iteration calculation 5 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Newton Backward difference method 5 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Backward difference table construct 2 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Backward difference table construct 3 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Mcqs mix hi thy
Papers Number 02
2 Questions of 2 marks
1 question was to find the value of p from the given table using Interpolation Newton's
forward Difference Formula.
Check the Topic ”Each Topic with One Solved Question.” On page number 7
2) we have to calculate the values of x, y and z for the given system of equations using
Gauss Siede iterative method taking initial solution vector as (0, 0, 0)t?
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Both of the questions were to find the value of theta for the given matrices using Jacobi's
method ? 2x3 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
1) Construct backward difference table for the given values ….5 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
2) Find the solution of following system of Equations using Jacobi's iterative method upto
three decimal places ? 5 marks
83x+11y-4z=95
7x+52y+13z=104
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.64
3x+8y+29z=71
taking the initial solution vector as (0, 0, 0)T
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Papers Number 03
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Regula Falsi method ka aik 5 marks ka sewal
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Orthogonal metrix ki defination 3 mark ki
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Crute's method ko briefly define kerna tha 3 marks
Check the Topic ”Each Topic with One Solved Question.” On page number 7
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.65
Papers Number 04_My_Paper_Spring_2011_by_Ishfaq
Use Gause Sediel Method to find the solution of linear system of equations,perform only
two iteration correct upto 3 decimal places. 5 marks
Very easy question tha .Just remember the process.
Prove the relation u=1/4( 2 +4) 5 marks.
See the average function.
What is the formula for Newton’s backward difference interpolation...3 marks
Answered in the file
A Matrix was given, which was Symmetric .Find the value of using Jacobi’s method……
3 marks
Answered in the file
What is the formula to find value of …………. 2 marks
Answered in the file
Define Extrapolation? ……… 2 mark.
Solved_Mid-Term Past Papers
Short Question (Set-1)
my paper was power method 5 marks ,
from vactor 3 marks,
Newton_Repshon's forwrd distnce formula 5 marks.,
define extra polation 2 marks ,, and most of MCQ's from ,Newton_Repshon's method
=========================================================>
Short Question (Set-2)
1.6
2x
1
dxx 2
Approximate the integral 1 using Simpson's 1/3 rule and calculate the percentage
error. (Take result up to 4 decimal places)
Note: In order to get full marks do all necessary steps.
Construct a forward difference table for the following values
x 0.1 0.3 0.5 0.7 0.9 1.1 1.3
y 0.003 0.067 0.148 0.248 0.37 0.518 0.697
Note : In order to get full marks do all necessary steps.
Solve the system
4 x 3 y 24
3x 4 y z 30
y 4 z 24
by Gauss Seidal Method, taking (0, 0, 0)t as initial approximation(Two iterations only and
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.66
take result up to 4 decimal places )
Note : In order to get full marks do all necessary steps
2
f ( x) x
Let x , use cubic Lagrange interpolation based on the
x 0.5, x1 1, x2 2 and x3 2.5
nodes 0 to approximate f(1.5) and f(1.3).
Note : In order to get full marks do all necessary stepsSolution
Approximate the Dominant Eigenvalue and corresponding Eigenvector for the matrix
0 11 5
2 17 7
4 26 10
X 0 1,1,1
t
by using Power Method. Start with . (Five iterations only and take result up to 4
decimal places)
Note : In order to get full marks do all necessary steps
=========================================================>
Short Question (Set-4)
=========================================================>
Short Question (Set-5)
Question No: 7 ( Marks: 10 )
Interpolate the value of 0.25 using Newton’s forward difference formula.
x 0.2 0.3 0.4 0.5 0.6
F(x) 0.2304 0.2788 0.3222 0.3617 0.3979
(Perform all the necessary calculation missing calculation and steps may deduct marks.)
=========================================================>
Solved_MCQz
MCQz (Set-1)
Bisection and false position methods are also known as bracketing method and are
always
Divergent
Convergent pg no. 26
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.67
Question No. 4
The Inverse of a matrix can only be found if the matrix is
Singular
Every square non-singular matrix will have an inverse.
Scalar
Diagonal
Question No. 5
If f (x) contains trigonometric, exponential or logarithmic functions then
this equation is known as
Algebraic
Polynomial
Linear
Question No. 6
In interpolation is used to represent the δ
Forward difference ∆
Central difference
Backward difference
Question No. 7
The base of the decimal system is _______
10
0
2
8
None of the above.
Question No. 8
Question No: 2 ( Marks: 2 ) - Please choose one
Bisection method is ……………….. method
► Open Method
► Bracketing Method
=========================================================>
MCQz (Set-2)
Exact solution of 2/3 is not exists.
TRUE
FALSE
The Jacobi’s method is a method of solving a matrix equation on a matrix that has ____
zeros along its main diagonal.
no
atleast one
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.68
A 3 x 3 identity matrix have three and __________eigen values.
same
different
Eigenvalues of a symmetric matrix are all _______ .
real
complex
zero
positive
The Jacobi iteration converges, if A is strictly diagonally dominant.
TRUE
FALSE
Below are all the finite difference methods EXCEPT _________.
jacobi’s method
newton's backward difference method
Stirlling formula
Forward difference method
If n x n matrices A and B are similar, then they have the same eigenvalues (with the
same multiplicities).
TRUE
FALSE
If A is a nxn triangular matrix (upper triangular, lower triangular) or diagonal matrix ,
the eigenvalues of A are the diagonal entries of A.
TRUE
FALSE
The characteristics polynomial of a 3x 3 identity matrix is __________, if x is the eigen
values of the given 3 x 3 identity matrix. where symbol ^ shows power.
(x-1)^3
(x+1)^3
x^3-1
x^3+1
Two matrices with the same characteristic polynomial need not be similar.
TRUE
FALSE
=========================================================>
=========================================================>
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.69
MCQz (Set-6)
Question # 1 of 10 ( Start time: 04:40:08 PM ) Total Marks: 1
The determinant of a diagonal matrix is the product of the diagonal elements.
True
1. False
Question # 2 of 10 ( Start time: 04:40:58 PM ) Total Marks: 1
Power method is applicable if the eigen vectors corresponding to eigen values are
linearly independent.
True
1. false
A 3 x 3 identity matrix have three and different eigen values.
1. True
False
If n x n matrices A and B are similar, then they have the different eigenvalues (with the
same multiplicities).
1. True
False
The Jacobi’s method is a method of solving a matrix equation on a matrix that has ____
zeros along its main diagonal.
No
1. At least one
An eigenvector V is said to be normalized if the coordinate of largest magnitude is
equal to ______.
Unity
1. zero
The Gauss-Seidel method is applicable to strictly diagonally dominant or symmetric
positive definite matrices A.
True
1. False
The determinant of a _______ matrix is the product of the diagonal elements.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.70
Diagonal
1. Upper triangular
2. Lower triangular
3. Scalar
Eigenvalues of a symmetric matrix are all _________.
Real
1. Zero
2. Positive
3. Negative
The Power method can be used only to find the eigen value of A that is largest in absolute
value—we call this eigen value the dominant eigen value of A.
True
1. False
The characteristics polynomial of a 3x 3 identity matrix is __________, if x is the eigen
values of the given 3 x 3 identity matrix. where symbol ^ shows power.
(x-1)^3
1. (x+1)^3
2. x^3-1
3. x^3+1
For differences methods we require the set of values.
True
1. False
If n x n matrices A and B are similar, then they have the different eigenvalues (with
the same multiplicities).
1. True
False
If x is an eigen value corresponding to eigen value of V of a matrix A. If a is any
constant, then x – a is an eigen value corresponding to eigen vector V is an of the
matrix A - a I.
True
1. False
Central difference method seems to be giving a better approximation, however it
requires more computations.
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.71
1. True
False
Iterative algorithms can be more rapid than direct methods.
True
1. False
Central Difference method is the finite difference method.
True
1. False
=========================================================>
MCQz (Set-8)
Question # 1 of 10 ( Start time: 10:22:24 AM ) Total Marks: 1
Eigenvalues of a symmetric matrix are all _________.
Select correct option:
real
zero
positive
negative
Question # 2 of 10 ( Start time: 10:23:07 AM ) Total Marks: 1
An eigenvector V is said to be normalized if the coordinate of largest magnitude is equal
to zero.
Select correct option:
TRUE
FALSE
Question # 3 of 10 ( Start time: 10:23:55 AM ) Total Marks: 1
Exact solution of 2/3 is not exists.
Select correct option:
TRUE
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.72
FALSE
Question # 4 of 10 ( Start time: 10:24:53 AM ) Total Marks: 1
The Gauss-Seidel method is applicable to strictly diagonally dominant or symmetric
________ definite matrices A.
Select correct option:
positive
negative
Question # 5 of 10 ( Start time: 10:26:04 AM ) Total Marks: 1
Differences methods find the ________ solution of the system.
Select correct option:
numerical
Analytical
Question # 6 of 10 ( Start time: 10:26:49 AM ) Total Marks: 1
The characteristics polynomial of a 3x 3 identity matrix is __________, if x is the eigen
values of the given 3 x 3 identity matrix. where symbol ^ shows power.
Select correct option:
(x-1)^3
(x+1)^3 i m not sure about this answer
x^3-1
x^3+1
Question # 7 of 10 ( Start time: 10:28:08 AM ) Total Marks: 1
The Power method can be used only to find the eigenvalue of A that is largest in absolute
value—we call this eigenvalue the dominant eigenvalue of A.
Select correct option:
TRUE
FALSE
Question # 8 of 10 ( Start time: 10:29:33 AM ) Total Marks: 1
The Jacobi’s method is a method of solving a matrix equation on a matrix that has no
zeros along its ________.
Select correct option:
main diagonal
last column
last row
first row i m not sure about this answer
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
lOMoARcPSD|42838227
MTH603-Numerical Analysis_ Muhammad Ishfaq Page No.73
Question # 9 of 10 ( Start time: 10:30:33 AM ) Total Marks: 1
If A is a nxn triangular matrix (upper triangular, lower triangular) or diagonal matrix ,
the eigenvalues of A are the diagonal entries of A.
Select correct option:
TRUE
FALSE
Question # 10 of 10 ( Start time: 10:31:28 AM ) Total Marks: 1
A 3 x 3 identity matrix have three and different eigen values.
Select correct option:
TRUE
FALSE
=========================================================>
MCQz (Set-12)
=========================================================>
MCQz (Set-13)
=========================================================>
Share your feedback/comments at pak.nchd@gmail.com to improve file|| Back to TOP || File Version v1.5.23 published for Midterm only
Downloaded by Kamran (kamran473473@gmail.com)
You might also like
- Lecture Notes For Math-CSE 451: Introduction To Numerical Computation100% (1)Lecture Notes For Math-CSE 451: Introduction To Numerical Computation102 pages
- NAVJOT SINGH (Roll No. - 2130014010224)No ratings yetNAVJOT SINGH (Roll No. - 2130014010224)132 pages
- MTH603 Final Term (GIGA FILE by Ishfaq V11.02.02) 2 PDFNo ratings yetMTH603 Final Term (GIGA FILE by Ishfaq V11.02.02) 2 PDF171 pages
- Ba-Bsc-Bcom-4-Sem-Computer-Science-Computer-Oriented-Numerical TechNo ratings yetBa-Bsc-Bcom-4-Sem-Computer-Science-Computer-Oriented-Numerical Tech46 pages
- Numerical Analysis - Spring 2021 (Final) - 2No ratings yetNumerical Analysis - Spring 2021 (Final) - 210 pages
- Course Title: Course Code: Program: Department: College: InstitutionNo ratings yetCourse Title: Course Code: Program: Department: College: Institution7 pages
- Computer Oriented Numerical Methods (CONM) 2620004No ratings yetComputer Oriented Numerical Methods (CONM) 26200043 pages
- Handout and Home AssignmentMCA2018BATCH2NDSEMNo ratings yetHandout and Home AssignmentMCA2018BATCH2NDSEM6 pages
- EL - 504 Numerical Methods - DC-II MathsNo ratings yetEL - 504 Numerical Methods - DC-II Maths3 pages
- Faculty of Eng. & Natural Sci.: InstructorNo ratings yetFaculty of Eng. & Natural Sci.: Instructor2 pages
- School of Science & Technology Department of MathematicsNo ratings yetSchool of Science & Technology Department of Mathematics3 pages
- (WWW - Entrance-Exam - Net) - PTU MCA 3rd Semester Sample Paper 2No ratings yet(WWW - Entrance-Exam - Net) - PTU MCA 3rd Semester Sample Paper 22 pages
- Numerical Methods-EEN332 & EEN271 Spring 2015: Midterm Exam-2No ratings yetNumerical Methods-EEN332 & EEN271 Spring 2015: Midterm Exam-24 pages
- Kanojo No Omoide (Selected Scans + An Official English Translation)No ratings yetKanojo No Omoide (Selected Scans + An Official English Translation)256 pages
- A Survey of Deep Learning and Foundation Models For Time Series ForecastingNo ratings yetA Survey of Deep Learning and Foundation Models For Time Series Forecasting35 pages
- EC Maths Grade 10 November 2020 P1 and MemoNo ratings yetEC Maths Grade 10 November 2020 P1 and Memo17 pages
- Differential Pressure Transmitter: Small Flange Remote Seal TypeNo ratings yetDifferential Pressure Transmitter: Small Flange Remote Seal Type16 pages
- Midterm Interim Assessment - Tourism PlanningNo ratings yetMidterm Interim Assessment - Tourism Planning4 pages
- Exp - No:9 Design and Implementation of Lag, Lead Compensator AimNo ratings yetExp - No:9 Design and Implementation of Lag, Lead Compensator Aim10 pages
- Design of Reverse Osmosis Desalination Plant in Suez CityNo ratings yetDesign of Reverse Osmosis Desalination Plant in Suez City10 pages
- Living Wax Museum Overall Research Project Rubric Requirement PointsNo ratings yetLiving Wax Museum Overall Research Project Rubric Requirement Points11 pages
- Blog Smarter, Not Harder: SEO, Blogging, and AI Strategies to Skyrocket Your TrafficFrom EverandBlog Smarter, Not Harder: SEO, Blogging, and AI Strategies to Skyrocket Your TrafficNo ratings yet