Numerical methods. Examples with solution
|
|
- Aileen McBride
- 5 years ago
- Views:
Transcription
1 Numerical methods Examples with solution
2 CONTENTS Contents. Nonlinear Equations 3 The bisection method Newton s method Linear Systems LU-factorization Jacobi iterative method Gauss-Seidel iterative method Interpolation and approximation of functions 7 Polynom interpolation The least square method Numerical integration Trapezoid rule Composite trapezoid rule Simpson s rule Composite Simpson s rule Initial value problems for ODEs 3 Euler s method Runge-Kutta method
3 . NONLINEAR EQUATIONS. Nonlinear Equations 3
4 . NONLINEAR EQUATIONS Example.. Find the roots of equation by bisection method: x 3 = ln(0 x) with the precision ε = 0 3. Solution: x 3 ln(0 x) x f(x) It seen the root (x-intercept) lies in the interval, (from the graph,). The function values are tabulated on this interval f(x) = x 3 ln(0 x) and is found the sign of functional values are changing between. and.3. The function is continuous on a given interval. The equation has one root in the interval.,.3. 4
5 . NONLINEAR EQUATIONS The root are computed in the interval.,.3, thus a 0 =., b 0 =.3. Compute the first approximation x : x = b0 + a 0 =.3 +. =.5. Compute the functional values of the function f(x) = x 3 ln(0 x) in points a 0, x, b 0 : f(a 0 ) = 0.446, f(x ) = 0.59, f(b 0 ) = and find the interval a, b : f(a 0 )f(x ) 0 a = x 0 =.5, b = b 0 =.3. Find the approximation error b0 a 0 = 0.05 ε = 0 3. The computation continues. Compute the next approximation x = a + b = =.75 and find the interval a, b : f(a ) = 0.59, f(x ) = , f(b ) = , f(a )f(x ) 0 a = x =.75 b = b =.3. Find the approximation error b a = 0.05 ε = 0 3. The computation continues. Compute the next approximation x 3 = a + b = =.875 and find the interval a 3, b 3 : f(a ) = , f(x 3 ) = , f(b ) = , f(a ) f(x 3 ) 0 a 3 = x 3 =.875, b 3 = b =.3. 5
6 . NONLINEAR EQUATIONS Find the approximation error b a = 0.05 ε = 0 3. The computation continues. Compute the next approximation x 4 = a3 + b 3 = =.938 and find the interval a 4, b 4 : f(a 3 ) = , f(x 3 ) = 0.004, f(b 3 ) = , f(a 3 )f(x 3 ) < 0 a 4 = a 3 =.875, b 4 = x 3 =.938 Find the approximation error. b 3 a 3 = ε = 0 3. The computation continues. Compute the next approximation x 4 : x 4 = a4 + b 4 = =.906 and find the interval a 5, b 5 : f(a 4 ) = f(x 4 ) = f(b 4 ) = f(a 4 ) f(x 4 ) 0 a 5 = x 4 =.906, b 5 = b 4 =.938 Find the approximation error. b 4 a 4 = ε = 0 3. The computation continues. Compute the next approximation x 5 : x 5 = a5 + b 5 = =.9 and find the interval a 6, b 6 : f(a 5 ) = f(x 5 ) = f(b 5 ) = f(a 5 ) f(x 5 ) > 0 a 6 = x 5 =.906, b 6 = b 5 =.938 Find the approximation error b4 a 4 = ε = 0 3. The computation continues. 6
7 . NONLINEAR EQUATIONS Compute the next approximation x 4 : x 6 = a6 + b 6 = =.93 Find the approximation error b6 a 6 = < ε = 0 3. The desired precision is satisfied. All information are written to the table: i a i f(a i ) x i f(x i ) b i f(b i ) b i a i / The root of the equation is: ξ =.93 ± 0.00 Calculation in Matlab: a=., b=.3 f=inline( x^3-log(0-x) ) x=(a+b)/ while abs(b-a)/ > e-3 if f(a)*f(x)<0;b=x;else a=x;end; x=(a+b)/ end given interval definition of a function f(x) initial approximation approximations in the calculation of cycle 7
8 . NONLINEAR EQUATIONS Example.. Find the roots of equation by Newton s method: x 3 = ln(0 x) with the pecision ε = 0 6. Solution: The equation has one root in the interval.,.3 (see the Example..). Compute the the first and second derivative: f(x) = x 3 ln(0 x), f (x) = 3 x + 0 x, f (x) = 6 x + (0 x), Verify assumptions of method: f(a) f (a) = = = b a and it is seen that the condition is not satisfied. Reduce the interval in which the search root is x 3 ln(0 x) It is seen that the root lies in the interval.5,.3 (from the graph) And again, verify the method conditions, this time in the interval.5,.3 : f(a) f (a) = 0.59 = < 0.05 = b a f(b) f (b) = = < 0.05 = b a
9 . NONLINEAR EQUATIONS x i f (x i ) f (x i ) It is seen from the first and second derivative values in the interval.5,.3 : Find the initial approximation x 0 = b =.3. Compute the next approximation x : f (x) > 0 na.5,.3 f (x) > 0 na.5,.3. The Newton method conditions are satisfied. x = x 0 f(x 0) =.3 f (x 0 ) = and the approximation error x x 0 = ε = 0 6. The computation continues. Compute the next approximation x : x = x f(x ) f (x ) = = and the approximation error x x = ε = 0 6. The computation continues. Compute the next approximation x 3 : x 3 = x f(x ) f (x ) = = and the approximation error x x 0 = = < ε = 0 6. The desired precision is satisfied. All information are written to the table: i x i x i x i
10 . NONLINEAR EQUATIONS The root of the equation is: ξ = ± 0 6 Calculation in Matlab: f=inline( x^3-log(0-x) ) input the function f(x) fd=inline( 3*x^+./(0-x) ) input the first derivative f (x) xs=.3 initial approximation x=xs-f(xs)/fd(xs) the next approximation computation while abs(xs-x)>e-6 approximations in the cycle computation xs=x; x=xs-f(xs)/fd(xs) end 0
11 . LINEAR SYSTEMS. Linear Systems
12 . LINEAR SYSTEMS Example.. The matrix is given. 3 A = Compute the LU-decomposition A = LU by using Gaussian elimination method with direct choice of main element. Using the LU-decomposition,compute the solution of a system of linear equations Ax = b,where b = ( 9, 59, 8), inverse matrix A and determinant deta. Solution: Compute by using Gaussian elimination method: 3 6, The matrix U is the result of the Gaussian elimination method, the matrix L we prepare of multipliers, which are written next to matrices : L = 6 0, U = In solving Ax = b first compute y from Ly = b and then from Ux = y compute x. y = 9 9 Weget : 6y +y = 59 y = 5, 3y +y 3 = 8 x 3x +x 3 = 9 x x 3 = 5 x = 3. x 3 =
13 . LINEAR SYSTEMS Inverse matrix is calculated by solving the system Aa i = e i, i =,, 3, where e i are the columns of the unit matrix I = (e,e,e 3 ). The results a i are the columns of inverse matrix, i.e. A = (a,a,a 3 ). The solution of systems is again decomposed into two steps. y =, 0, y +y = 0,, 0 Y = 6 0, 3y +y 3 = 0, 0, 3 0 x 3x +x 3 =, 0, x x 3 = 6,, 0 A = 0. x 3 =, 0, 3 0 Determinant det A be calculated as the product of the diagonal elements of the matrix L and U: deta = detldetu = ( ) (( ) ( ) ) =. Example.. Solve the system of linear equations x +x 3x 3 = 4, x +5x +x 3 = 5, 4x x +x 3 =, using Jacobi iterative method with the pecision ǫ = 0. Solution: To ensure convergence, modify the system in form with strong diagonal dominant matrix. Replace the first and third equation 4x x +x 3 =, x +5x +x 3 = 5, x +x 3x 3 = 4. 3
14 . LINEAR SYSTEMS The modified system is rewrite to iterative form x = 4 ( + x x 3 ), x = 5 (5 x x 3 ), x 3 = 3 ( 4 x x ) The recurrent formulas for Jacobi method are obtained by iteration indices. x (k+) = ( + x(k) 4 x (k) 3 ), x (k+) = (5 x(k) 5 x (k) 3 ), x (k+) 3 = ( 4 x(k) 3 x (k) ). The computation converges for any initial approximation x (0) = (x (0), x (0), x (0) 3 ), so the following initial approximation is choosen x (0) = ( 0, 0, 0 ). Substitute x (0) into the right side of the recurrent formulas. Compute continue equally x () = ( 3,, 4 3 ). x () = ( 4, 9 5, 3 ), x (3) = (.8500,.333, ), atd. Find that the true is x (k) x (k ) R 0 3, compute thus x () x (0) R = max{ 3 0, 0, 4 0 } = 3, 3 x () x () R = max{ 4 + 3, 9, } =, x (3) x () R = , etc. The results are written to the table: Stop after the eleventh iteration, because x () x (0) R = , and the result is x = 3.00 ± 0, x =.00 ± 0, x 3 =.00 ± 0. 4
15 . LINEAR SYSTEMS k x (k) x (k) x (k) 3 x (k) x (k ) R Example.3. Solve the system of linear equations x +x 3x 3 = 4, x +5x +x 3 = 5, 4x x +x 3 =, using the Gauss-Seidel iterative methods with the pecision ǫ = 0. Solution: To ensure convergence, modify the system in form with strong diagonal dominant matrix. Replace the first and third equation 4x x +x 3 =, x +5x +x 3 = 5, x +x 3x 3 = 4. The modified system is rewrite to iterative form x = ( + x 4 x 3 ), x = (5 x 5 x 3 ), x 3 = ( 4 x 3 x ) 5
16 . LINEAR SYSTEMS The recurrent formulas for Jacobi method are obtained by iteration indices. x (k+) = ( + x(k) 4 x (k) 3 ), x (k+) = (5 x(k+) 5 x (k) 3 ), x (k+) 3 = ( 4 x(k+) 3 x (k+) ). The computation converges for any initial approximation x (0) = (x (0), x (0), x (0) 3 ), so the following initial approximation is choosen x (0) = ( 0, 0, 0 ). Substitute x (0) into the right side of the recurrent formulas. Compute x () = ( ,.000,.0667 ). continue equally x () = (.9833,.9800, ), atd. Find that the true is x (k) x (k ) R 0, compute thus x () x (0) R = max{ 3 0,. 0, } = 3, x () x () R = max{ ,.98., } = 0.00, atd. The results are written to the table: k x (k) x (k) x (k) 3 x (k) x (k ) R Stop after the four iteration, because x (4) x (3) R = , and the result is x = 3.00 ± 0, x =.00 ± 0, x 3 =.00 ± 0. 6
17 3. INTERPOLATION AND APPROXIMATION OF FUNCTIONS 3. Interpolation and approximation of functions 7
18 3. INTERPOLATION AND APPROXIMATION OF FUNCTIONS Example 3.. The nodes x i and function values f i are given below. i 0 x i 3 f i Find the interpolation polynom for theese datas: a) in basic form, b) in Lagrang form, c) in Newton form. Solution: a)interpolating polynomial is finding in form p(x) = a 0 + a x + a x. Replace the interpolation conditions p(x i ) = f i, i = 0,, p( 3) = = a 0 3a + 9a =, p( ) = = a 0 a + a =, p() = = a 0 + a + a =. These equalities represent system of linear equations, which can be written in the form 3 9 a 0 a = a. The results are a 0 =, a 8 = 3, a = 5, the interpolation polynomial has the 8 form p(x) = x x. b) Interpolation polynomial in Lagrange form is prescribed p(x) = f 0 ϕ 0 (x) + f ϕ (x) + f ϕ (x), where ϕ 0 (x), ϕ (x), ϕ (x) are Lagrange basis polynomials: ϕ 0 (x) = (x + )(x ) ( 3 + )( 3 ) = (x + )(x ), 8 8
19 3. INTERPOLATION AND APPROXIMATION OF FUNCTIONS ϕ (x) = ϕ (x) = (x + 3)(x ) ( + 3)( ) = (x + 3)(x ), 4 (x + 3)(x + ) = (x + 3)(x + ). ( + 3)( + ) 8 It is obtained p(x) = 8 (x + )(x ) + 4 (x + 3)(x ) + (x + 3)(x + ). 4 c) Interpolation polynomial in Newton form is prescribed p(x) = f 0 + f[x, x 0 ](x x 0 ) + f[x, x, x 0 ](x x 0 )(x x ), where f[x, x 0 ] a f[x, x, x 0 ] are the relative difference st and nd order. The results are written to the table: i x i f i st order nd order Using the computed values in the first row of the table, find interpolation polynomial: p(x) = (x + 3) + 5 (x + 3)(x + ). 8 9
20 3. INTERPOLATION AND APPROXIMATION OF FUNCTIONS Example 3.. Aproximate the following data x i y i by the line ϕ(x) = c x + c by the least squares method. Solution: Find the system of linear equations for the unknown coefficients c, c. 5 5 c x i + c x i = i= i= 5 5 c x i + c = i= i= 5 y i x i i= 5 y i i= Compute individual totals: 5 x i = = 90.5 i= 5 x i = = 8.5 i= 5 y i x i = = 57 i= 5 y i = = 8.5 i= Replace the system of linear equations: 90.5c + 8.5c = c + 5c = 8.5 Solve the system and the coefficients c =.3647, c = are obtained. The line ϕ(x) =.3647x is the result. 0
21 3. INTERPOLATION AND APPROXIMATION OF FUNCTIONS 6 4 zadane body nalezena funkce Calculation in Matlab: x=[,.5,3,5,7] y=[0,,6.5,6,5] n=length(x) imput x i imput y i number of points A=[sum(x.^) sum(x);sum(x) n] matrix system b=[sum(y.*x);sum(y)] c=a\b f=c()*x+c() plot(x,y, o,x,f) vector of right-hand sides calculation of the coefficients line values in points x i graph
22 4. NUMERICAL INTEGRATION 4. Numerical integration
23 4. NUMERICAL INTEGRATION Example 4.. Compute the integral by composite trapezoid tule for n = 8: 3 e x dx Solution: Compute the step h in the interval a, b, h = b a n = 3 ( ) 8 = 0.5. Write to the table division points x 0,...,x 8 and the values of integrated functions in these points. i x i y i = e x i Finally, compute the approximate value of the integral: I 0.5 ( ) ( ) n y0 + y n y0 + y 8 7 = h + y i = y i = i= i= ( = ( ) ) = = The result is: I 0.5 =
24 4. NUMERICAL INTEGRATION Calculation in Matlab: a=-,b=3 interval bounds a, b n=8 number of segments on dividing n h=(b-a)/n x=a:h:b computation of step h computation of x i y=exp(x.^) computation of y i I=h*((y()+y(n+))/+sum(y(:n))) computation of I = h ( y 0 +y n + ) n i= y i In Matlab is indexed from, and we index points from 0, so the indexes in the formula in Matlab are shifted by. 4
25 4. NUMERICAL INTEGRATION Example 4.. Compute the integral by composite trapezoid rule with the precision ε = 0 3 : 3 e x dx Solution: Compute the integral composite trapezoid formula for n =, 4, 8, 6,..., until the error is less than the precision, i.e. I h I h < ε. n h I h I h I h ε = 0 3 = < 0.00 The result is: 3 e x = ±
26 4. NUMERICAL INTEGRATION Example 4.3. Compute the integral by composite Simpson rule for n = 6: 0 arctg (x) + x + dx Solution: Compute the step h in the interval a, b, h = b a n = 0 = Write to the table division points x 0,...,x 6 and the value of integrated functions in these points. i x i y i = arctg (x i ) + x i i x i y i = arctg (x i ) + x i Finally, compute the approximate value of the integral: I 0.5 = h 3 (y 0 + y n + 4 S L + S S ) = (y 0 + y (y + y 3 + y 5 + +y 7 + y 9 + y + y 3 + y 5 ) + (y + y 4 + y 6 + y 8 + y 0 + y + y 4 ) = = ( ) =
27 4. NUMERICAL INTEGRATION The result is: I 0.5 = 5.00 Calculation in Matlab: a=0,b= interval bounds a, b n=6 number of segments on dividing n h=(b-a)/n computation of steph x=a:h:b computation of x i y=atan(x).^+x+ computation of y i SL=sum(y(::n)) computation of S L SS=sum(y(3::n-)) computation of S S I=h/3*(y()+y(n+)+4*SL+*SS) computation of I In Matlab is indexed from, and we index points from 0, so the indexes in the formula in Matlab are shifted by. 7
28 4. NUMERICAL INTEGRATION Example 4.4. For the function f(x) = sin x + arctanx + x 4 compute approximately the first derivative on the interval 0, by the formula f (x) f(x + h) f(x h) h with a step h = 0.5 a h = 0.. Compare the results with values of exact calculated derivative. Solution: The exact derivative is found by the formula: f (x) = x cosx + +x 4 +x 4x 3 arctan x ( + x 4 ). Find the nodes x i = ih, i = 0,..., 8 and compute the functional values in them f i = f(x i ) to compute approximate derivative with step h = 0.5. There are the first three columns in the table: i x i f i f[x i+,x i ] f (x i ) e i The formula for calculation of approximate derivative is written in form f (x i ) f i+ f i h = f[x i+, x i ], 8
29 4. NUMERICAL INTEGRATION. Compute the approximate derivative values in all internal nodes, eg f (0.75) =.008;. This is the fourth column of the table. The fifth column contains the exact values of derivatives that are in the last column compared with the approximate values: e i = f[x i+, x i ] f (x i ). For h = 0. the procedure is analogous, we obtain more accurate values, as shown in the following images: h= derivace aproximace derivace h= derivace aproximace derivace
30 4. NUMERICAL INTEGRATION Calculation in Matlab: a=0,b= interval bounds a, b h=0.5 step h x=[a:h:b] computation of nodes x i n=length(x) number of nodes n + df= *x.*cos(x.^) +((+x.^4)./(+x.^)-... 4*x.^3.*atan(x))./(+x.^4).^ computation of f (x i ) f=sin(x.^)+atan(x)./(+x.^4) computation of f i = f(x i ) adf=[nan;(f(3:n)-f(:n-))/(*h);nan] approximate derivative tab=[[0:n-],x,f,adf,df,abs(adf-df)] table plot(x,df, :,x,adf, - ) graph 30
31 5. INITIAL VALUE PROBLEMS FOR ODES 5. Initial value problems for ODEs 3
32 5. INITIAL VALUE PROBLEMS FOR ODES Example 5.. Solve the differential equation with the initial condition: y = 0. y, y( ) =, x + by the Euler method in the interval, 3 with the step h = 0.5. Solution: First, compute the number of dividing segments in the interval a, b : n = b a h = 3 ( ) 0.5 = 0. The initial condition is x 0 =,, y 0 =. Compute x, y by the Euler metod: x = x 0 + h = =.5 ( ) ( ) y = y 0 + h x y 0 = ( ) = 0.85 ( ) + Compute x, y, x 3, y 3 : x y = x + h = = ( ) ( ) = y + h x + 0. y = ( 0.85) = (.5) + = x 3 = x + h = = 0.5 ( ) ( ) y 3 = y + h x + 0. y = ( ) = ( ) + = 0.37 Other approximate values of solution x 4, y 4,..., x 0, y 0 compute accordingly. 3
33 5. INITIAL VALUE PROBLEMS FOR ODES The approximate values of the solution are in the table and graph: i x i y i y Priblizne reseni x Calculation in Matlab: a=-,b=3 h=0.5 n=(b-a)/h x()=-,y()=- for i=:n x(i+)=x(i)+h y(i+)=y(i)+h*(/(x(i)^+)-0.*y(i)) end [x y ] plot(x,y) interval bounds a, b input step h number of dividing segments n input initial condition computation of x i, y i in the cycle extract into the table graph 33
34 5. INITIAL VALUE PROBLEMS FOR ODES Example 5.. Solve the differential equation with the initial condition: y = sin(x) y x 3 + y, y() = 0, by the 4th order Runge Kutta method in the interval, with the step h = 0.. Solution: First, compute the number of dividing segments in the interval a, b : n = b a h The initial condition is x 0 =,, y 0 = 0. = 0. = 5. Compute x, y by the 4th order Runge Kutta method: k = h f (x 0, y 0 ) = h (sin (x 0 ) y 0 x y 0) = 0. ( sin () ) = k k 3 = 0. ( = h f x 0 + h, y 0 + k ) = ( = h sin x 0 + h ) ( y 0 + k ) ( x 0 + h ) 3 + y 0 + k = ( ( = 0. sin + 0. ) ( ) ( + 0. ) ) = ( = h f x 0 + h, y 0 + k ) = ( = h sin x 0 + h ) ( y 0 + k ) ( x 0 + h ) 3 + y 0 + k = ( ( = 0. sin + 0. ) ( ) ( + 0. ) ) = = k 4 = h f (x 0 + h, y 0 + k 3 ) = = h (sin ) (x 0 + h) (y 0 + k 3 ) (x 0 + h) 3 + y 0 + k 3 = = 0. (sin ( + 0.) ( ) ( + 0.) ) = x = x 0 + h = + 0. =. 34
35 5. INITIAL VALUE PROBLEMS FOR ODES y = y (k + k + k 3 + k 4 ) = = 0 + ( 0. + ( 0.304) + ( 0.337) ) = Compute x, y : k = h f (x, y ) = h (sin (x ) y x 3 + y ) = = 0. (sin (.) ( 0.3) ) = = ( k = h f x + h, y + k ) = ( = h sin x + h ) ( y + k ) ( x + h ) 3 + y + k = ( ( = 0. sin ) ( ) ( ) ) = ( k 3 = h f x + h, y + k ) = ( = h sin x + h ) ( y + k ) ( x + h ) 3 + y + k = ( ( = 0. sin ) ( ) ( + 0. ) ) = k 4 = h f (x + h, y + k 3 ) = h (sin ) (x + h) (y + k 3 ) (x + h) 3 + y + k 3 = = 0. (sin ( + 0.) ( ) ( + 0.) ) = 0.95 x = x + h = =.4 y = y + 6 (k + k + k 3 + k 4 ) = = ( ( ) + ( ) 0.95) = Other approximate values of solution x 3, y 3,..., x 5, y 5 compute accordingly. 35
36 5. INITIAL VALUE PROBLEMS FOR ODES The approximate values of the solution are in the table and graph: 0 Priblizne reseni i x i y i y x Calculation in Matlab: f=inline( sin(x)*y-x^3+y ) right sidef(x, y) a=,b= h=0. n=(b-a)/h x()=,y()=0 for i=:n k=h*f(x(i),y(i)) k=h*f(x(i)+h/,y(i)+k/) k3=h*f(x(i)+h/,y(i)+k/) k4=h*f(x(i)+h,y(i)+k3) x(i+)=x(i)+h y(i+)=y(i)+/6*(k+*k+*k3+k4) end [x y ] plot(x,y) interval bounds a, b input step h number of dividing segments n input initial condition computation of x i, y i in the cycle extract into the table graph 36
Hence a root lies between 1 and 2. Since f a is negative and f(x 0 ) is positive The root lies between a and x 0 i.e. 1 and 1.
The Bisection method or BOLZANO s method or Interval halving method: Find the positive root of x 3 x = 1 correct to four decimal places by bisection method Let f x = x 3 x 1 Here f 0 = 1 = ve, f 1 = ve,
More informationM.SC. PHYSICS - II YEAR
MANONMANIAM SUNDARANAR UNIVERSITY DIRECTORATE OF DISTANCE & CONTINUING EDUCATION TIRUNELVELI 627012, TAMIL NADU M.SC. PHYSICS - II YEAR DKP26 - NUMERICAL METHODS (From the academic year 2016-17) Most Student
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationPrevious Year Questions & Detailed Solutions
Previous Year Questions & Detailed Solutions. The rate of convergence in the Gauss-Seidal method is as fast as in Gauss Jacobi smethod ) thrice ) half-times ) twice 4) three by two times. In application
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More informationToday s class. Linear Algebraic Equations LU Decomposition. Numerical Methods, Fall 2011 Lecture 8. Prof. Jinbo Bi CSE, UConn
Today s class Linear Algebraic Equations LU Decomposition 1 Linear Algebraic Equations Gaussian Elimination works well for solving linear systems of the form: AX = B What if you have to solve the linear
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationMTH603 FAQ + Short Questions Answers.
Absolute Error : Accuracy : 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
More informationVirtual University of Pakistan
Virtual University of Pakistan File Version v.0.0 Prepared For: Final Term Note: Use Table Of Content to view the Topics, In PDF(Portable Document Format) format, you can check Bookmarks menu Disclaimer:
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More information5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns
5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns (1) possesses the solution and provided that.. The numerators and denominators are recognized
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationMathematical Methods for Numerical Analysis and Optimization
Biyani's Think Tank Concept based notes Mathematical Methods for Numerical Analysis and Optimization (MCA) Varsha Gupta Poonam Fatehpuria M.Sc. (Maths) Lecturer Deptt. of Information Technology Biyani
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationTMA4125 Matematikk 4N Spring 2017
Norwegian University of Science and Technology Institutt for matematiske fag TMA15 Matematikk N Spring 17 Solutions to exercise set 1 1 We begin by writing the system as the augmented matrix.139.38.3 6.
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationExamination paper for TMA4125 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA45 Matematikk 4N Academic contact during examination: Anne Kværnø a, Louis-Philippe Thibault b Phone: a 9 66 38 4, b 9 3 0 95 Examination date:
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationMath 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018
1 Linear Systems Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March, 018 Consider the system 4x y + z = 7 4x 8y + z = 1 x + y + 5z = 15. We then obtain x = 1 4 (7 + y z)
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationNumerical and Statistical Methods
F.Y. B.Sc.(IT) : Sem. II Numerical and Statistical Methods Time : ½ Hrs.] Prelim Question Paper Solution [Marks : 75 Q. Attempt any THREE of the following : [5] Q.(a) What is a mathematical model? With
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationNumerical Marine Hydrodynamics Summary
Numerical Marine Hydrodynamics Summary Fundamentals of Digital Computing Error Analysis Roots of Non-linear Equations Systems of Linear Equations Gaussian Elimination Iterative Methods Optimization, Curve
More informationNumerical and Statistical Methods
F.Y. B.Sc.(IT) : Sem. II Numerical and Statistical Methods Time : ½ Hrs.] Prelim Question Paper Solution [Marks : 75 Q. Attempt any THREE of the following : [5] Q.(a) What is a mathematical model? With
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationINSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad
1 P a g e INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 04 Name : Mathematics-II Code : A0006 Class : II B. Tech I Semester Branch : CIVIL Year : 016 017 FRESHMAN ENGINEERING
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad MECHANICAL ENGINEERING TUTORIAL QUESTION BANK
Course Name Course Code Class Branch INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 500 043 Mathematics-II A30006 II-I B. Tech Freshman Engineering Year 016 017 Course Faculty MECHANICAL ENGINEERING
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More informationNUMERICAL METHODS USING MATLAB
NUMERICAL METHODS USING MATLAB Dr John Penny George Lindfield Department of Mechanical Engineering, Aston University ELLIS HORWOOD NEW YORK LONDON TORONTO SYDNEY TOKYO SINGAPORE Preface 1 An introduction
More informationExamination paper for TMA4130 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA4130 Matematikk 4N Academic contact during examination: Morten Nome Phone: 90 84 97 83 Examination date: 13 December 2017 Examination time (from
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationFundamental Numerical Methods for Electrical Engineering
Stanislaw Rosloniec Fundamental Numerical Methods for Electrical Engineering 4y Springei Contents Introduction xi 1 Methods for Numerical Solution of Linear Equations 1 1.1 Direct Methods 5 1.1.1 The Gauss
More informationTutorial 1, B. Tech. Sem III, 24 July, (Root Findings and Linear System of Equations)
Tutorial, B. Tech. Sem III, 24 July, 206 (Root Findings and Linear System of Equations). Find the root of the equation e x = 3x lying in [0,] correct to three decimal places using Bisection, Regula-Falsi
More informationA = (a + 1) 2 = a 2 + 2a + 1
A = (a + 1) 2 = a 2 + 2a + 1 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 A = ( (a + b) + 1 ) 2 = (a + b) 2 + 2(a + b) + 1 = a 2 + 2ab + b 2 + 2a + 2b + 1 3 A = (
More informationSolutions to Exam 2, Math 10560
Solutions to Exam, Math 6. Which of the following expressions gives the partial fraction decomposition of the function x + x + f(x = (x (x (x +? Solution: Notice that (x is not an irreducile factor. If
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationSchool of Sciences Indira Gandhi National Open University Maidan Garhi, New Delhi (For January 2012 cycle)
MTE-0 ASSIGNMENT BOOKLET Bachelor's Degree Programme Numerical Analysis (MTE-0) (Valid from st January, 0 to st December, 0) School of Sciences Indira Gandhi National Open University Maidan Garhi, New
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationNUMERICAL ANALYSIS I. MARTIN LOTZ School of Mathematics The University of Manchester. May 2016
NUMERICAL ANALYSIS I by MARTIN LOTZ School of Mathematics The University of Manchester May 06 Contents Contents ii Week. Computational Complexity....................... Accuracy...............................
More informationSTATISTICS AND NUMERICAL METHODS
STATISTICS AND NUMERICAL METHODS Question IV November / December 2011 Part-A 1. The heights of college students in Chennai are normally distributed with standard deviation 6 cm and sample of 100 students
More informationChapter 5: Numerical Integration and Differentiation
Chapter 5: Numerical Integration and Differentiation PART I: Numerical Integration Newton-Cotes Integration Formulas The idea of Newton-Cotes formulas is to replace a complicated function or tabulated
More informationUSHA RAMA COLLEGE OF ENGINEERING & TECHNOLOGY
Code No: R007/R0 Set No. I B.Tech I Semester Supplementary Examinations, Feb/Mar 04 MATHEMATICAL METHODS ( Common to Civil Engineering, Electrical & Electronics Engineering, Computer Science & Engineering,
More informationNumerical Methods with MATLAB
Numerical Methods with MATLAB A Resource for Scientists and Engineers G. J. BÖRSE Lehigh University PWS Publishing Company I(T)P AN!NTERNATIONAL THOMSON PUBLISHING COMPANY Boston Albany Bonn Cincinnati
More informationPage No.1. MTH603-Numerical Analysis_ Muhammad Ishfaq
Page No.1 File Version v1.5.3 Update: (Dated: 3-May-011) This version of file contains: Content of the Course (Done) FAQ updated version.(these must be read once because some very basic definition and
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationUnit I (Testing of Hypothesis)
SUBJECT NAME : Statistics and Numerical Methods SUBJECT CODE : MA645 MATERIAL NAME : Part A questions REGULATION : R03 UPDATED ON : November 07 (Upto N/D 07 Q.P) Unit I (Testing of Hypothesis). State level
More informationNumerical Analysis: Solutions of System of. Linear Equation. Natasha S. Sharma, PhD
Mathematical Question we are interested in answering numerically How to solve the following linear system for x Ax = b? where A is an n n invertible matrix and b is vector of length n. Notation: x denote
More informationApplied Numerical Analysis (AE2220-I) R. Klees and R.P. Dwight
Applied Numerical Analysis (AE0-I) R. Klees and R.P. Dwight February 018 Contents 1 Preliminaries: Motivation, Computer arithmetic, Taylor series 1 1.1 Numerical Analysis Motivation..........................
More informationSuggested solutions, TMA4125 Calculus 4N
Suggested solutions, TMA5 Calculus N Charles Curry May 9th 07. The graph of g(x) is displayed below. We have b n = = = 0 [ nπ = nπ ( x) nπx dx nπx dx cos nπx ] x nπx dx [ nπx x cos nπ ] ( cos nπ + cos
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationCHAPTER 8: Matrices and Determinants
(Exercises for Chapter 8: Matrices and Determinants) E.8.1 CHAPTER 8: Matrices and Determinants (A) means refer to Part A, (B) means refer to Part B, etc. Most of these exercises can be done without a
More informationExamination paper for TMA4122/TMA4123/TMA4125/TMA4130 Matematikk 4M/N
Department of Mathematical Sciences Examination paper for TMA4122/TMA4123/TMA4125/TMA4130 Matematikk 4M/N Academic contact during examination: Markus Grasmair Phone: 97 58 04 35 Code C): Basic calculator.
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationPractice Problems for the Final Exam
Math 114 Spring 2017 Practice Problems for the Final Exam 1. The planes 3x + 2y + z = 6 and x + y = 2 intersect in a line l. Find the distance from the origin to l. (Answer: 24 3 ) 2. Find the area of
More information1 Lecture 8: Interpolating polynomials.
1 Lecture 8: Interpolating polynomials. 1.1 Horner s method Before turning to the main idea of this part of the course, we consider how to evaluate a polynomial. Recall that a polynomial is an expression
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationASSIGNMENT BOOKLET. Numerical Analysis (MTE-10) (Valid from 1 st July, 2011 to 31 st March, 2012)
ASSIGNMENT BOOKLET MTE-0 Numerical Analysis (MTE-0) (Valid from st July, 0 to st March, 0) It is compulsory to submit the assignment before filling in the exam form. School of Sciences Indira Gandhi National
More informationRoot Finding: Close Methods. Bisection and False Position Dr. Marco A. Arocha Aug, 2014
Root Finding: Close Methods Bisection and False Position Dr. Marco A. Arocha Aug, 2014 1 Roots Given function f(x), we seek x values for which f(x)=0 Solution x is the root of the equation or zero of the
More informationSOLUTION OF EQUATION AND EIGENVALUE PROBLEMS PART A( 2 MARKS)
CHENDU COLLEGE OF ENGINEERING AND TECHNOLOGY (Approved by AICTE New Delhi, Affiliated to Anna University Chennai. Zamin Endathur Village, Madurntakam Taluk, Kancheepuram Dist.-603311.) MA6459 - NUMERICAL
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 17 Numerical Integration Formulas PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationSubbalakshmi Lakshmipathy College of Science. Department of Mathematics
ॐ Subbalakshmi Lakshmipathy College of Science Department of Mathematics As are the crests on the hoods of peacocks, As are the gems on the heads of cobras, So is Mathematics, at the top of all Sciences.
More informationLösning: Tenta Numerical Analysis för D, L. FMN011,
Lösning: Tenta Numerical Analysis för D, L. FMN011, 090527 This exam starts at 8:00 and ends at 12:00. To get a passing grade for the course you need 35 points in this exam and an accumulated total (this
More informationBasics of Calculus and Algebra
Monika Department of Economics ISCTE-IUL September 2012 Basics of linear algebra Real valued Functions Differential Calculus Integral Calculus Optimization Introduction I A matrix is a rectangular array
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More informationNUMERICAL ANALYSIS. K.D. Anderson J.S.C. Prentice C.M. Villet
NUMERICAL ANALYSIS KD Anderson JSC Prentice CM Villet This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 4 International License To view a copy of this license, visit
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationReview of matrices. Let m, n IN. A rectangle of numbers written like A =
Review of matrices Let m, n IN. A rectangle of numbers written like a 11 a 12... a 1n a 21 a 22... a 2n A =...... a m1 a m2... a mn where each a ij IR is called a matrix with m rows and n columns or an
More informationWeek 1: need to know. November 14, / 20
Week 1: need to know How to find domains and ranges, operations on functions (addition, subtraction, multiplication, division, composition), behaviors of functions (even/odd/ increasing/decreasing), library
More informationLecture Notes on Introduction to Numerical Computation 1
Lecture Notes on Introduction to Numerical Computation 1 Wen Shen Fall 2012 1 These notes are provided to students as a supplement to the textbook. They contain mainly examples that we cover in class.
More informationALGEBRA II WITH TRIGONOMETRY EXAM
First Round : March 7, 00 Second Round: April, 00 at The University of Alabama ALGEBRA II WITH TRIGONOMETRY EXAM Construction of this test directed by Congxiao Liu, Alabama A&M University and Alzaki Fadl
More informationECONOMICS 207 SPRING 2006 LABORATORY EXERCISE 5 KEY. 8 = 10(5x 2) = 9(3x + 8), x 50x 20 = 27x x = 92 x = 4. 8x 2 22x + 15 = 0 (2x 3)(4x 5) = 0
ECONOMICS 07 SPRING 006 LABORATORY EXERCISE 5 KEY Problem. Solve the following equations for x. a 5x 3x + 8 = 9 0 5x 3x + 8 9 8 = 0(5x ) = 9(3x + 8), x 0 3 50x 0 = 7x + 7 3x = 9 x = 4 b 8x x + 5 = 0 8x
More information. (a) Express [ ] as a non-trivial linear combination of u = [ ], v = [ ] and w =[ ], if possible. Otherwise, give your comments. (b) Express +8x+9x a
TE Linear Algebra and Numerical Methods Tutorial Set : Two Hours. (a) Show that the product AA T is a symmetric matrix. (b) Show that any square matrix A can be written as the sum of a symmetric matrix
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationNote: Final Exam is at 10:45 on Tuesday, 5/3/11 (This is the Final Exam time reserved for our labs). From Practice Test I
MA Practice Final Answers in Red 4/8/ and 4/9/ Name Note: Final Exam is at :45 on Tuesday, 5// (This is the Final Exam time reserved for our labs). From Practice Test I Consider the integral 5 x dx. Sketch
More informationMATH 163 HOMEWORK Week 13, due Monday April 26 TOPICS. c n (x a) n then c n = f(n) (a) n!
MATH 63 HOMEWORK Week 3, due Monday April 6 TOPICS 4. Taylor series Reading:.0, pages 770-77 Taylor series. If a function f(x) has a power series representation f(x) = c n (x a) n then c n = f(n) (a) ()
More information