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 Methods for the Solution of First Order Ordinary Differential Equations N. S. Yakusak *, S. Emmanuel and M. O. Ogunniran Department of Mathematics Federal University Lokoja, Kogi State, Nigeria Abstract: We adopted the method of interpolation of the approximation and collocation of its differential equation and with Legendre polynomial of the first kind as basis function to yield a continuous Linear Multistep Method with constant step size. The methods are verified to be consistent and satisfies the stability condition. Our methods was tested on first order Ordinary Differential Equation (ODE) and found to give better result when compared with the analytical solution. Keywords: Collocations, Legendre polynomial, Interpolation, Continuous scheme. I. Introduction We consider a numerical method for solving general first order Initial Value Problems (IVPs) of Ordinary Differential Equations (ODEs) of the form: y = f(x, y), y x 0 = y 0 (1) Where f is a continuous function and satisfies Lipschitz condition of the existence and uniqueness of solution. A differential equation in which the unknown function is a function of two or more independent variable is called a Partial Differential Equations (PDEs). Those in which the unknown function is function of only one independent variable are called Ordinary Deferential Equations (ODEs). Many scholars have worked extensively on the solution of (1) in literatures [1-6]. These authors proposed different method ranging from predictor corrector method to block method using different polynomials as basis functions, evaluated at some selected points. In this paper we proposed three-step and four-step hybrid block method with two off- step points, using Legendre polynomials evaluated at grids and off-grids points to give a discrete scheme. II. Derivation Of The Method In this section, we intend to develop the Linear Multistep Method (LMM), by interpolating and collocating at some selected points. We consider a Legendre approximation of the form; s+r 1 y x = a j p j x () j =0 Where r and s are interpolation and collocation point Substituting (3) into (1), we have s+r 1 y x = ja j p j 1 x (3) j =0 s+r 1 f(x, y) = ja j p j 1 x (4) j =0.1 Three-Step Method With Two Off-Step Points. Interpolating () at x n and collocating (4) at x n+s, s = 0,1, 3,, 5, 3 give the system of polynomial equation in the form. AX = U (5) DOI: 10.990/58-11140914 www.iosrjournals.org 9 Page
Solving the above matrix by Gaussian elimination method we obtained the following results, a 0 = y n + 3 8 hf n +3 18 35 hf n+ 5 58 35 hf n+ 3 + 1 80 hf n + 11 0 hf n +1 + 459 80 hf n + a 1 = 3 56 hf n + 3 56 hf n +3 1 35 hf n + 3 + 43 80 hf n +1 + 43 80 hf n + + 459 80 hf n + a = 1 56 hf n +3 + 18 35 hf n+ 5 + 6 hf n + 3 19 80 hf n 45 56 hf n +1 56 hf n + a 3 = 1 16 hf n + 1 16 hf n +3 4 5 hf n+ 3 a 4 = 11 6160 hf n +3 + 54 385 hf n+ 5 + 18 hf n+ 3 61 6160 hf n + 351 6160 hf n +1 511 6160 hf n + + 80 hf n +1 + 80 hf n + a 5 = 3 140 hf n +3 + 1 35 hf n+ 3 + 3 140 hf n 140 hf n +1 140 hf n + a 6 = 9 308 hf n +3 54 385 hf n + 5 18 hf n+ 3 9 1540 hf n + 308 hf n +1 + 81 308 hf n + (6) Substituting (6) into () we obtained, the LMM as, 3 y z = α 0 z y n + j=0 β 0 (z)f n+1 + β3 z f n+ 3 + β5(z)f 5 n+ where α 0 z and β 0 (z) are continuous coefficients obtained as α 0 z = 1 β 0 z = z 31 z4 9 0 z + 9 z3 + 4 45 z5 1 135 z6 β3 β 1 z = 119 4 z4 + 15 z 19 z3 6 5 z5 + 1 9 z6 z = 104 9 z4 40 3 z + 536 z3 + 136 45 z5 8 z6 β z = 91 8 z4 + 45 4 z 18z 3 16 5 z5 + 1 3 z6 β5 z = 4 5 z 16 3 z4 + 8z 3 8 45 z6 + 8 5 z5 β 3 z = 1 z4 + 5 6 z 54 z3 14 45 z5 + 1 z6 () Where z = x x n h Equation () is known as the continuous scheme. Evaluating () at x n +s, s = 1, 1,, 3, 3 we obtained the following discrete scheme DOI: 10.990/58-11140914 www.iosrjournals.org 10 Page
y n+1 = y n + h 11 40 f n + 63 360 f n +1 104 45 f n+ 3 + 11 10 f n + 3 45 f n + 5 + 43 360 f n +3 y 3 n+ = y n + h 35 18 f n + 133 640 f n +1 40 f n + 3 + 1053 640 f n + 40 f n+ 5 + 3 640 f n +3 y n+ = y n + h 3 135 f n + 9 45 f n +1 4 135 f n + 3 + 9 15 f n + 3 45 f n+ 5 + 16 145 f n +3 y 5 n+ = y n + h 35 18 f n + 35 115 f n +1 15 f n+ 3 + 85 384 f n + 35 f n + 5 + 15 115 f n +3 y n+3 = y n + h 11 40 f n + 81 40 f n +1 8 5 f n+ 3 + 81 40 f n + + 11 40 f n +3 (8) Equation (8) is known as the required block method.. Four-Step Method With Two Off-Step Points Following similar procedure above, we interpolate () at x n and collocation (4) at x n +s, s =0, 1, 5,,, 4 and Evaluating at different grids and off-grid point, we obtain the method as y n+1 = y n + h 039 056 f n + 04 151 f n +1 958 315 f n + + 4636 946 f n + 5 3134 05 f n + 3359 1510 f n +4 y n+ = y n + h 83 94 f n + 1598 945 f n +1 1 105 f n + + 339 945 f n+ 5 30 105 f n +3 + 83 38 f n + 341 1890 f n +4 y n+3 = y n + h 110 390 f n + 43 80 f n +1 51 35 f n + + 148 35 f n + 5 13 89 f n +3 + 6 45 f n+ 101 560 f n +4 68 f n + 5831 3456 f n +1 343 40 f n + + 4459 4459 190 f n +3 + 161 10 f n+ 651 34560 f n +4 y n+ = y n + h 1 y 5 n+ = y n + h 1080 f n+ 5 31895 11896 f n + 4085 419 f n +1 85 3016 f n + + 5965 151 f n+ 5 385 8064 f n +3 + 405 358 f n+ 885 48384 f n +4 (9) (10) y n+4 = y n + h 6 f 05 n + 30 f 189 n +1 4 f 315 n + + 4096 f 946 n + 5 83 f 315 n +3 + 4096 f 06 n + Equation (9) and (10) together form the block method. 6 946 f n +4 III. Analysis Of The Methods. In this section we discuss the local truncation error and order, consistency and zero stability of the scheme generated. 3.1 Order And Error Constant Let the linear operator L y x ; h associated with the block formula be as k L y x ; h = j =0 α j y x n + jh hβ j y x n + jh (3.1) Expanding in Taylor series expansion and comparing the coefficients of h gives L y x ; h = C 0 y x + C 1 hy x + C 1 h y x C p h p y p x + C p+1 h p+1 y p+1 x Definition 3.1 The Linear operator L and the associated continuous LMM (3.1) are said to be of order p if C 0 = C 1 = C = C 3 = C 3 = C p = 0, C p+1 0 is called the error constant. DOI: 10.990/58-11140914 www.iosrjournals.org 11 Page
Table.1: Features of the block method (8) Order y n +1 y n+ 3 y n + y n+ 5 y n +3 Error Constant 4 419 14336 9 1510 5 380 9 4480 Table.: Features of the method (10) and (11) Order Error constant y n +4 16 6615 y n +3 31 1544 y n + 9 1160 y n +1 965 33860 y 5 n + 105 4334064 4 y n + 3. Consistent And Zero Stability Definition 3. The LMM is said to be consistent if it has order P 1. Definition 3.3 The block method is said to be zero stable if the roots Z s, s = 1,,, N of the characteristic p(z) defined by p Z = det ZA 0 A Satisfied Z s 1 have the multiplicity not exceeding the order of the differential r μ (z 1)μ equation, as h 0, p Z = Z Where μ is the order of the differential equation r is the order of the matrix A 0 and A []. Putting (8) in matrix form we obtain 98304 DOI: 10.990/58-11140914 www.iosrjournals.org 1 Page
Normalizing the matrix as h 0 ρ Z = [ZA 0 A ] Z 4 Z 1 = 0 Z = 0, Z = 1 The block method (8) is observed to be zero-stable. Equation (9) and (10) in matrix form, and following the same procedure above is found to be zero-stable. [8] 3.3 Convergence The convergence of the continuous hybrid block method is considered in the light of the basic properties discussed above in conjunction with the fundamental theorem of Dahlquist [8] for LMM; we state the Dahlquist theorem without proof. Theorem 3.1 The necessary and sufficient condition for a linear multistep method to be convergent is for it to be consistent and zero stable. Following the theorem 3.1 above shows that both methods are convergent. IV. Numerical Examples We now implement our derived block methods, on first order initial value problems. In order to test the efficiency of the methods, we employed the following notations in our tables below: 3SHBM: 3 Step two off - step Hybrid Block Method. 4SHBM: 4 Step two off - step Hybrid Block Method. Example 4.1 y = 5y, y 0 = 1, h = 0.01 Exact solution y(x) = e 5x Example 4. y + y = 0, y 0 = 1, h = 0.1 Exact solution y x = e x Table 4.1: Numerical results for example 4.1 X EXACT 3SHBM 4SHBM 0.0 1.000000000 1.000000000 1.000000000 0.01 1.0511096 1.0511096 1.0511096 0.0 1.10510918 1.10510918 1.10510918 0.03 1.16183443 1.16183443 1.16183443 0.04 1.14058 1.14058 1.14058 0.05 1.840541 1.840541 1.840541 0.06 1.349858808 1.349858808 1.34985880 0.0 1.41906549 1.41906549 1.41906548 0.08 1.49184698 1.49184698 1.4918469 0.09 1.56831185 1.56831186 1.56831185 0.1 1.64811 1.64811 1.64810 Table 4.: Comparison of error for example 4.1 X 3SHBM 4SHBM [] 0.01 0.0 0.0 6. E -10 0.0 0.0 0.0 1.1 E -09 DOI: 10.990/58-11140914 www.iosrjournals.org 13 Page
0.03 0.0 0.0 1.3 E -09 0.04 0.0 0.0 1.6 E -09 0.05 0.0 0.0 1.3 E -09 0.06 0.0 1.0 E -09 4. E -10 0.0 0.0 1.0 E -09.4 E -09 0.08 0.0 1.0 E -09.4 E -09 0.09 1.0 E -09 0.0.5 E -09 0.1 0.0 1.0 E -09.3 E -09 Table 4.3: Numerical results for example 4. X EXACT 3SHBN 4SHBM 0.0 1.000000000 1.000000000 1.000000000 0.1 0.90483418 0.90483418 0.90483418 0. 0.8183053 0.8183053 0.8183053 0.3 0.408180 0.408180 0.408180 0.4 0.6030046 0.6030046 0.6030046 0.5 0.606530659 0.606530660 0.606530660 0.6 0.548811636 0.548811636 0.548811636 0. 0.49655303 0.496585304 0.496583304 0.8 0.44938964 0.44938964 0.44938965 0.9 0.406569590 0.406569659 0.406569660 1.0 0.3689441 0.3689441 0.3659441 Table 4.4: Comparison of Errors for example 4. X 3SHBM 4SHBM 0.0 0.0 0.0 0.1 0.0 0.0 0. 0.0 0.0 0.3 1.0 E -09 1.0 E -09 0.4 0.0 1.0 E -09 0.5 1. E -03 1. E -03 0.6 0.0 1.0 E -09 0. 1.0 E -09 1.0 E -09 0.8 0.0 1.0 E-09 0.9 0.0.1 E -08 0.0 0.0 1.0 E -09 V. Conclusion The desirable property of a numerical solution is to behave like the theoretical solution of the problem which can be seen in the above result. The implementation of the scheme is done with the aid of maple software. The method are tested and found to be consistent, zero stable and convergent. We implement the methods on two numerical examples and the numerical evidences shows that the methods are accurate and effective and therefore favourable. References [1]. Adeniyi, R.B., Alabi, M.O. and Folaranmi, R.O. (008): A Chebyshev collocation approach for a continuous formulation of hybrid methods for Initial Value Problems in ordinary differential equations. Journal of the Nigerian Association of Mathematical Physics, 1: 369-38 []. Awoyemi, D.O. (199). On some continuous linear multistep methods for initial value problem Unpublished Phd thesis, University of Ilorin, Ilorin Nigeria. [3]. Areo, E.A., Ademiluyi, R.A. and Babatola, P.O. (008). Accurate collocation multistep method for integration of first order ordinary differential equation, Int.J.Comp.maths,(1):15-. [4]. Ehigie, J.O., Okunuga, S.A. and Sofoluwa, A.B. (011). A class of -step continuous hybrid implicit linear multistep method for IVP's, Journal of Nigerian Mathematical Society, 30 145-161. [5]. Odejide, S.A and Adeniran, A.O. (01). A Hybrid linear collocation multistep scheme for solving first order initial value problems. Journal of Nigerian mathematical society. 31:9-41. [6]. Bolaji, B. (014). A Three-Step First derivative numerical integration scheme for IVP in ODEs. Journal of Mathematical Association Of Nigeria (ABACUS), 41(A):89-98. []. Abdullahi, A., Chollom, J. P and Yahaya, Y. A. (014). A family of Two-Step block generalized Adams Methods for the solution of Non-Stiff IVP in ODEs. Journal of Mathematical Association Of Nigeria (ABACUS), 41(A):6-. [8]. Henrici, P. (196). Discrete Variable Methods for ODEs, John Wiley and Sons, New York, USA. DOI: 10.990/58-11140914 www.iosrjournals.org 14 Page