Two Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly
|
|
- Ethan William Goodwin
- 5 years ago
- Views:
Transcription
1 Global Journal of Pure and Applied Matematics. ISSN Volume 2, Number 2 (206), pp Researc India Publications ttp:// Two Step Hybrid Block Metod wit Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly Raed, Abdalreem Department of Pysics, College of Art and Sciences, Jordan University of Science and Tecnology. Moammad, Aldalal a Department of matematics, College of Art and Sciences, Univeristi Utara Malaysia. Abstract In tis paper, te new two step ybrid block metod wit two generalized off-step points for solving second order ordinary differential equation directly as been proposed. In te derivation of te metod, power series of order six is used as basis function to obtain te continuous sceme troug collocation and interpolation tecnique. As required by all numerical metods, te numerical properties of te developed block metod wic include consistent, zero stability, convergent and stability region are also establised. te new metod was found to compare favourably wit te existing metods in term of error. AMS subject classification: 65L05, 65L06, 65L20. Keywords: Block metod, Hybrid metod, Second order differential equation, Collocation and Interpolation, Two off step points.
2 520 Raed, Abdalreem and Moammad, Aldalal a. Introduction Tis article considered te solution to te general second order initial value problem (IVPs) of te form y f(x,y,y ), y(a) δ 0,y (a) δ.x [a,b]. () Equations () often arises in several areas of science and engineering an suc as biology, pysics and cemical. Generally, tese equations ave not exact solution. Tus, numerical metods become very crucial. In literature, numerous numerical metods for solving Equation () ave been proposed, for example, Euler metod, linear multistep metod, Runge-Kutta, predictor - corrector metod, ybrid metod and block metod (see [3], [5] and [2]). However, tese metods ave teir setbacks wic affects on teir accuracy and efficiency. Recently, ybrid block metods for solving equations () directly ave been proposed. In te latter, te researcers ave tried to combine advantages of direct, block and ybrid metods (see [], [4], [5], [6] and [7]) to overcome te zero stability problem in linear multistep as well as to avoid setbacks in reduction metods and generating numerical results concurrently. 2. Metodology In tis part, two step ybrid block metod wit two generalized off-step points i.e x n+s and x n+r for solving () is derived. Let te approximate solution of () to be te power series polomial of te form: were, y(x) q+d i0 i x [x n,x n+ ] for n 0,, 2,...,N, ( ) x i xn a i. (2) ii q denotes of te number of interpolation points wic is equal to te order of differential equation, iii d represents te number of collocation points, iv x n x n is constant step size of partition of interval [a,b] wic is given by a x 0 <x < <x N <x N b. Differentiating (2) twice gives y (x) f(x,y,y ) q+d i2 i(i ) 2 a i ( ) x i 2 xn. (3)
3 Two Step Hybrid Block Metod wit Two Generalized Off-step Points 52 Interpolating (2) at x n+s, x n+r and collocating (3) at all points in te selected interval produces seven equations wic can be written in matrix of te form: s s 2 s 3 s 4 s 5 s 6 r r 2 r 3 r 4 r 5 r s 2s 2 20s 3 30s r 2r 2 20r 3 30r a 0 a a 2 a 3 a 4 a 5 a 6 y n+s y n+r f n f n+s f n+r f n+ f n+2 (4) Employing Gaussian elimination metod to (4) gives te unknown values of a i s, i 0()6 wic are ten substituted back into equation (2) to produce a continuous implicit sceme of te form y(x) α i (x)y n+i + is,r Te first derivative of equation (5) gives were y (x) α s (x n x + r) (r s) α r (x x n s) (r s) is,r x α i(x)y n+i + i0,s,r, i0,s,r, β i (x)f n+i (5) x β i(x)f n+i (6) β 0 (x n x + s)(x n x + r) (20rs 4 ) ( 4 r 4 4 r 3 s 6 4 r 3 4 r 2 s r 2 s +0 4 r 2 4 rs rs rs + 4 s s s r 3 x 3 r 3 x n 3 r 2 sx + 3 r 2 sx n 6 3 r 2 x r 2 x n 3 rs 2 x + 3 rs 2 x n rsx rx
4 522 Raed, Abdalreem and Moammad, Aldalal a 9 3 rsx n 0 3 rx n + 3 s 3 x 3 s 3 x n 6 3 s 2 x s 2 x n sx 0 3 sx n 6 2 rx rxx n 6 2 rx 2 n + 2 s 2 x s 2 xx n + 2 s 2 x 2 n 62 sx sxx n xx n 0 2 x 2 n + rx3 3rx 2 x n + 3rxx 2 n rx3 n +sx 3 3sx 2 x n + 3sxx 2 n sx3 n 62 sx 2 n 02 x r 2 x r 2 xx n + 2 r 2 x 2 n 2 rsx rsxx n 2 rsx 2 n + 9x3 27x 2 x n +27xx 2 n 9x3 n 2x4 + 8x 3 x n 2x 2 x 2 n + 8xx3 n 2x4 n ) β s (x n x + s)(x n x + r) (60 4 s(s )(s 2)(r s)) (4 r r 3 s 6 4 r r 2 s r 2 s +0 4 r rs rs rs 2 4 s s s r 3 x 3 r 3 x n + 3 r 2 sx 3 r 2 sx n 6 3 r 2 x r 2 x n + 3 rs 2 x 3 rs 2 x n 6 3 rsx rsx n rx 0 3 rx n 2 3 s 3 x s 3 x n s 2 x 9 3 s 2 x n 0 3 sx sx n + 2 r 2 x r 2 x 2 n + 2 rsx rsxx n + 2 rsx 2 n 62 rx rxx n 6 2 rx 2 n 22 s 2 x sx sxx n +9 2 sx 2 n 02 x xx n 0 2 x 2 n + rx3 3rx 2 x n +3rxx 2 n rx3 n 2sx3 + 6sx 2 x n 6sxx 2 n + 2sx3 n + 9x3 27x 2 x n +8xx 3 n 22 r 2 xx n s 2 xx n 2 2 s 2 x 2 n + 27xx2 n 9x3 n 2x4 +8x 3 x n 2x 2 x 2 n 2x4 n ) β r (x n x + s)(x n x + r) (60 4 r(r )(r 2)(r s)) (24 r 4 4 r 3 s 9 4 r 3 4 r 2 s r 2 s r 2 4 rs rs rs 4 s s s rsx sx n 9 2 rx 2 n + 23 r 3 x 2 3 r 3 x n 3 r 2 sx + 3 r 2 sx n 9 3 r 2 x r 2 x n 3 rs 2 x + 3 rs 2 x n 6 3 rsx n rx 0 3 rx n 3 s 3 x + 3 s 3 x n s 2 x 6 3 s 2 x n 0 3 sx r 2 x r 2 xx n r 2 x 2 n 2 rsx rsxx n 2 rsx 2 n 92 rx 2 8xx 3 n 2 s 2 x s 2 xx n 2 s 2 x 2 n +6 2 sx sxx n sx 2 n + 02 x xx n x 2 n + 2rx3 6rx 2 x n + 6rxx 2 n 2rx3 n sx3 + 3sx 2 x n 3sxx 2 n + sx3 n +8 2 rxx n 9x x 2 x n 27xx 2 n + 9x3 n + 2x4 8x 3 x n +2x 2 x 2 n + 2x4 n ) β (x n x + s)(x n x + r) (60 4 (s )(r )) ( 4 r 4 4 r 3 s 4 4 r 3 4 r 2 s r 2 s 4 rs rs s s r 3 x 3 r 3 x n 3 r 2 sx + 3 r 2 sx n 4 3 r 2 x r 2 x n
5 Two Step Hybrid Block Metod wit Two Generalized Off-step Points rs 2 x + 3 rs 2 x n rsx 6 3 rsx n + 3 s 3 x 3 s 3 x n 4 3 s 2 x +4 3 s 2 x n 2 2 r 2 xx n + 2 r 2 x 2 n 2 rsx rsxx n 2 rsx 2 n 4 2 rx rxx n 4 2 rx 2 n 22 s 2 xx n + 2 s 2 x 2 n 4 2 sx sxx n 4 2 sx 2 n + rx3 3rx 2 x n + 3rxx 2 n +sx 3 3sx 2 x n + 3sxxn 2 sx3 n + 6x3 8x 2 x n + 2 s 2 x 2 rxn 3 +8xxn r 2 x 2 6xn 3 2x4 + 8x 3 x n 2x 2 xn 2 + 8xx3 n 2x4 n ) β 2 (x n x + s)(x n x + r) (20 4 ( 4 r 4 4 r 3 s 2 4 r 3 4 r 2 s 2 (s 2)(r 2)) +3 4 r 2 s 4 rs rs s s r 3 x 3 r 3 x n 3 r 2 sx + 3 r 2 sx n 2 3 r 2 x r 2 x n 3 rs 2 x + 3 rs 2 x n rsx 3 3 rsx n + 3 s 3 x 3 s 3 x n 2 3 s 2 x s 2 x n + 2 r 2 x r 2 xx n + 2 r 2 x 2 n 2 rsx rsxx n 2 rsx 2 n 22 rx rxx n 2 2 rx 2 n 22 s 2 xx n + 2 s 2 x 2 n 22 sx sxx n 2 2 sx 2 n + rx3 3rx 2 x n + 3rxx 2 n +sx 3 3sx 2 x n + 3sxx 2 n sx3 n + 3x3 9x 2 x n + 9xx 2 n 3x3 n 2x 4 + 8x 3 x n + 2 s 2 x 2 rx 3 n 2x2 x 2 n + 8xx3 n 2x4 n ) Equation (5) is evaluated at te non-interpolatig point x n+ and x n+2 wile Equation (6) is evaluated at all points to give te discrete scemes and its derivative. Te discrete sceme and its derivatives are combined i a matrix form as below [ ] A [2] 2 Y M B [2] (7) A [2] 2 r (r s) (r ) (r s) (r 2) (r s) ((r s)) ((r s)) ((r s)) ((r s)) ((r s)) s (r s) (s ) (r s) (s 2) (r s) ((r s)) ((r s)) ((r s)) ((r s)) ((r s)) ,Y M y n+s y n+r y n+ y n+2 y n+s y n+r y n+ y n+2,
6 524 Raed, Abdalreem and Moammad, Aldalal a ( y n B [2] 2 ), R [] , 2 ( ) f n, f n+s f n+r f n+ f n+2
7 Two Step Hybrid Block Metod wit Two Generalized Off-step Points (r4 r 3 s 6r 3 r 2 s 2 + 9r 2 s + 0r 2 rs 3 + 9rs 2 30rs + s 4 6s 3 + 0s 2 ) 2 2 (r )(s ) (r 4 r 3 s 5r 3 r 2 s 2 + 8r 2 s + 5r 2 rs 3 + 8rs 2 20rs) 22rs + 5r + s 4 5s 3 + 5s 2 + 5s 3) 3 2 (r 2)(s 2) (r 4 r 3 s 4r 3 r 2 s 2 + 7r 2 s + 2r 2 rs 3 + 7rs 2 20rs 4 5 6rs + 4r + s 4 4s 3 + 2s 2 + 4s) 20rs (r5 r 4 s 6r 4 r 3 s 2 + 9r 3 s + 0r 3 r 2 s 3 + 9r 2 s 2 30r 2 s rs 4 + 9rs 3 30rs 2 + s 5 6s 4 + 0s 3 ) 20rs (r5 r 4 s 6r 4 r 3 s 2 + 9r 3 s + 0r 3 r 2 s 3 + 9r 2 s 2 30r 2 s 8rs 5 + 4rs 4 2rs rs 2 8s 6 68s 5 + 9s 4 0s 3 ) 6 20rs (8r6 + 8r 5 s + 68r 5 4r 4 s 9r 4 + r 3 s 2 + 2r 3 s + 0r 3 + r 2 s 3 9r 2 s 2 30r 2 s + rs 4 9rs rs 2 s 5 + 6s 4 0s 3 ) rs (r5 r 4 s 6r 4 r 3 s 2 + 9r 3 s + 0r 3 r 2 s 3 + 9r 2 s 2 30r 2 s rs 4 + 9rs 3 30rs rs 33r + s 5 6s 4 + 0s 3 33s 59) 20rs (r5 r 4 s 6r 4 r 3 s 2 + 9r 3 s + 0r 3 r 2 s 3 + 9r 2 s 2 30r 2 s rs 4 + 9rs 3 30rs rs 576r + s 5 6s 4 + 0s 3 576s 228) ( 2 r (60(r s)(s )(s 2)) (r4 + r 3 s 6r 3 + r 2 s 2 6r 2 s + 0r 2 + rs 3 6rs 2 + 0rs 2s 4 + 9s 3 0s 2 )) ( 2 s (60(r s)(r )(r 2)) (2r4 r 3 s 9r 3 r 2 s 2 + 6r 2 s + 0r 2 rs 3 + 6rs 2 0rs s 4 + 6s 3 0s 2 )) ( 2 rs (60(r )(s )) (r4 r 3 s 4r 3 r 2 s 2 + 6r 2 s rs 3 +6rs 2 + s 4 4s 3 )) ( 2 rs (20(s 2)(r 2)) (r4 r 3 s 2r 3 r 2 s 2 + 3r 2 s rs 3 +3rs 2 + s 4 2s 3 ))
8 526 Raed, Abdalreem and Moammad, Aldalal a 2 22 ( 2 (r ) (60s(r s)(s 2)) (r4 + r 3 s 5r 3 + r 2 s 2 5r 2 s + 5r 2 + rs 3 5rs 2 + 5rs + 5r 2s 4 + 7s 3 3s 2 3s 3)) ( 2 (s ) (60r(r s)(r 2)) (2r4 r 3 s 7r 3 r 2 s 2 + 5r 2 s + 3r 2 rs 3 +5rs 2 5rs + 3r s 4 + 5s 3 5s 2 5s + 3)) 23 (2 60 (r4 r 3 s 3r 3 r 2 s 2 + 5r 2 s 3r 2 rs 3 + 5rs 2 + 5rs 3r +s 4 3s 3 3s 2 3s + 4)) 24 (2 (s )(r ) (20(r 2)(s 2)) (r4 r 3 s r 3 r 2 s 2 + 2r 2 s r 2 rs 3 +2rs 2 + 2rs r + s 4 s 3 s 2 s + )) D [2] ( 2 (r 2) 2 3 (60s(r s)(s )) (r4 + r 3 s 4r 3 + r 2 s 2 4r 2 s + 2r 2 + rs 3 4rs 2 + 2rs + 4r 2s 4 + 5s 3 )) D [2] ( 2 (s 2) 2 32 (60r(r s)(r )) (2r4 r 3 s 5r 3 r 2 s 2 + 4r 2 s rs 3 + 4rs 2 2rs s 4 + 4s 3 2s 2 4s)) 33 (2 (r 2)(s 2) (60(r )(s )) (r4 r 3 s 2r 3 r 2 s 2 + 4r 2 s 4r 2 rs 3 34 ( rs 2 + 8rs 8r + s 4 2s 3 4s 2 8s + 6)) 20 (r4 r 3 s r 2 s 2 + r 2 s rs 3 + rs 2 + 2rs + s 4 8)) (60s(r s)(s )(s 2)) (r5 + r 4 s 6r 4 + r 3 s 2 6r 3 s + 0r 3 +r 2 s 3 6r 2 s 2 + 0r 2 s + rs 4 6rs 3 + 0rs 2 2s 5 + 9s 4 0s 3 )) (60r(r s)(r )(r 2)) (2r5 r 4 s 9r 4 r 3 s 2 + 6r 3 s + 0r 3 r 2 s 3 +6r 2 s 2 0r 2 s rs 4 + 6rs 3 0rs 2 s 5 + 6s 4 0s 3 ) (60(r )(s )) r5 r 4 s 4r 4 r 3 s 2 + 6r 3 s r 2 s 3 + 6r 2 s 2 rs 4 + 6rs 3 + s 5 4s 4 ) (20(s 2)(r 2)) (r5 r 4 s 2r 4 r 3 s 2 + 3r 3 s r 2 s 3 + 3r 2 s 2 rs 4 +3rs 3 + s 5 2s 4 ))
9 Two Step Hybrid Block Metod wit Two Generalized Off-step Points (60s(r s)(s )(s 2)) (r5 + r 4 s 6r 4 + r 3 s 2 6r 3 s + 0r 3 + r 2 s 3 6r 2 s 2 + 0r 2 s 8rs 5 4rs rs 3 50rs 2 56s 5 36s s 3 ) (60r(r s)(r )(r 2)) (2r5 r 4 s 9r 4 r 3 s 2 + 6r 3 s + 0r 3 r 2 s 3 +6r 2 s 2 0r 2 s rs 4 + 6rs 3 0rs 2 + 8s s 5 9s 4 + 0s 3 ) (60(r )(s )) (r5 r 4 s 4r 4 r 3 s 2 + 6r 3 s r 2 s 3 + 6r 2 s 2 8rs 5 +4rs 4 4rs 3 8s 6 50s 5 + 6s 4 ) ) (20(s 2)(r 2) (r5 r 4 s 2r 4 r 3 s 2 + 3r 3 s r 2 s 3 + 3r 2 s 2 8rs 5 +4rs 4 7rs 3 8s 6 32s 5 + 3s 4 ) (60s(r s)(s )(s 2)) (8r6 + 68r 5 r 4 s 9r 4 r 3 s 2 + 6r 3 s + 0r 3 r 2 s 3 + 6r 2 s 2 0r 2 s rs 4 + 6rs 3 0rs 2 + 2s 5 9s 4 + 0s 3 ) ) (60r(r s)(r )(r 2)) (8r5 s + 56r 5 + 4r 4 s + 36r 4 r 3 s 2 54r 3 s 30r 3 r 2 s 3 + 6r 2 s r 2 s rs 4 + 6rs 3 0rs 2 s 5 +6s 4 0s 3 ) (60(r )(s )) (8r6 + 8r 5 s + 50r 5 4r 4 s 6r 4 + r 3 s 2 + 4r 3 s +r 2 s 3 6r 2 s 2 + rs 4 6rs 3 s 5 + 4s 4 ) (20(s 2)(r 2)) (8r6 + 8r 5 s + 32r 5 4r 4 s 3r 4 + r 3 s 2 + 7r 3 s +r 2 s 3 3r 2 s 2 + rs 4 3rs 3 s 5 + 2s 4 ) (60s(r s)(s )(s 2)) (r5 + r 4 s 6r 4 + r 3 s 2 6r 3 s + 0r 3 + r 2 s 3 6r 2 s 2 + 0r 2 s + rs 4 6rs 3 + 0rs 2 33r 2s 5 + 9s 4 0s 3 59) (60r(r s)(r )(r 2)) (2r5 r 4 s 9r 4 r 3 s 2 + 6r 3 s + 0r 3 r 2 s 3 +6r 2 s 2 0r 2 s rs 4 + 6rs 3 0rs 2 s 5 + 6s 4 0s s + 59) (60(r )(s )) (r5 r 4 s 4r 4 r 3 s 2 + 6r 3 s r 2 s 3 + 6r 2 s 2 rs 4 +6rs 3 40rs + 7r + s 5 4s 4 + 7s 66)
10 528 Raed, Abdalreem and Moammad, Aldalal a (20(s 2)(r 2)) (r5 r 4 s 2r 4 r 3 s 2 + 3r 3 s r 2 s 3 + 3r 2 s 2 rs 4 +3rs 3 0rs 3r + s 5 2s 4 3s 33) (60s(r s)(s )(s 2)) (r5 + r 4 s 6r 4 + r 3 s 2 6r 3 s + 0r 3 + r 2 s 3 6r 2 s 2 + 0r 2 s + rs 4 6rs 3 + 0rs 2 576r 2s 5 + 9s 4 0s 3 228) (60r(r s)(r )(r 2)) (2r5 r 4 s 9r 4 r 3 s 2 + 6r 3 s + 0r 3 r 2 s 3 +6r 2 s 2 0r 2 s rs 4 + 6rs 3 0rs 2 s 5 + 6s 4 0s s + 228) (60(r )(s )) (r5 r 4 s 4r 4 r 3 s 2 + 6r 3 s r 2 s 3 + 6r 2 s 2 rs 4 +6rs 3 80rs 496r + s 5 4s 4 496s 632) (20(s 2)(r 2)) (r5 r 4 s 2r 4 r 3 s 2 + 3r 3 s r 2 s 3 + 3r 2 s 2 rs 4 +3rs rs 656r + s 5 2s 4 656s 86) Multiplying Equation (7) by te inverse of A [] 2 gives IY m B [ ] [2] (8) were I is 8 8 identity matrix and B [] 2 s r , , Ē []
11 Two Step Hybrid Block Metod wit Two Generalized Off-step Points 529 wit (2 s 2 (40r 0s 5rs + 2rs 2 + 6s 2 s 3 )) (20r) 2 (2 r 2 (0r 40s + 5rs 2r 2 s 6r 2 + r 3 )) (20s) 3 (2 (35rs 8s 8r + 3)) (20rs) 4 (22 (5rs s r)) (5rs) 5 (s(20s 60r + 30rs 5rs2 + 8rs 3 5s s 3 + 8s 4 )) (20r) 6 (r(20r 60s + 30rs 5r2 s + 8r 3 s 5r r 3 + 8r 4 )) (20s) 7 8 ((50rs 33s 33r 59)) 20rs) ((5rs 72s 72r 266)) (5rs) (2 s 2 (20r 0s 5rs + 3rs 2 + 9s 2 2s 3 )) (60(s )(s 2)(r s)) 2 ( 2 s 4 (s 2 6s + 0)) (60r(r )(r 2)(r s)) 3 (2 s 4 (4s 0r + 2rs s 2 )) (60(r )(s )) 4 (2 s 4 (2s 5r + 2rs s 2 )) (20(r 2)(s 2)) 2 ( 2 r 4 (r 2 6r + 0)) (60s(s )(s 2)(r s)) 22 (2 r 2 (0r 20s + 5rs 3r 2 s 9r 2 + 2r 3 )) (60(r )(r 2)(r s)) 23 (2 r 4 (0s 4r 2rs + r 2 )) (60(s )(r )) 24 (2 r 4 (5s 2r 2rs + r 2 )) (20(s 2)(r 2))
12 530 Raed, Abdalreem and Moammad, Aldalal a 3 32 ( 2 (8r 3)) (60s(s )(s 2)(r s)) ( 2 (8s 3)) (60r(r )(r 2)(r s)) 33 (2 (5rs 7s 7r + 4)) (60(s )(r )) 34 (2 (5rs 2s 2r + )) (20(s 2)(r 2)) 4 42 (4 2 r) (5s(s )(s 2)(r s)) (4 2 s) (5r(r )(r 2)(r s)) 43 (42 (5rs 4s 4r + 4)) (5(s )(r )) 44 (22 (r + s 2)) (5(s 2)(r 2)) 5 (s(60r 40s 60rs + 5rs2 + 8rs s s 3 )) (60(s )(s 2)(r s)) 52 (s3 (8s s 2 5s + 20)) (60r(r )(r 2)(r s)) 53 (s3 (20r 0s 5rs + 8rs 2 + 5s 2 + 8s 3 )) (60(s )(r )) 54 (s3 (0r 5s 5rs + 8rs s 2 + 8s 3 )) (20(s 2)(r 2)) 6 (r3 (8r r 2 5r + 20)) (60s(s )(s 2)(r s)) 62 (r(60s 40r 60rs + 5r2 s + 8r 3 s + 45r r 3 )) (60(r )(r 2)(r s)) 63 (r3 (20s 0r 5rs + 8r 2 s + 5r 2 + 8r 3 )) (60(s )(r ))
13 Two Step Hybrid Block Metod wit Two Generalized Off-step Points (r3 (0s 5r 5rs + 8r 2 s + 33r 2 + 8r 3 )) (20(s 2)(r 2)) ((33r + 59)) (60s(s )(s 2)(r s)) ((33s + 59)) (60r(r )(r 2)(r s)) ((40rs 7s 7r + 66)) (60(s )(r )) ((3r + 3s + 0rs + 33)) (20(s 2)(r 2)) (4(36r + 33)) (5s(s )(s 2)(r s)) (4(36s + 33)) (5r(r )(r 2)(r s)) (4(3r + 3s + 5rs + 02)) (5(s )(r )) ((5rs 82s 82r 02)) (5(s 2)(r 2)) 3. Analysis of te Metod 3.. Order of te Metod Te linear difference operator L associated wit (8) is defined as L[y(x); ] IY M B [ ] [2] (9) were y(x) is an arbitrary test function continuously differentiable on [a,b]. Y M and 3 components are expanded in Taylor s series respectively and its terms are collected in powers of to give L[y(x),] C 0 y(x) + C y (x) + C 2 y (x) + (0) Definition 3.. Hybrid block metod (8) and associated linear operator (9) are said to be of order p, ifc 0 C 2 C p+2 0 and C p+2 0 wit error vector constants C p+2.
14 532 Raed, Abdalreem and Moammad, Aldalal a Expanding (8) in Taylor series about x n gives (s) j j j y n (s)y n (s2 (40r 0s 5rs + 2rs 2 + 6s 2 s 3 )) (20r) j0 (s2 (20r 0s 5rs + 3rs 2 + 9s 2 2s 3 )) (s) j (60(s )(s 2)(r s)) j0 (s 4 (s 2 6s + 0)) (60r(r )(r 2)(r s)) j0 + (s4 (4s 0r + 2rs s 2 )) (r) j (60(r )(s )) j0 j0 j (s4 (2s 5r + 2rs s 2 )) (20(r 2)(s 2)) j0 y n yj n y n ()y n (r2 (0r 40s + 5rs 2r 2 s 6r 2 + r 3 )) (20s) (r 4 (r 2 6r + 0)) (s) j (60s(s )(s 2)(r s)) j0 (r2 (0r 20s + 5rs 3r 2 s 9r 2 + 2r 3 )) (60(r )(r 2)(r s)) j0 (r4 (0s 4r 2rs + r 2 )) (r) j (60(s )(r )) (r4 (5s 2r 2rs + r 2 )) (20(s 2)(r 2)) j0 (r) j j j y n (r)y ((35rs 8s 8r + 3)) n (20rs) j0 (r 4 (r 2 6r + 0)) (s) j (60s(s )(s 2)(r s)) j0 ( 2 (8s 3)) (60r(r )(r 2)(r s)) j0 (2 (5rs 7s 7r + 4)) (r)j (60(s )(r )) j0 (2 (5rs 2s 2r + )) (20(s 2)(r 2)) j0 (2) j j j y n (2)y n (22 (5rs s r)) y n (5rs) j0 (4 2 r) (s) j (5s(s )(s 2)(r s)) j0 (4 2 s) (5r(r )(r 2)(r s)) j0 (42 (5rs 4s 4r + 4)) (r)j (22 (r + s 2)) (5(s )(r )) (5(s 2)(r 2)) j0 j0 y n j0 y n y n y n
15 Two Step Hybrid Block Metod wit Two Generalized Off-step Points 533 By comparing te coefficient of, we obtain te order of te metod to be [5, 5, 5, 5] T wit error constant vector s 4 (70r 28s 42rs + 7rs 2 + 2s 2 4s 3 ) r 4 (70s 28r 42rs + 7r 2 s + 2r r 3 ) C (56rs 2s 2r + 60) (7rs + 396) Zero Stability Te ybrid block metod (8) is said to be zero stable if te first caracteristic polomial π(x) aving roots suc tat x z, and if x z, ten, te multiplicity of x z must not greater tan two. In order to find te zero-stability of main block ˆB [2] 2 [y n+s,y n+r,y n+,y n+2 ] T in (8), we only consider te solution of te first caracteristic polomial were I is 4 4 identity matrix, tat is (r) xi ˆB [2] 2 x x 3 (x ) wic implies x 0, 0,,. Hence, our metod is zero stable for all s, r (0, ) 3.3. Consistency Te one step ybrid block metod (8) is said to be consistent if its order greater tan or equal one i.e. P. Tis proves tat our metod is consistent for all s, r (0, ) Convergence Teorem 3.2. [Henrici, 962] Consistency and zero stability are sufficient conditions for a linear multistep metod to be convergent. Since te metod is consistent and zero stable, it implies te metod is convergent for all s, r (0, ) Numerical Results In finding te accuracy of our metods, te following second order ODEs are examined. Te new block metods solved te same problems te existing metods solved in order to compare results in terms of error.
16 534 Raed, Abdalreem and Moammad, Aldalal a Table : Comparison of te new metod wit Kayode and Adeyeye (203) for solving Problem, 0 x exact solution computed solution in new metod error in our metod, s 4, r 3 errors in [2] e e e e e e e e e e 5 Table 2: Comparison of te new metod wit Awoyemi et al. (20) for solving problem 2, 320 x exact solution computed solution in new metod s error in our metod,s 3 4, r 3 errors in [3] e e e e e e e e e e 8 Problem 2: y y 0, y(0) 0, y (0), 0.. Exact solution: y(x) e x Problem : y x(y ) 2 0, y(0), y (0) 2, 320. Exact solution: y(x) + ( ) 2 + x 2 ln 2 x 4. Conclusion A two step ybrid block metod wit two generalized off-step points was developed. Te metod was tested to be convergent wit order five for all general off-step point belong to selected interval. Te derived metod was applied to solve bot non-linear and linear second ODEs problems witout converting to te equivalents system of first order ODEs. New generated results confirm te accuracy of te new metods in terms of error.
17 Two Step Hybrid Block Metod wit Two Generalized Off-step Points 535 References [] S. N. Jator, Solving second order initial value problems by a ybrid multistep metod witout predictors, Applied Matematics and Computation, (200), 27(8), [2] D. O.Awoyemi, E.A.Adebile, A. O.Adesanya & T.A.Anake, Modified block metod for te direct solution of second order ordinary differential equations, International Journal of Applied Matematics and Computation, (20)., 3(3), [3] Kayode, & O. Adeyey, Two-step two-point ybrid metods for general second order differential equations, African Journal of Matematics and Computer Science Researc, (203), 6(0), [4] T. A. Anake, D. O. Awoyemi and A. A. Adesanya, A one step metod for te solution of general second order ordinary differential equations, International Journal of science and Tecnology, (202b), 2(4), [5] J.D. Lambert, Computational metods in ordinary differential equations, (973). [6] T. A Anake, D. O. Awoyemi and A. A. Adesanya, One-step implicit ybrid block metod for te direct solution of general second order ordinary differential equations, IAENG International Journal of Applied Matematics, (202a), 42(4), [7] A. Sagir, An accurate computation of block ybrid metod for solving stiff ordinary differential equations, Journal of Matematics, (202), 4, 8 2. [8] Henrici, P. (962). Discrete variable metods in ordinary differential equations.
18
One-Step Hybrid Block Method with One Generalized Off-Step Points for Direct Solution of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 10, 2016, no. 29, 142-142 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.6127 One-Step Hybrid Block Method with One Generalized Off-Step Points for
More informationZ. Omar. Department of Mathematics School of Quantitative Sciences College of Art and Sciences Univeristi Utara Malaysia, Malaysia. Ra ft.
International Journal of Mathematical Analysis Vol. 9, 015, no. 46, 57-7 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.015.57181 Developing a Single Step Hybrid Block Method with Generalized
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationContinuous Block Hybrid-Predictor-Corrector method for the solution of y = f (x, y, y )
International Journal of Mathematics and Soft Computing Vol., No. 0), 5-4. ISSN 49-8 Continuous Block Hybrid-Predictor-Corrector method for the solution of y = f x, y, y ) A.O. Adesanya, M.R. Odekunle
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 informationOn One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations
Pure and Applied Matematics Journal 7; 6(5: 74 ttp://wwwsciencepublisinggroupcom/j/pamj doi: 648/jpamj765 ISSN: 6979 (Print; ISSN: 698 (Online On One Justiication on te Use o Hybrids or te Solution o First
More informationAn Accurate Self-Starting Initial Value Solvers for. Second Order Ordinary Differential Equations
International Journal of Contemporar Matematical Sciences Vol. 9, 04, no. 5, 77-76 HIKARI Ltd, www.m-iari.com ttp://dx.doi.org/0.988/icms.04.4554 An Accurate Self-Starting Initial Value Solvers for Second
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationBlock Algorithm for General Third Order Ordinary Differential Equation
ICASTOR Journal of Mathematical Sciences Vol. 7, No. 2 (2013) 127-136 Block Algorithm for General Third Order Ordinary Differential Equation T. A. Anake 1, A. O. Adesanya 2, G. J. Oghonyon 1 & M.C. Agarana
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationModified block method for the direct solution of second order ordinary differential equations
c Copyright, Darbose International Journal of Applied Mathematics and Computation Volume 3(3),pp 181 188, 2011 http://ijamc.psit.in Modified block method for the direct solution of second order ordinary
More informationContinuous Hybrid Multistep Methods with Legendre Basis Function for Direct Treatment of Second Order Stiff ODEs
American Journal o Computational and Applied Matematics 06, 6(): 8-9 DOI: 0.59/j.ajcam.06060.0 Continuous Hybrid Multistep Metods wit Legendre Basis Function or Direct Treatment o Second Order Sti ODEs
More informationMA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE
MA64 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS B Dr.T.Kulandaivel Department of Applied Matematics SVCE Numerical ordinar differential equations is te part
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationRunge Kutta Collocation Method for the Solution of First Order Ordinary Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 1, 17-26 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.5823 Runge Kutta Collocation Method for the Solution of First
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationOne Step Continuous Hybrid Block Method for the Solution of
ISSN 4-86 (Paper) ISSN -9 (Online) Vol.4 No. 4 One Step Continuous Hybrid Block Method for the Solution of y = f ( x y y y ). K. M. Fasasi * A. O. Adesanya S. O. Adee Department of Mathematics Modibbo
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationUniform Order Legendre Approach for Continuous Hybrid Block Methods for the Solution of First Order Ordinary Differential Equations
IOSR Journal of Mathematics (IOSR-JM) e-issn: 8-58, p-issn: 319-65X. Volume 11, Issue 1 Ver. IV (Jan - Feb. 015), PP 09-14 www.iosrjournals.org Uniform Order Legendre Approach for Continuous Hybrid Block
More informationRunge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More information(5.5) Multistep Methods
(5.5) Mulstep Metods Consider te inial-value problem for te ordinary differenal equaon: y t f t, y, a t b, y a. Let y t be te unique soluon. In Secons 5., 5. and 5.4, one-step numerical metods: Euler Metod,
More informationSolving Second Order Linear Dirichlet and Neumann Boundary Value Problems by Block Method
IAENG International Journal o Applied Matematics 43: IJAM_43 04 Solving Second Order Linear Diriclet and Neumann Boundar Value Problems b Block Metod Zanaria Abdul Majid Mod Mugti Hasni and Norazak Senu
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
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 informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationCONVERGENCE PROPERTIES OF THE MODIFIED NUMEROV RELATED BLOCK METHODS
Vol 4 No 1 January 213 CONVERGENCE PROPERTIES OF THE MODIFIED NUMEROV RELATED BLOCK METHODS YS, Awari 1 and AA, Abada 2 1 Department of Mathematics & 2 Department of Statistics, BinghamUniversity, Karu,
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationCollege of Art and Sciences Universiti Utara Sintok Kedah, 0060, MALAYSIA
Iteratioal Joural of Pure ad Applied Mathematics Volume 2 No. 3 207, 497-57 ISSN: 3-8080 (prited versio); ISSN: 34-3395 (o-lie versio) url: http://www.ijpam.eu doi: 0.2732/ijpam.v2i3.5 PAijpam.eu GENERALIZED
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationPacking polynomials on multidimensional integer sectors
Pacing polynomials on multidimensional integer sectors Luis B Morales IIMAS, Universidad Nacional Autónoma de México, Ciudad de México, 04510, México lbm@unammx Submitted: Jun 3, 015; Accepted: Sep 8,
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationAn Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method
Int. J. Contemp. Mat. Sciences, Vol. 2, 27, no. 2, 983-989 An Approximation to te Solution of te Brusselator System by Adomian Decomposition Metod and Comparing te Results wit Runge-Kutta Metod J. Biazar
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationDe-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationA Zero-Stable Block Method for the Solution of Third Order Ordinary Differential Equations.
A Zero-Stable Block Method for the Solution of Third Order Ordinary Differential Equations. K. Rauf, Ph.D. 1 ; S.A. Aniki (Ph.D. in view) 2* ; S. Ibrahim (Ph.D. in view) 2 ; and J.O. Omolehin, Ph.D. 3
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationVariance Estimation in Stratified Random Sampling in the Presence of Two Auxiliary Random Variables
International Journal of Science and Researc (IJSR) ISSN (Online): 39-7064 Impact Factor (0): 3.358 Variance Estimation in Stratified Random Sampling in te Presence of Two Auxiliary Random Variables Esubalew
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationOn Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
Nonlinear Analysis and Differential Equations, Vol. 6, 08, no., 5-64 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.08.8 On Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More informationPAijpam.eu NEW SELF-STARTING APPROACH FOR SOLVING SPECIAL THIRD ORDER INITIAL VALUE PROBLEMS
International Journal of Pure and Applied Mathematics Volume 118 No. 3 218, 511-517 ISSN: 1311-88 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 1.12732/ijpam.v118i3.2
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationResearch Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation
Advances in Numerical Analysis Volume 204, Article ID 35394, 8 pages ttp://dx.doi.org/0.55/204/35394 Researc Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic
More informationModule 4 Linear Boundary value problems
Module 4 Linear Boundary value problems Lecture Content Hours 1 Finite Difference Metods: Diriclet type boundary 1 condition Finite Difference Metods: Mixed boundary condition 1 3 Sooting Metod 1 4 Sooting
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More information1 Lecture 13: The derivative as a function.
1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More information160 Chapter 3: Differentiation
3. Differentiation Rules 159 3. Differentiation Rules Tis section introuces a few rules tat allow us to ifferentiate a great variety of functions. By proving tese rules ere, we can ifferentiate functions
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a
More informationPECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER
Materials Pysics and Mecanics (5) 5- Received: Marc 7, 5 PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER Т.I. Belyankova *, V.V. Kalincuk Soutern Scientific Center of Russian
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More information-Stable Second Derivative Block Multistep Formula for Stiff Initial Value Problems
IAENG International Journal of Applied Mathematics, :3, IJAM 3_7 -Stable Second Derivative Bloc Multistep Formula for Stiff Initial Value Problems (Advance online publication: 3 August ) IAENG International
More informationResearch Article Discrete Mixed Petrov-Galerkin Finite Element Method for a Fourth-Order Two-Point Boundary Value Problem
International Journal of Matematics and Matematical Sciences Volume 2012, Article ID 962070, 18 pages doi:10.1155/2012/962070 Researc Article Discrete Mixed Petrov-Galerkin Finite Element Metod for a Fourt-Order
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationResearch Article An Accurate Block Hybrid Collocation Method for Third Order Ordinary Differential Equations
Applied Mathematics, Article ID 549597, 9 pages http://dx.doi.org/1.1155/14/549597 Research Article An Accurate Bloc Hybrid Collocation Method for Third Order Ordinary Differential Equations Lee Ken Yap,
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationLecture 2: Symplectic integrators
Geometric Numerical Integration TU Müncen Ernst Hairer January February 010 Lecture : Symplectic integrators Table of contents 1 Basic symplectic integration scemes 1 Symplectic Runge Kutta metods 4 3
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
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 informationDerivatives. if such a limit exists. In this case when such a limit exists, we say that the function f is differentiable.
Derivatives 3. Derivatives Definition 3. Let f be a function an a < b be numbers. Te average rate of cange of f from a to b is f(b) f(a). b a Remark 3. Te average rate of cange of a function f from a to
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationESTIMATION OF A POPULATION MEAN OF A SENSITIVE VARIABLE IN STRATIFIED TWO-PHASE SAMPLING
Pak. J. Stati. 06 Vol. 3(5), 393-404 ESTIMATION OF A POPUATION MEAN OF A SENSITIVE VARIABE IN STRATIFIED TWO-PHASE SAMPING Nadia Muaq, Muammad Noor-ul-Amin and Muammad Hanif National College of Business
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationProblem Solving. Problem Solving Process
Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and
More informationSKP Engineering College
S.K.P. Engineering College, Tiruvannamalai SKP Engineering College Tiruvannamalai 666 A Course Material on By S.Danasekar Assistant Professor Department of Matematics S.K.P. Engineering College, Tiruvannamalai
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4. Strict Convexity, Smootness, and Gateaux Di erentiablity Definition 4... Let X be a Banac space wit a norm denoted by k k. A map f : X \{0}!X \{0}, f 7! f x is called
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationMath 242: Principles of Analysis Fall 2016 Homework 7 Part B Solutions
Mat 22: Principles of Analysis Fall 206 Homework 7 Part B Solutions. Sow tat f(x) = x 2 is not uniformly continuous on R. Solution. Te equation is equivalent to f(x) = 0 were f(x) = x 2 sin(x) 3. Since
More information