9 J. Mah. Fund. Sci., Vol. 8, No.,, 9-5 New Seven-Sep Numerical Mehod for Direc Soluion of Fourh Order Ordinary Differenial Equaions Zurni Omar & John Olusola Kuboye Deparmen of Mahemaics, School of Quaniaive Sciences, Universii Uara Malaysia,, Sinok, Kedah, Malaysia E-mail: kubbysholly7@yahoo.com Absrac: A new numerical mehod for solving fourh order ordinary differenial equaions direcly is proposed in his paper. Inerpolaion and collocaion were employed in developing his mehod using seven seps. The use of he approimaed power series as an inerpolaion equaion was adoped in deriving he mehod. The basic properies of he new mehod such as zero-sabiliy, consisency, convergence and order are esablished. The numerical resuls indicae ha he new mehod gives beer accuracy han he eising mehods when i is applied o fourh ordinary differenial equaions. Keywords: block mehod; collocaion; direc soluion; inerpolaion; ordinary differenial equaions; seven-sep. Inroducion The general fourh order ordinary differenial equaions ODEs of he form,,,,, 3, [, ] s y f y y y y y ys s b are considered in his paper. Eq. can be solved by reducing i o is equivalen firs order sysem as menioned in [-9]. However, his approach suffers some sebacks, such as evaluaion of oo many funcions and heavy compuaion see [-]. Direc mehods of solving Eq. have been eamined by several researchers [,-]. They developed linear mulisep mehods using inerpolaion and collocaion whereby he use of he approimaed power series as a basis funcion was considered. Kayode [] developed an efficien zero-sable numerical mehod wih sep number k = and 5 for fourh order iniial value problems IVPs, which was implemened in predicor-correcor mode. Furhermore, a five-sep block mehod for solving fourh order ODEs direcly is presened in []. In addiion, in [3] and [] si-sep block mehods are developed for solving fourh order ODEs using a mulisep collocaion approach whereby collocaion poins are seleced a some grid poins. These mehods, however, have low accuracy. Received June h, 5, Revised December 5 h, 5, Acceped for publicaion February 8 h,. Copyrigh Published by ITB Journal Publisher, ISSN: 337-57, DOI:.5/j.mah.fund.sci..8..
Numerical Mehod for Direc Soluion of Fourh Order ODEs 95 In order o improve he accuracy of he eising mehods, his aricle proposes a new block mehod for direcly solving general fourh order IVPs of ODEs by increasing sep number k. Mehodology Le he approimae soluion o Eq. be a power series of he form k j y a j j where k is he sep number. The fourh derivaive of Eq. gives k j j j y j j j j3 a f, y, y, y, y 3 Eq. is inerpolaed a, i k5 k and Eq. 3 is collocaed a ni, i k. This gives a sysem of equaions in he form where X ni AX B a y nk5 a y nk a y nk3, B a3 ynk a y k nk, 3 k nk 5 nk5 nk 5 nk 5 nk 5... nk 5 3 k nk nk nk nk nk... nk 3 k nk3 nk 3 nk 3 nk3 nk3... nk3 3 k nk nk nk nk nk... nk A n... k k3 k k k n... k k3 k k n k nk... k k3 k k k n nk
9 Zurni Omar & John Olusola Kuboye By using Gaussian eliminaion, he values of, j as j k in Eq. are obained and hen subsiued ino Eq. o give a coninuous linear mulisep mehod in he form 5 k k j n j j n j j k j y y h f 5 For k = 7, we have 3 5 3 3 3 3 9 7 7, 9 8 7 5 3 7 5 3 G where 3 798 37 38 98 55 3 8 33 3533 9958 83 7 335 35 95 88 73 378 559 83 33 587 5 7 7 353 79 w w w w w w w w w w w w w w w w w w w w u u u u u u u u u u u v G 8 533 77 33 38 98 89 55 5 5873 33 73797 9 989 58 778 37 9 59 3575 3 85 89 33 989 3 5 v v v v v v v v v v v v v v v v v v v v v u u u u u u u 75 9 7 39 959 787 73533 989 98 8 7 55 87 39 338 3878 85 5 593 395 375 99 3 u u u u w w w w w w w w w w w w w w w w w w w w w w w.
Numerical Mehod for Direc Soluion of Fourh Order ODEs 97 The values of w, u and v are w = 975, u = 335, v = 3958, for n. h Eq. 5 is evaluaed a he non-inerpolaing poins, i.e. = -, -5, and, o produce he discree schemes. The firs, second and hird derivaives of Eq. 5 are evaluaed a all he poins wihin he inerval, i.e. = -, -5, -, -3, -, -, and, o give he derivaives of he discree schemes. These schemes are combined in a mari, whereby boh he y and he f funcion are muliplied by he inverse of he coefficiens of yn j, j k. This yields a block of he form where AY AY hay h BY hby h E F EF 3 N N N N N N N yn ynk y nk y nk yn yn k y n k y nk YN, YN Y N, Y N, ynk yn y n y n y nk fn fnk y nk fn fnk Y N, FN, FN. y n fnk fn If k = 7, we obain A, A,
98 Zurni Omar & John Olusola Kuboye 3 A, 5 7 3 9 B 3, 3 5 3 33 9 B 8, 5 8 9 E 37 975 35 7775 577 98 38, 7775 585 79 5 55 8993 77 37 55 55 39785 87898 5735 95 975 975 975 975 975 975 975 35 38 3335 335 9899 37 35 7775 7775 7775 7775 7775 7775 7775 3 9537 38 98535 5588 797 555 98 98 98 98 98 98 98 99 3 838 5998 8339 779 97 E 7775 7775 7775 7775 7775 7775 7775 9895 7335 775 995 388 5875 795 79 79 79 79 79 79 79 8 53 379 99 358 797 3 55 55 55 55 55 55 55 93958 37985 8878 3855 8793 353 83 77 77 77 77 77 77 77
Numerical Mehod for Direc Soluion of Fourh Order ODEs 99 The corresponding derivaives of Eq. are given by where y n y n 9 y 3 y n5 5 5 y n y 8 n7 7 9 n3 3 y n y n hy 8 h y h T f f f f f fn fn fn n n n n3 n 5 7 335799 58 757 93 575 39 583 8 388 38 38 38 38 38 38 38 337 387 383 373 53 93 377 398 835 835 835 835 835 835 835 835 537 835 538 355 98955 88 879 5 8 8 8 8 8 8 8 8 997 83 8 888 585 77 73 98 T 75 75 75 75 75 75 75 75 5875 55 5775 375 89875 575 5375 55 55 55 55 55 55 55 55 55 9 3 5 8 8 78 7 359 597738 39787 337 95 33 5573 5 58 58 58 58 58 58 58 58 y n n fn y n f n y n 3 3 f n3 yn yn hyn h M f n y n5 5 fn5 y n f n y n7 7 f n7 f
Zurni Omar & John Olusola Kuboye for 3855 79 75977 7858 598 8387 73 89 378 378 378 378 378 378 378 378 939 55 398 3995 95 39 383 835 835 835 835 835 835 835 835 9773 3 5997 398 98535 5 7899 555 8 8 8 8 8 8 8 8 58 75 9 57 99 73 7 9 M 75 75 75 75 75 75 75 75 5 75 5 5 35 9 875 35 757 757 757 757 757 757 757 757 597 8 8 7 95 98 8 35 35 35 35 35 35 35 35 85 3 3 8 553 38 38 38 38 38 38 38 38 fn y n fn y n f n y n 3 fn3 yn yn hn f n y n5 fn5 y n f n y n7 f n7 for 3799 3989 797 333 8857 99 35 375 9 9 9 9 9 9 9 9 5535 93 395 935 8 35 89 89 89 89 89 89 89 89 5 33975 885 935 55 885 85 5 78 8 78 8 N 95 95 95 95 95 95 95 95 755 375 975 5 35 755 75 75 9 9 9 9 9 9 9 9 7 7 7 557 539 9 93 93 9 539 557 78 78 78 78 78 78 78 78
Numerical Mehod for Direc Soluion of Fourh Order ODEs 3 Analysis of he Properies of he Block Mehod 3. Order of he Mehod The linear operaor associaed wih Eq. can be defined as L{ y : h} AY AY hay h BY h BY h E FN E FN 3 N N N N N Eq. 7 is epanded in Taylor series, which gives Ly [ : h] Cy Chy Ch y C h y p p p p p p The Eq. and he associaed linear operaor are said o have order p if C C C... C C C C, C. Therefore, our mehod s Eq. p p p p3 p has order p [8,8,8,8,8,8,8] T wih error consans 7 73 53 5 33 37 T p p C p [,,,,,, ].. Cp h y n 959 77 75 78 39 33 79 as he principal local runcaion error a poin n. is known 3. Zero Sabiliy of he Mehod Block Eq. is said o be zero-sable if he roos z,,..., N s of he firs characerisic polynomial z de za A saisfies z and he roo z has mulipliciy no greaer han he order of he differenial equaion, which is. Now, z de za A implies ha z z z. Hence z,,,,,,. Therefore, he new mehod s Eq. is convergen because i is zerosable and has order greaer han one. Numerical Eperimen The accuracy of he new mehod is eamined by solving he following differenial problems. iv Problem : y, y, y, y y, h. Eac soluion: 5 y
Zurni Omar & John Olusola Kuboye Problem : iv. y y,, y, y, 7 5. y, y, h. cos.sin Eac soluion: y iv 3 Problem 3: y 9 3 e, y y, y, y,. Eac Soluion: y e iv Problem : y y y y y y, Eac Soluion: y e The resuls generaed afer solving he above problems are shown in Tables. The following noaions are used in Tables 3 and : SPEB: Sequenial implemenaion of he -Poin Eplici Block Mehod. PPEB: Parallel implemenaion of he -Poin Eplici Block Mehod. S3PEB: Sequenial implemenaion of he 3-Poin Eplici Block Mehod. P3PEB: Parallel implemenaion of he 3-Poin Eplici Block Mehod. Table Comparison of new mehod wih [3] and [] for solving problem. Eac Soluion Compued Soluion Error in [3], Error in [], Error in New k = k = Mehod, k = 7..833333333.8333335 7.E-.7E-.87E-..9.9 8.99999-3.3333335E-.E+.3.35.35.599993E-9 5.9999999E-.E+..8533333333335.8533333333335 5.33E-9 7.75E-.E+.5.55.555 7.799979E-9 9.333333E-.87E-..87.7999997.89E-8.9E-9.75597E-.7.7583333333.75833983.3E-8.7E-9 3.573E-.8.8737.873359.E-8.533333E-9 3.573E-.9.99755.997999585.88E-8.99999E-9.7559E-..83333333333333.8333333385733.8335E-8.87E-9.75997E-
Numerical Mehod for Direc Soluion of Fourh Order ODEs 3 Table Comparison of new mehod wih [] and [] for solving problem. Eac Soluion Compued Soluion Error in Error in [], Error in New [], k = 5 k = 5 Mehod k = 7..8995837.8995837.55E-9.8355E-7.7859E-..5739533.5739533.3E-8.3933E- 5.E-.3.3859579797.3859579797.77E-8.893E-.755E-9..538839795.538839795.737E-7.9E-5.8E-9.5.79878537.79878537.338E-7.73E-5.3389E-9..79878537.79878537 9.59E-7 9.87E-5 8.737E-9.7.8939598.89395987.87E-.9E-.8E-8.8.975877.975879 3.57E-.57E-.7373E-8.9.389778533.389778535 5.88E- 3.88E-.889E-8..3785.3785 8.9E- 5.5E- 3.97E-8 Table 3 Comparison of new mehod wih [7] for solving problem 3. h-values New Mehod - 7-Sep Mehod -3 7-Sep Mehod - 7-Sep Mehod -5 7-Sep Mehod h-values Omar In [7] Number Error in New Mehod, Error in [7] of Seps k = 7 k = 8 SPEB 5 3.579E-.778E- PPEB 5 3.579E-.778E- S3PEB 39.5555E-.778E- P3PEB 39.5555E-.778E- SPEB 5.E-.778E-3 PPEB 5.E-.778E- S3PEB 339.7E-.778E-3 P3PEB 339.7E-.778E- SPEB 5.3973E-9.8E- PPEB 5.3973E-9.8E- S3PEB 3339 3.39E-.8E- P3PEB 3339 3.39E-.8E- SPEB 5.95E-9.E-5 PPEB 5.95E-9.E-5 S3PEB 33339 9.79E-.E-5 P3PEB 33339 9.79E-.E-5 Table Comparison of new mehod wih [7] for solving problem. New Mehod - 7-Sep Mehod -3 7-Sep Mehod - 7-Sep Mehod -5 7-Sep Mehod Omar In [7] Number Error in New Mehod, Error in [7] of Seps k = 7 k = 8 SPEB 5 3.75589E- 8.37E- PPEB 5 3.75589E- 8.37E- S3PEB 39.7E- 8.375E- P3PEB 39.7E- 8.375E- SPEB 5 3.7587E-3 8.3E-5 PPEB 5 3.7587E-3 8.3E-5 S3PEB 339.755E- 8.3E-5 P3PEB 339.755E- 8.3E-5 SPEB 5 7.977E- 8.3353E- PPEB 5 7.977E- 8.3353E- S3PEB 3339.3555E-3 8.3353E- P3PEB 3339.3555E-3 8.3353E- SPEB 5.8E- 8.33E-7 PPEB 5.8E- 8.33E-7 S3PEB 33339.977E- 8.333E-7 P3PEB 33339.977E- 8.333E-7
Zurni Omar & John Olusola Kuboye 5 Conclusion A seven-sep block mehod for he soluion of fourh order ODEs is proposed in his paper. The new mehod was used o solve fourh order IVPs of ODEs. The numerical resuls were compared wih he eising mehods. The new mehod performed beer han he mehods in [,-], which employed 5 and seps refer o Tables and. This implies ha beer accuracy can be achieved when sep number k is increased. The accuracy of he new mehod was also found o be beer han he mehod in [7], which was developed hrough numerical inegraion using 8 seps refer o Tables 3 and. Thus, based on he numerical resuls, he new mehod ouperformed he eising mehods in erms of accuracy and should be considered as a viable alernaive for direcly solving fourh order iniial value problems. References [] Awari, Y.S. & Abada, A.A., A Class of Seven Poin Zero Sable Coninuous Block Mehod for Soluion of Second Order Ordinary Differenial Equaion, Inernaional Journal and Mahemaics and Saisics Invenion IJMSI,, pp. 7-5,. [] Anake, T.A., Awoyemi D.O., & Adesanya, A.O., One-Sep Implici Hybrid Block Mehod for he Direc Soluion of General Second Order Ordinary Differenial Equaions, IEANG Inernaional Journal of Applied Mahemaics,, IJAM,. [3] Bun, R.A. & Vasil yer., A Numerical Mehod for Solving Differenial Equaion of Any Orders, Comp. Mah. Phys., 33, pp. 37-33, 99. [] Brugnano, L. & Trigiane, D., Solving Differenial Problems by Mulisep IVPs and BVP Mehods, Gordon and Breach Science Publishers, 998. [5] Lamber, J.D., Compuaional Mehods in Ordinary Differenial Equaions, John Wiley & Sons Inc. 973. [] Senu, N., Suleiman, M., Ismail, F., & Ohman, M., A Singly Diagonally Implici Runge-kua-Nysrom Mehod for Solving Oscillaory Problems, IEANG Inernaional Journal of Applied Mahemaics,, IJAM,. [7] Omar, Z., Parallel Block Mehods for Solving Higher Order Ordinary Differenial Equaion Direcly. Docor Thesis in he Faculy of Science and Environmenal Sudies Universii Pura Malaysia, Unpublished, 999. [8] Omar, Z. & Kuboye, J.O., Compuaion of an Accurae Implici Block Mehod for Solving Third Order Ordinary Differenial Equaions Direcly, Global Journal of Pure and Applied Mahemaics,, pp. 77-8, 5.
Numerical Mehod for Direc Soluion of Fourh Order ODEs 5 [9] Omar, Z. & Kuboye, J.O., Developing Block Mehod of Order Seven for Solving Third Order Ordinary Differenial Equaions Direcly using Mulisep Collocaion Approach. Inernaional Journal of Applied Mahemaics and Saisics, 533, pp. 5-73, 5. [] Adesanya, A.O., Momoh, A.A., Alkali, M.A., & Tahir, A., A Five Seps Block Mehod for he Soluion of fourh Order Ordinary Differenial Equaions, Inernaional Journal of Engineering Research and Applicaions IJERA,, pp. 99-998,. [] Bucher, J.C, Numerical Mehods for Ordinary Differenial Equaions, John Wiley & Sons, New York, 3. [] Kayode, S.J., An Efficien Zero Sable Mehod for Fourh Order Ordinary Differenial Equaions, In. J. Mah. Sci, pp. -, 8. d.doi.org/.55/8/3 [3] Mohammed, U., A Si Sep Block Mehod for Soluion of Fourh Order Ordinary Differenial Equaions, The Pacific Journal of Science and Technology,, pp. 59-5,. [] Olabode, B.T., A Si Sep Scheme for he Soluion of Fourh Order Ordinary Differenial Equaions, The Pacific Journal of Science and Technology,, pp. 3-8, 9.