2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l.
|
|
- Magnus Goodman
- 5 years ago
- Views:
Transcription
1 Numerical Methods II UNIT.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS.1.1 Runge-Kutta Method of Fourth Order 1. Let = f x,y,z, = gx,y,z be the simultaneous first order ordinary differential equations with the initial values yx = y, zx = z. The modified form of expression for computation are as follows. k 1 = hfx,y,z ; l 1 = hgx,y,z k = hf x h,y k 1,z l 1 ; l = hg x h,y k 1,z l 1 k 3 = hf x h,y k,z l ; l 3 = hg x h,y k,z l k 4 = hf x h,y k 3,z l 3 ; l4 = hg x h,y k 3,z l 3 Then according to second order Runge-Kutta method yx h = y k zx h = z l According to fourth order Runge-Kutta method yx h = k1 k 6 k 3 k 4 zx h = l1 l 6 l 3 l 4
2 46 Engineering Mathematics-IV =.1[ =.11 l 4 = hg t h, x k 3, y l 3 =.1[ =.15 xt h = x 1 [ k1 k 6 k 3 k 4 yt h = y 1 6 = 1 1 [ = [ l1 l l 3 l 4 Thus, = 1 1 [ = x.1 = y.1 = Example : equations Use fourth order Runge-Kutta method to solve the system of = 1 xz; = xy given y =, z = 1 at x =.3. Solution: Given simultaneous equations are = 1 xz, = xy y =, z = 1 or, f x,y,z = 1 xz, gx,y,z = xy Here x =, y =, z = 1 given h =.3.
3 Numerical Methods II 47 Now, we shall compute the following. k 1 = hfx, y, z =.3[1 1 =.3 l 1 = hgx, y, z =.3[ = k = hf x h, y k 1, z l 1 =.3[ =.345 l = hg x h, y k 1, z l 1 =.3[ =.675 k 3 = hf x h, y k, z l =.3[ = l 3 = hf yx h = y 1 6 x h, y k, z l =.3[ =.7765 k 4 = hf x h, y k 3, z l 3 =.3[ =.3893 l 4 = hf x h, y k 3, z l 3 =.3[ =.314 [ k1 k k 3 k 4 y.3 = 1 [ =.34483
4 48 Engineering Mathematics-IV zx h = z 1 6 [ l1 l l 3 l 4 z.3 = Picard s Method = 1 1 [ Let the simultaneous differential equations be = f x,y,z = gx,y,z with initial conditions yx = y zx = z. Then = y y = y y 3 = y f x,y,z, = z f x,,, z = z f x,y,z, z 3 = z gx,y,z gx,, gx,y,z so on. Continuing this process, a sequence of functions of x, i.e.,,y,y 3...,z,z 3... are obtained each giving a better approximations of the desired solution. Example 1: Using Picard s method find approximate values of y z corresponding to x =.1 given that y =, z = 1 = x z, = x y. Solution: Given equations are = f x,y,z = x z = gx,y,z = x y y = y f x,y,z z = z gx,y,z x x
5 Numerical Methods II 49 First approximation = y = z x f x,y,z = Second approximation y = y x gx,y,z = 1 x f x,, = x 1 = x x x 4 = 1 4x 1 x x 1 4x 1 x = x 3 x x3 6 z = z Third approximation so on. x gx,, = 1 = 1 4x 3 x x 3 x4 4 x5 y 3 = y x f x,y,z [ x = x 3 x 1 x3 1 4 x4 1 x5 z 3 = z gx,y,z x 1 x 1 1 x6 = 1 4x x 3 x3 1 x4 6 x5 1 x6 5 x7
6 5 Engineering Mathematics-IV when x =.1 =.15, y =.8517 y 3 =.8447 =.65, z = z 3 =.5867 Thus, y.1 =.845, z.1 =.5867 correct to four decimal places..1.3 Taylor s Series Method Let the simultaneous differential equations be with initial conditions yx = y zx = z. If h be the step-size, = yx h = zx h Then Taylor s algorithm for 1 gives = f x,y, = gx,y,z = y hy h! y h3 = z hz h! z h3 3! y 3! z Differentiating 1 successively we get y,z, etc So the values y,y,y... z,z,z... are known. Substituting these in 3 4, we obtain, for the next step. Similarly, y = hy 1 h! y 1 h3 z = hz 1 h! z 1 h3 3! y 1 3! z Since are known, we calculate y, 1 y... 1 z, 1 z... substituting 1 these in 5 6 we get y z. Proceeding further, we can calculate the other values of y z step by step.
7 Numerical Methods II 51 Example 1: Given = z = xz y with y = 1, z =, obtain y z for x =.1,.,.3 by Taylor s series method. Hint: We have y = z z = xz y We use Taylor s series method to find y z.. NUMERICAL SOLUTION OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS..1 Picard s Method Let the second order differential equation be d y = f x,y, 1 with yx = y y x = y. Now, let = z then d y =. Substituting in 1 = f x,y,z with yx = y zx = z The problem reduces to solving the simultaneous equations = z = f 1 x,y,z = f x,y,z subject to yx = y zx = z. By Picard s method we have = y f 1 x,y,z, = z x y = y f 1 x,,, z = z x x x f x,y,z f x,,
8 Numerical Methods II 53 =.5.1. x.5x, z =.1 x.1. x.5x.5.1x =.5.1x. x3 6.5x, =.1.x.x4 4 1x3 3.5x.1 x so on. when x =.1 =.51 y =.574 = y 1 =.49 y = z =.488 y.1 =.574 y.1 =.488 Example : solution of Using Picard s method, obtain the second approximation to the Solution: d y = x3 x3 y so that y = 1, y 1 =. Given the second order differential equation is d y = x3 x3 y with y = 1, y 1 = Let = z so that d y = = z = x3 z x 3 y with y = 1, =. By Picard s methods, = f x,y,z = f 1 x,y,z
9 54 Engineering Mathematics-IV First approximation = y = 1 x f x y z = z = 1 x Second approximation y = y = 1 1 = 1 x = 1 3x4 8 x f x,, z = z 1 3x4 8 x = 1 x = 1 x 3x5 x = 1 f x, y z 1 x3 1x = x3 x f x,, x 3 1 3x4 x 3 1 x 8 x x7 x 3 x4 = 1 x4 8 3x8 64 x4 4 x5 1 Example 3: 1 = 3x4 8 x5 1 3x8 64. Using Picard s method find the solution of the equation y xy y =, with y = 1 y = for x =.1.3. Solution: Given equation is y xy y = with y = 1, y =. Let = y = z so that d y = y = = z = xz y with y = 1, z = = f x,y,z, = f 1 x,y,z.
10 56 Engineering Mathematics-IV when x =.1.1 = 1 y.1 =.955, y 3.1 = =.1 z.1 =.99, z 3.1 =.995 y.1 =.9951 z.1 = y.1 =.995 when x =.3.3 = 1 y.3 =.955 y 3.3 = =.3 z.3 =.91, z 3.3 =.911 y.3 =.9556 z.3 = y.3 = Runge-Kutta Method Let the second order differential equation be d y = f x, y,...1 with yx = y, y x = y. Now, let = z so that d y = = z = f x,y,z = f 1 x,y,z, = f x,y,z with yx = y, zx = z. The problem reduces to solving the simultaneous equations = z = f 1 x,y,z = f x,y,z with yx = y zx = z. Starting at x,y,z taking the step-sizes for x, y, z to be h, k, l respectively, the Runge-Kutta method gives
11 Numerical Methods II 57 K 1 = hf 1 x,y,z l 1 = hf x,y,z K = hf 1 x 1 h, y k 1, z l 1 k 3 = hf 1 x h, y k, z l k 4 = hf 1 x h, y k 3, z l 3 l = hf x h, y k 1, z l 1 l 3 = hf x h, y k, z l l 4 = hf x h, y h 3, z l 3 = y 1 6 k1 k k 3 k 4 = z h 6 l1 l l 3 l 4 Example 1: Using Runge-Kutta method solve y = xy y for x =. correct to 4 decimal places initial conditions are x =, y = 1, y =. Solution: Let Given equation of second order is = z = f 1 x,y,z, then y = xy y with x =, y = 1, y = = xz y = f x,y,z with x =, y = 1, z = we compute k 1, k, k 3, k 4 for f 1 x, y, z l 1, l, l 3, l 4 for f x, y, z. by Runge-Kutta formulae, we have k 1 = hf 1 x,y,z l 1 = hf x,y,z =. = =. 1 =. k = hf 1 x h, y k 1, z l 1 l = hf x h, y k 1, z l 1 =..1 =. =..999 =.1998 k 3 = hf 1 x h, y k, z l l 3 = hf x h, y k, z l =..999 =. = =.1958
12 Numerical Methods II 59 Let = z then d y =. Substituting in 1 = f x,y,z with yx = y zx = z The problem reduces to solving the simultaneous equations. = z = f 1 x,y,z = f x,y,z subject to yx = y, zx = z. The predictor corrector formulae can be written as y P n 1 = y 4h [ n 3 y 3 n y n 1 y n z P n 1 = z 4h [ n 3 z 3 n z n 1 z n y c n 1 = y h [ n 1 y 3 n 1 4y n y n 1 z c n 1 = z h [ n 1 z 3 n 1 4z n z n 1 In particular when n = 3 y P 4 = y 4h 3 z P 4 = z 4h 3 [ y 1 y y 3 [ z 1 z z 3 y c 4 = y h [ y 3 4y 3 y 4 z c 4 = z h [ z 3 4z 3 z 4 3 is called Milne s predictor formula 3 is called Milne s corrector formula.
13 6 Engineering Mathematics-IV Example 1: Given d y x y =, y = 1, y = obtain y for x=,.1,.,.3 by any method. Further compute y.4 by Milne s method. Solution: Let y = z then y = z The given equation is y xy y =...1 with y = 1, y = 1 z xz y = y = z = y = z = f 1 x,y,z with y = 1 = z = f x,y,z = xz y with z = We use Taylor s series method to find y. Differentiating 1 n times we get At x =, y n = n 1y n y n xy n 1 ny n y n = y = 1 gives y = 1, y 4 = gives y 3 = y 5 =... =. Exp yx by Taylor s series = 3, y 6 = 5 3 yx = y x x! y x3 3! y 3... yx = 1 x 3 5 3! 4! x4 x ! From we have zx = y x = x 1 x3 1 8 x5... = xy 3 y.1 = =.995 y. = 1. y.3 = = =.956.
14 6 Engineering Mathematics-IV y.4 = y. h 3 =.98 =.93 [ y. 4y.3 y.4 Hence, y.4 =.93 y.4 = [ Example : Given d y x y = y =.5,y =.1 obtain y for x =,.1,.,.3 by any suitable method. Hence, apply Milne s method to compute y.4. Example 3: Given y = xy y y = 1, y = obtain y for x =,.1,.,.3 by Taylor s series method. Hence, apply Milne s method to compute y.4. Ans: y.4= Taylor s Series Method Example 1: Given y xy y =, y = 1, y =. Obtain y for x =,.1,.,.3 by Taylor s series method. Solution: Let y = z then y = z Given second order differential equation is We use Taylor s series method to find y. y xy y = 1 z xz y = ; y = z Differentiating the given equation n times we get At x =, y n xy n 1 ny n y n = y n = n 1y n y = 1 gives y = 1, y 4 = 3, y 6 = 5 3 gives y 3 = y 5 =... =
Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles.
» ~ $ ) 7 x X ) / ( 8 2 X 39 ««x» ««! «! / x? \» «({? «» q «(? (?? x! «? 8? ( z x x q? ) «q q q ) x z x 69 7( X X ( 3»«! ( ~«x ««x ) (» «8 4 X «4 «4 «8 X «x «(» X) ()»» «X «97 X X X 4 ( 86) x) ( ) z z
More informationMANY BILLS OF CONCERN TO PUBLIC
- 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -
More informationLOWELL WEEKI.Y JOURINAL
/ $ 8) 2 {!»!» X ( (!!!?! () ~ x 8» x /»!! $?» 8! ) ( ) 8 X x /! / x 9 ( 2 2! z»!!»! ) / x»! ( (»»!» [ ~!! 8 X / Q X x» ( (!»! Q ) X x X!! (? ( ()» 9 X»/ Q ( (X )!» / )! X» x / 6!»! }? ( q ( ) / X! 8 x»
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
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 informationLOWELL WEEKLY JOURNAL
G $ G 2 G ««2 ««q ) q «\ { q «««/ 6 «««««q «] «q 6 ««Z q «««Q \ Q «q «X ««G X G ««? G Q / Q Q X ««/«X X «««Q X\ «q «X \ / X G XX «««X «x «X «x X G X 29 2 ««Q G G «) 22 G XXX GG G G G G G X «x G Q «) «G
More informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
More informationLOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort
- 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [
More information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More information1 Partial differentiation and the chain rule
1 Partial differentiation and the chain rule In this section we review and discuss certain notations and relations involving partial derivatives. The more general case can be illustrated by considering
More informationLOWELL. MICHIGAN. WEDNESDAY, FEBRUARY NUMllEE 33, Chicago. >::»«ad 0:30am, " 16.n«l w 00 ptn Jaekten,.'''4snd4:4(>a tii, ijilwopa
4/X6 X 896 & # 98 # 4 $2 q $ 8 8 $ 8 6 8 2 8 8 2 2 4 2 4 X q q!< Q 48 8 8 X 4 # 8 & q 4 ) / X & & & Q!! & & )! 2 ) & / / ;) Q & & 8 )
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More information2x (x 2 + y 2 + 1) 2 2y. (x 2 + y 2 + 1) 4. 4xy. (1, 1)(x 1) + (1, 1)(y + 1) (1, 1)(x 1)(y + 1) 81 x y y + 7.
Homework 8 Solutions, November 007. (1 We calculate some derivatives: f x = f y = x (x + y + 1 y (x + y + 1 x = (x + y + 1 4x (x + y + 1 4 y = (x + y + 1 4y (x + y + 1 4 x y = 4xy (x + y + 1 4 Substituting
More informationPanHomc'r I'rui;* :".>r '.a'' W"»' I'fltolt. 'j'l :. r... Jnfii<on. Kslaiaaac. <.T i.. %.. 1 >
5 28 (x / &» )»(»»» Q ( 3 Q» (» ( (3 5» ( q 2 5 q 2 5 5 8) 5 2 2 ) ~ ( / x {» /»»»»» (»»» ( 3 ) / & Q ) X ] Q & X X X x» 8 ( &» 2 & % X ) 8 x & X ( #»»q 3 ( ) & X 3 / Q X»»» %» ( z 22 (»» 2» }» / & 2 X
More informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationA. H. Hall, 33, 35 &37, Lendoi
7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationInitial Value Problems
Numerical Analysis, lecture 13: Initial Value Problems (textbook sections 10.1-4, 10.7) differential equations standard form existence & uniqueness y 0 y 2 solution methods x 0 x 1 h h x 2 y1 Euler, Heun,
More informationLOWELL WEEKLY JOURNAL
Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q
More informationd A L. T O S O U LOWELL, MICHIGAN. THURSDAY, DECEMBER 5, 1929 Cadillac, Nov. 20. Indignation
) - 5 929 XXX - $ 83 25 5 25 $ ( 2 2 z 52 $9285)9 7 - - 2 72 - - 2 3 zz - 9 86 - - - - 88 - q 2 882 q 88 - - - - - - ( 89 < - Q - 857-888 - - - & - - q - { q 7 - - - - q - - - - - - q - - - - 929 93 q
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 informationMATH 19520/51 Class 5
MATH 19520/51 Class 5 Minh-Tam Trinh University of Chicago 2017-10-04 1 Definition of partial derivatives. 2 Geometry of partial derivatives. 3 Higher derivatives. 4 Definition of a partial differential
More information13, Applications of molecular symmetry and group theory
Subject Paper No and Title Module No and Title Module Tag Chemistry 13, Applications of molecular symmetry and group theory 27, Group theory and vibrational spectroscopy: Part-IV(Selection rules for IR
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationQ SON,' (ESTABLISHED 1879L
( < 5(? Q 5 9 7 00 9 0 < 6 z 97 ( # ) $ x 6 < ( ) ( ( 6( ( ) ( $ z 0 z z 0 ) { ( % 69% ( ) x 7 97 z ) 7 ) ( ) 6 0 0 97 )( 0 x 7 97 5 6 ( ) 0 6 ) 5 ) 0 ) 9%5 z» 0 97 «6 6» 96? 0 96 5 0 ( ) ( ) 0 x 6 0
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 informationII&Ij <Md Tmlaiiiiiit, aad once in Ihe y a w Teataa m i, the vmb thatalmta oot Uiaapirit world. into as abode or wotld by them- CooTBOtioa
382 4 7 q X
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
More informationCHAPTER 2 BOOLEAN ALGEBRA
CHAPTER 2 BOOLEAN ALGEBRA This chapter in the book includes: Objectives Study Guide 2.1 Introduction 2.2 Basic Operations 2.3 Boolean Expressions and Truth Tables 2.4 Basic Theorems 2.5 Commutative, Associative,
More informationLB 220 Homework 4 Solutions
LB 220 Homework 4 Solutions Section 11.4, # 40: This problem was solved in class on Feb. 03. Section 11.4, # 42: This problem was also solved in class on Feb. 03. Section 11.4, # 43: Also solved in class
More informationMath 321 Final Exam. 1. Do NOT write your answers on these sheets. Nothing written on the test papers will be graded.
Math 31 Final Exam Instructions 1 Do NOT write your answers on these sheets Nothing written on the test papers will be graded Please begin each section of questions on a new sheet of paper 3 Do not write
More informationHigher order derivative
2 Î 3 á Higher order derivative Î 1 Å Iterated partial derivative Iterated partial derivative Suppose f has f/ x, f/ y and ( f ) = 2 f x x x 2 ( f ) = 2 f x y x y ( f ) = 2 f y x y x ( f ) = 2 f y y y
More informationA Memorial. Death Crash Branch Out. Symbol The. at Crossing Flaming Poppy. in Belding
- G Y Y 8 9 XXX G - Y - Q 5 8 G Y G Y - - * Y G G G G 9 - G - - : - G - - ) G G- Y G G q G G : Q G Y G 5) Y : z 6 86 ) ; - ) z; G ) 875 ; ) ; G -- ) ; Y; ) G 8 879 99 G 9 65 q 99 7 G : - G G Y ; - G 8
More informationi r-s THE MEMPHIS, TENN., SATURDAY. DEGfMBER
N k Q2 90 k ( < 5 q v k 3X3 0 2 3 Q :: Y? X k 3 : \ N 2 6 3 N > v N z( > > :}9 [ ( k v >63 < vq 9 > k k x k k v 6> v k XN Y k >> k < v Y X X X NN Y 2083 00 N > N Y Y N 0 \ 9>95 z {Q ]k3 Q k x k k z x X
More informationQuestion: 1. Use suitable identities to find the following products:
CH-2 Polynomial Question: 1. Use suitable identities to find the following products: (i) (x + 4) (x + 10) Solution:- (x+4)(x+10) = x 2 +10x+4x+4 x 10 = x 2 +14x+40 (ii) (x + 8) (x 10) Solution: x 2-10x+8x-80
More informationP A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.
? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>
More informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
More informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
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 information. ^e Traveler in taesnok. i the IHilty.-^ifStiiart. BbUaaoa aad WalL.""ras 'crossing a mountain»h ch w e are A«ply inteiwted. Add
x 8[ x [qqq xq F x & R FX G NR F XN R X ( F R Y
More informationg(t) = f(x 1 (t),..., x n (t)).
Reading: [Simon] p. 313-333, 833-836. 0.1 The Chain Rule Partial derivatives describe how a function changes in directions parallel to the coordinate axes. Now we shall demonstrate how the partial derivatives
More informationM E M P H I S, T E N N., S A T U E D A Y, OCTOBER 8, 1870.
5 L V 8 5 x - L : L Q ) L - \ \ Q Q - V 84 z < L L 4 Y z ( (
More information100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX
100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX.1 Norms If we have an approximate solution at a given point and we want to calculate the absolute error, then we simply take the magnitude
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationGovernor Green Triumphs Over Mudslinging
; XXX 6 928 - x 22 5 Q 0 x 2- Q- & & x 30 - x 93000000 95000000 50 000 x 0:30 7 7 2 x q 9 0 0:30 2;00 7:30 9 ( 9 & ( ( - ( - 225000 x ( ( 800 ) - 70000 200000 - x ; 200-0: 3333 0850; 778: 5-38 090; 002;
More informationIntegration of Ordinary Differential Equations
Integration of Ordinary Differential Equations Com S 477/577 Nov 7, 00 1 Introduction The solution of differential equations is an important problem that arises in a host of areas. Many differential equations
More informationCOSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods
COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationLESSON 7.1 FACTORING POLYNOMIALS I
LESSON 7.1 FACTORING POLYNOMIALS I LESSON 7.1 FACTORING POLYNOMIALS I 293 OVERVIEW Here s what you ll learn in this lesson: Greatest Common Factor a. Finding the greatest common factor (GCF) of a set of
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationDemonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES
Demonstration of the Coupled Evolution Rules 163 APPENDIX F: DEMONSTRATION OF THE COUPLED EVOLUTION RULES Before going into the demonstration we need to point out two limitations: a. It assumes I=1/2 for
More informationFORTRAN 77 Lesson 7. Reese Haywood Ayan Paul
FORTRAN 77 Lesson 7 Reese Haywood Ayan Paul 1 Numerical Integration A general all purpose Numerical Integration technique is the Trapezoidal Method. For most functions, this method is fast and accurate.
More informationUnit 2 Boolean Algebra
Unit 2 Boolean Algebra 1. Developed by George Boole in 1847 2. Applied to the Design of Switching Circuit by Claude Shannon in 1939 Department of Communication Engineering, NCTU 1 2.1 Basic Operations
More information1 Error Analysis for Solving IVP
cs412: introduction to numerical analysis 12/9/10 Lecture 25: Numerical Solution of Differential Equations Error Analysis Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationDIFFERENCE SCHEMES FOR ODE. by Michael E. Taylor
DIFFERENCE SCHEMES FOR ODE by Michael E. Taylor Contents: 1. Difference schemes based on power series. 2. Numerical integration and improved Euler methods. 3. The Runge-Kutta scheme. 1 Typeset by AMS-TEX
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationJUST THE MATHS UNIT NUMBER ORDINARY DIFFERENTIAL EQUATIONS 3 (First order equations (C)) A.J.Hobson
JUST THE MATHS UNIT NUMBER 15.3 ORDINARY DIFFERENTIAL EQUATIONS 3 (First order equations (C)) by A.J.Hobson 15.3.1 Linear equations 15.3.2 Bernouilli s equation 15.3.3 Exercises 15.3.4 Answers to exercises
More informationMathematics 222a Quiz 2 CODE 111 November 21, 2002
Student s Name [print] Student Number Mathematics 222a Instructions: Print your name and student number at the top of this question sheet. Print your name and your instructor s name on the answer sheet.
More informationSolution of First Order Initial Value Problem by Sixth Order Predictor Corrector Method
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 2277 2290 Research India Publications http://www.ripublication.com/gjpam.htm Solution of First Order Initial
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS ( WI3097 TU AESB0 ) Thursday April 6
More informationChapter 7: Exponents
Chapter 7: Exponents Algebra 1 Chapter 7 Notes Name: Algebra Homework: Chapter 7 (Homework is listed by date assigned; homework is due the following class period) HW# Date In-Class Homework Section 7.:
More information25. Chain Rule. Now, f is a function of t only. Expand by multiplication:
25. Chain Rule The Chain Rule is present in all differentiation. If z = f(x, y) represents a two-variable function, then it is plausible to consider the cases when x and y may be functions of other variable(s).
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationChapter 7: Exponents
Chapter : Exponents Algebra Chapter Notes Name: Algebra Homework: Chapter (Homework is listed by date assigned; homework is due the following class period) HW# Date In-Class Homework M / Review of Sections.-.
More informationLOWHLL #WEEKLY JOURNAL.
# F 7 F --) 2 9 Q - Q - - F - x $ 2 F? F \ F q - x q - - - - )< - -? - F - - Q z 2 Q - x -- - - - 3 - % 3 3 - - ) F x - \ - - - - - q - q - - - - -z- < F 7-7- - Q F 2 F - F \x -? - - - - - z - x z F -
More informationMATH443 PARTIAL DIFFERENTIAL EQUATIONS Second Midterm Exam-Solutions. December 6, 2017, Wednesday 10:40-12:30, SA-Z02
1 MATH443 PARTIAL DIFFERENTIAL EQUATIONS Second Midterm Exam-Solutions December 6 2017 Wednesday 10:40-12:30 SA-Z02 QUESTIONS: Solve any four of the following five problems [25]1. Solve the initial and
More informationPolynomials. In many problems, it is useful to write polynomials as products. For example, when solving equations: Example:
Polynomials Monomials: 10, 5x, 3x 2, x 3, 4x 2 y 6, or 5xyz 2. A monomial is a product of quantities some of which are unknown. Polynomials: 10 + 5x 3x 2 + x 3, or 4x 2 y 6 + 5xyz 2. A polynomial is a
More informationChem8028(1314) - Spin Dynamics: Spin Interactions
Chem8028(1314) - Spin Dynamics: Spin Interactions Malcolm Levitt see also IK m106 1 Nuclear spin interactions (diamagnetic materials) 2 Chemical Shift 3 Direct dipole-dipole coupling 4 J-coupling 5 Nuclear
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationPithy P o i n t s Picked I ' p and Patljr Put By Our P e r i p a tetic Pencil Pusher VOLUME X X X X. Lee Hi^h School Here Friday Ni^ht
G G QQ K K Z z U K z q Z 22 x z - z 97 Z x z j K K 33 G - 72 92 33 3% 98 K 924 4 G G K 2 G x G K 2 z K j x x 2 G Z 22 j K K x q j - K 72 G 43-2 2 G G z G - -G G U q - z q - G x) z q 3 26 7 x Zz - G U-
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
More informationNUMERICAL ANALYSIS SYLLABUS MATHEMATICS PAPER IV (A)
NUMERICAL ANALYSIS SYLLABUS MATHEMATICS PAPER IV (A) Unit - 1 Errors & Their Accuracy Solutions of Algebraic and Transcendental Equations Bisection Method The method of false position The iteration method
More information1. A polynomial p(x) in one variable x is an algebraic expression in x of the form
POLYNOMIALS Important Points 1. A polynomial p(x) in one variable x is an algebraic expression in x of the form p(x) = a nx n +a n-1x n-1 + a 2x 2 +a 1x 1 +a 0x 0 where a 0, a 1, a 2 a n are constants
More informationJUST THE MATHS UNIT NUMBER PARTIAL DIFFERENTIATION 1 (Partial derivatives of the first order) A.J.Hobson
JUST THE MATHS UNIT NUMBER 14.1 PARTIAL DIFFERENTIATION 1 (Partial derivatives of the first order) by A.J.Hobson 14.1.1 Functions of several variables 14.1.2 The definition of a partial derivative 14.1.3
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationHomework 1/Solutions. Graded Exercises
MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
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 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 informationoenofc : COXT&IBCTOEU. AU skaacst sftwer thsa4 aafcekr will be ehat«s«ai Bi. C. W. JUBSSOS. PERFECT THBOUGH SDFFEBISG. our
x V - --- < x x 35 V? 3?/ -V 3 - ) - - [ Z8 - & Z - - - - - x 0-35 - 3 75 3 33 09 33 5 \ - - 300 0 ( -? 9 { - - - -- - < - V 3 < < - - Z 7 - z 3 - [ } & _ 3 < 3 ( 5 7< ( % --- /? - / 4-4 - & - % 4 V 2
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 informationIOAN ŞERDEAN, DANIEL SITARU
Romanian Mathematical Magazine Web: http://www.ssmrmh.ro The Author: This article is published with open access. TRIGONOMETRIC SUBSTITUTIONS IN PROBLEM SOLVING PART IOAN ŞERDEAN, DANIEL SITARU Abstract.
More informationOptimization: Problem Set Solutions
Optimization: Problem Set Solutions Annalisa Molino University of Rome Tor Vergata annalisa.molino@uniroma2.it Fall 20 Compute the maxima minima or saddle points of the following functions a. f(x y) =
More informationCounty Council Named for Kent
\ Y Y 8 9 69 6» > 69 ««] 6 : 8 «V z 9 8 x 9 8 8 8?? 9 V q» :: q;; 8 x () «; 8 x ( z x 9 7 ; x >«\ 8 8 ; 7 z x [ q z «z : > ; ; ; ( 76 x ; x z «7 8 z ; 89 9 z > q _ x 9 : ; 6? ; ( 9 [ ) 89 _ ;»» «; x V
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 informationExistence and uniqueness: Picard s theorem
Existence and uniqueness: Picard s theorem First-order equations Consider the equation y = f(x, y) (not necessarily linear). The equation dictates a value of y at each point (x, y), so one would expect
More informationECEN 248: INTRODUCTION TO DIGITAL SYSTEMS DESIGN. Week 2 Dr. Srinivas Shakkottai Dept. of Electrical and Computer Engineering
ECEN 248: INTRODUCTION TO DIGITAL SYSTEMS DESIGN Week 2 Dr. Srinivas Shakkottai Dept. of Electrical and Computer Engineering Boolean Algebra Boolean Algebra A Boolean algebra is defined with: A set of
More informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
More informationAMC Lecture Notes. Justin Stevens. Cybermath Academy
AMC Lecture Notes Justin Stevens Cybermath Academy Contents 1 Systems of Equations 1 1 Systems of Equations 1.1 Substitution Example 1.1. Solve the system below for x and y: x y = 4, 2x + y = 29. Solution.
More informationMultivariable Calculus and Matrix Algebra-Summer 2017
Multivariable Calculus and Matrix Algebra-Summer 017 Homework 4 Solutions Note that the solutions below are for the latest version of the problems posted. For those of you who worked on an earlier version
More informationNORTH MAHARASHTRA UNIVERSITY JALGAON.
NORTH MAHARASHTRA UNIVERSITY JALGAON. Syllabus for S.Y.B.Sc. (Mathematics) With effect from June 013. (Semester system). The pattern of examination of theory papers is semester system. Each theory course
More information(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere.
MATH 4 FINAL EXAM REVIEW QUESTIONS Problem. a) The points,, ) and,, 4) are the endpoints of a diameter of a sphere. i) Determine the center and radius of the sphere. ii) Find an equation for the sphere.
More information13 Applications of molecular symmetry and group theory
Subject Chemistry Paper No and Title Module No and Title Module Tag 13 Applications of molecular symmetry and 26 and and vibrational spectroscopy part-iii CHE_P13_M26 TABLE OF CONTENTS 1. Learning Outcomes
More informationE S T A B L IS H E D. n AT Tnn G.D.O. r.w.-bal'eu. e d n e s d a y. II GRANVILLE HOUSE. GATJDICK ROAD. MEADS. EASTBOUENk
K q X k K 5 ) ) 5 / K K x x) )? //? q? k X z K 8 5 5? K K K / / $8 ± K K K 8 K / 8 K K X k k X ) k k /» / K / / / k / ] 5 % k / / k k? Z k K ] 8 K K K )» 5 ) # 8 q»)kk q»» )88{ k k k k / k K X 8 8 8 ]
More information