BLOCK PROCEDURE WITH IMPLICIT SIXTH ORDER LINEAR MULTISTEP METHOD USING LEGENDRE POLYNOMIALS FOR SOLVING STIFF INITIAL VALUE PROBLEMS
DOI:
https://doi.org/10.4314/jfas.v11i1.1Keywords:
Collocation, Interpolation, Legendre Polynomials, Linear Multistep Method, StiffAbstract
In this paper, a discrete implicit linear multistep method in block form of uniform step size for the solution of first-order ordinary differential equations is presented using the power series as a basis function. To improve the accuracy of the method, a perturbation term is added to the approximated solution. The method is based on collocation of the differential equation and interpolation of the approximate solution using power series at the grid points. The procedure yields four linear multistep schemes which are combined as simultaneous numerical integrators to form block method. The method is found to be consistent and zero-stable, and hence convergent. The accuracy of the method is tested with some standard stiff first order initial value problems. The results are compared with the two point implicit block backward difference formula and two point implicit block extended backward difference formula. All numerical examples show that our proposed method has a better accuracy than the existing literature.
Downloads
References
[2] Suleiman, M.B., Musa, M., Ismail, F., An Implicit 2-point Block Extended Backward Differentiation Formula for Integration of Stiff Initial Value Problems. Malaysian Journal of Mathematical Science, 2015, 9(1):33-51
[3] Alt, R., A-stable one-step methods with step-size control for stiff systems of ordinary differential equations. J Comput Appl Math, 1978, 4, 29-35
[4] Cash. J., On the integration of stiff systems of ODEs using extended backward differentiation formulae. Numer Math, 1980, 34,235-246
[5] Alvarez, J., Rojo, J., An improved class of generalized Runge-Kutta methods for stiff problems. Part I: The scalar case. Appl Math Comput, 2002, 130, 537-560
[6] Ibrahim, Z., Othman, K., Suleiman, M., Fixed coefficients block backward differentiation formulas for the numerical solution of stiff ordinary differential equations. Eur J Sci Res, 2008, 21, 508-520
[7] Yatim SAM, Ibrahim ZB, Othman KI, Suleiman MB, A quantitative comparison of numerical method for solving stiff ordinary differential equations. Math Probl Eng, unkown, 2011, 193691
[8] Zawawi, I.S.., Ibrahim, Z.B., Othman, K.I., Derivation of diagonally implicit block backward differentiation formulas for solving stiff IVPs. Math Probl Eng,unknown 179231.
[9] Abasi, N., Suleiman M., Abbasi, N., Musa, H. , 2-point block BDF method with off-step points for solving stiff ODEs. Journal of Soft Computing and Applications, 2014, 1-15,.
[10] Hirayama, H., Arbitrary order and A-stable numerical method for solving algebraic ordinary differential equation by power series. In 2nd International Conference on Mathematics and Computers in Physics, 2000
[11] Abualnaja, K., A block procedure with LMMs using Legendre polynomials for Solving ODEs. Applied mathematics series, 2015, 6(1), 717-723
[12] Lambert, J., Computational methods in ordinary differential equations. New York: Willey and Sons, 1973.
[13] Henrici, P., Discrete variable methods for ODEs. New York: John Wiley, 1962.
[14] Fatunla, S.O., A class of block methods for second order IVPs. Int. J. Comput. Math., 1994, 55, 119-133
[15] Suli,E. and Mayers, D.F., An Introduction To Numerical Analysis (1st ed.). New York: Cambridge University Press, 2003.
[16] Shampine, L., Watts, H., Block implicit one-step methods. Journal of Computer Maths., 1969, 23, 731-740
[17] Fatunla, S., Numerical Methods for initial value problems for ordinary differential equations. Boston: U.S.A Academy Press, 1988.
[18] Randall, J. L., (2004). Finite Difference Methods for differential equations. University of Washington.
[19] Naghmeh Abasi, Mohamed Suleiman, Neda Abbasi, & Hamisu Musa, 2-Point Block BDF Method with Off-Step Points for Solving Stiff ODEs. Journal of Soft Computing and Applications, 2014, 1-15, 2014
[20] Ibrahim, Z.B., Block multistep methods for solving ordinary differential equations, PhD Thesis, Universiti Putra Malysia, (2006).
Submissions
