# Exponentially Fitted and Trigonometrically Fitted Explicit Modified Runge-Kutta Type Methods for Solvingy′′′x=fx,y,y′

Authors: N. Ghawadri; N. Senu; F. Ismail; Z. B. Ibrahim
Journal: Journal of Applied Mathematics
 (2018)

---

## Abstract

Exponentially fitted and trigonometrically fitted explicit modified Runge-Kutta type (MRKT) methods for solvingy′′′x=fx,y,y′ are derived in this paper. These methods are constructed which exactly integrate initial value problems whose solutions are linear combinations of the set functions eωx and e-ωx for exponentially fitted and sin⁡ωx and cos⁡ωx for trigonometrically fitted with ω∈R being the principal frequency of the problem and the frequency will be used to raise the accuracy of the methods. The new four-stage fifth-order exponentially fitted and trigonometrically fitted explicit MRKT methods are called EFMRKT5 and TFMRKT5, respectively, for solving initial value problems whose solutions involve exponential or trigonometric functions. The numerical results indicate that the new exponentially fitted and trigonometrically fitted explicit modified Runge-Kutta type methods are more efficient than existing methods in the literature.


---

## Body

## 1. Introduction

This work deals with exponentially fitted and trigonometrically fitted modified Runge-Kutta type methods for solving third-order ordinary differential equations (ODEs)(1)y′′′x=fx,yx,y′x,yx0=y0,y′x0=y0′,y′′x0=y0″,x≥x0.This sort of problems is often found in numerous physical problems like thin film flow, gravity-driven flows, electromagnetic waves, and so on. In the past and recent years many researchers constructed exponentially fitted and trigonometrically fitted explicit Runge-Kutta methods for solving first-order and second-order ordinary differential equations. Paternoster [1] developed Runge-Kutta-Nyström methods for ODEs with periodic solutions based on trigonometric polynomials. Vanden Berghe et al. [2] developed exponentially fitted Runge-Kutta methods. Simos [3] extended exponentially fitted Runge-Kutta methods for the numerical solution of the Schrodinger equation and related problems. Kalogiratou et al. [[4, 5]] constructed trigonometrically and exponentially fitted Runge-Kutta-Nyström methods for the numerical solution of the Schrodinger equation and related problems which is eighth algebraic order. Next Simos et al. [6] constructed exponentially fitted Runge-Kutta-Nyström method for the numerical solution of initial value problems with oscillating solutions. Sakas et al. [7] developed a fifth algebraic order trigonometrically fitted modified Runge-Kutta Zonneveld method for the numerical solution of orbital problems. Van de Vyver [8] in 2005 constructed Runge-Kutta-Nyström pair for the numerical integration of perturbed oscillators. Then Yang et al. [9] constructed trigonometrically fitted adapted Runge-Kutta-Nyström methods for perturbed oscillators. Recently, Demba et al. [10] constructed an explicit trigonometrically fitted Runge-Kutta-Nyström method using Simos technique.In this paper we construct explicit exponentially fitted and trigonometrically fitted modified Runge-Kutta type methods with four-stage fifth-order, called EFMRKT5 and TFMRKT5, respectively. Section2 discussed the oscillatory and nonoscillatory properties of the third-order linear differential equation. In Section 3, the necessary conditions and the derivation for exponentially fitted and trigonometrically fitted modified Runge-Kutta type methods for solving third-order ODEs are given. The error analysis of the new EFMRKT5 and TFMRKT5 methods was discussed in Section 4, respectively. The effectiveness of the new methods when compared with existing methods is given in Section 5. The thin film flow problem is discussed in Section 6.


## 2. Third-Order Linear Differential Equation with Oscillating and Nonoscillating Solutions

This section discusses the oscillatory and nonoscillatory properties of the third-order linear differential equation(2)y′′′x+pxy′+qxy=0.A solution of (2) will be said to be oscillatory if it changes signs for arbitrarily large values of x. The other solutions will be said to be nonoscillatory.Ifpx<0 and qx<0 are constants, then it is easy to show that if (2) has an oscillatory solution, then there are two linearly independent oscillatory solutions of (2) whose zeroes separate and such that any oscillatory solution of (2) is a linear combination of them. Assuming that px,p′x, and qx are continuous on 0,+∞ the following will be established (see [11–14]).Definition 1.
A solution of (2) will be called oscillatory iff it has an infinity of zeroes in(0,+∞) and nonoscillatory iff it has but a finite number of zeroes in this interval. Equation (2) is said to be oscillatory iff it has at least one (nontrivial) oscillatory solution and nonoscillatory iff all of its (nontrivial) solutions are nonoscillatory.Particularly, this paper deals with two cases based on (2) when qx=0, as follows:(i)
y′′′(x)=py′,(p>0); it is clear that the characteristic roots equations are real and one of them is zero; then solutions will consist of exponential functions.(ii)
y′′′(x)=-py′,(p>0); one of the characteristic roots equations is zero and another two are conjugate roots and the solutions are in oscillatory form,wherep is constant.


## 3. Exponentially Fitted and Trigonometrically Fitted MRKT Methods

In this section, we will determine the conditions and develop exponentially fitted and trigonometrically fitted MRKT methods. In order to construct the exponentially fitted and trigonometrically fitted MRKT methods, the extraγi and γ^i are absolutely necessary to insert at each stage and the MRKT methods is given as follows:(3)yn+1=yn+hyn′+h22yn′′+h3∑i=1sbiki,(4)yn+1′=yn′+hyn′′+h2∑i=1sbi′ki,(5)yn+1′′=yn′′+h∑i=1sbi′′ki,where(6)k1=fxn,yn,yn′,(7)ki=fxn+cih,γiyn+hciyn′+h22ci2yn′′+h3∑j=1saijkj,yn′+γ^ihciyn′′+h2∑j=1sa^ijkjfori=2,3,…,s.The parameters of the MRKT methods areci,aij,a^ij,bi,bi′,bi′′,γi and γ^i for i=1,2,…,s and j=1,2,…,s are assumed to be real. If aij=0 and a^ij=0 for i⩽j, it is an explicit method and otherwise implicit method.The MRKT method can be expressed in Butcher notation using the table of coefficients as follows (see Table1).Table 1
The Butcher tableau MRKT method.
c γ γ ^ A A ^ b b ′ b ′ ′


### 3.1. Exponentially Fitted MRKT Method

To construct the exponentially fitted Runge-Kutta type four-stage fifth-order method the functionseωx and e-ωx need to integrate exactly at each stage; therefore the following four equations are obtained:(8)e±civ=γi±civ+12ci2v2±v3∑j=1saije±cjv,(9)e±civ=1±γ^iciv±v2∑j=1sa^ije±cjv,and six more equations corresponding toy, y′, and y′′:(10)e±v=1±v+12v2±v3∑i=1sbie±civ,(11)e±v=1±v+v2∑i=1sbi′e±civ,(12)e±v=1±v∑i=1sbi′′e±civ,where v=ωh, ω∈R. The relations cosh⁡v=ev+e-v/2 and sinh⁡v=ev-e-v/2 will be used in the derivation process. The following order conditions are obtained: (13)cosh⁡vci=γi+12v2ci2+v3∑j=1i-1aijsinh⁡vcj,(14)sinh⁡vci=vci+v3∑j=1i-1aijcosh⁡vcj,(15)cosh⁡vci=1+v2∑j=1i-1a^ijcosh⁡vcj,(16)sinh⁡vci=γ^iciv+v2∑j=1i-1a^ijsinh⁡vcj,and six equations corresponding to y, y′, and y′′:(17)cosh⁡v=1+12v2+v3∑i=1sbisinh⁡vci,(18)sinh⁡v=v+v3∑i=1sbicosh⁡vci,(19)cosh⁡v=1+v2∑i=1sbi′cosh⁡vci,(20)sinh⁡v=v+v2∑i=1sbi′sinh⁡vci,(21)cosh⁡v=1+v∑i=1sbi′′sinh⁡vci,(22)sinh⁡v=v∑i=1sbi′′cosh⁡vci.Solving (13) to (16), we find ai,i-1, a^i,i-1,γi, and γ^i.(23)γi=cosh⁡vci-12v2ci2-v3∑j=1i-1ai,jsinh⁡vcj,(24)ai,i-1=sinh⁡vci-vci-v3∑j=1i-2ai,jcosh⁡vcjv3cosh⁡vci-1,(25)a^i,i-1=cosh⁡vci-1-v2∑j=1i-2a^i,jcosh⁡vcjv2cosh⁡vci-1,(26)γ^i=sinh⁡vci-v2∑j=1i-2a^i,jsinh⁡vcjvci,i=2,…,s.Referring to the following fifth-order four-stage method developed by Fawzi et al. [15]:(27)c1=0,c2=15,c3=23,c4=1,a21=0,a31=-494860,a41=750,a42=-150,a^31=-127,a^41=310,a^42=-235,b3=3112,b4=0,b3′=956,b4′=0,b3′′=756,b4′=548,we solve (23) to (26) and let a^21, a^32, a^43, a32, a43, γ2, γ3, γ4, γ^2, γ^3, and γ^4 be free parameters and yields.(28)a^21=cosh⁡v/5-1v2,a^32=cosh⁡2v/3-1+1/27v2v2cosh⁡v/5,a^43=cosh⁡v-1-v23/10-2/35cosh⁡v/5v2cosh⁡2v/3,a32=sinh⁡2v/3-2v/3+49/4860v3v3cosh⁡c2v,a43=sinh⁡v-v-7/50v3+1/50v3cosh⁡v/5v3cosh⁡2v/3,γ2=cosh⁡v5-v250,γ3=cosh⁡2v3-2v29-301v34860sinh⁡v5,γ4=cosh⁡v-v22-v3-150sinh⁡v5+125sinh⁡2v3,γ^2=5vsinh⁡v5,γ^3=23vsinh⁡2v3-2cosh⁡2v/3-1+v2/27sinh⁡v/53cosh⁡v/5v2,γ^4=1vsinh⁡v+2v35sinh⁡v5-cosh⁡v-1-3v2/10-2v2/35cosh⁡v/5sinh⁡2v/3vcosh⁡2v/3.Next, we solve (17) to (22) and use the above coefficients to find b1, b2, b1′, b2′, b1′′, and b2′′.(29)b1=3112cosh⁡v/5sinh⁡2v/3-cosh⁡2v/3sinh⁡v/5sinh⁡v/5+-2cosh⁡v/5cosh⁡v+2cosh⁡v/5+cosh⁡v/5v2+2sinh⁡vsinh⁡v/5-2vsinh⁡v/52v3sinh⁡v/5,b2=-3112sinh⁡2v/3sinh⁡v/5--2cosh⁡v+2+v22v3sinh⁡v/5,b1′=956cosh⁡v/5sinh⁡2v/3-cosh⁡2v/3sinh⁡v/5sinh⁡v/5-cosh⁡v/5sinh⁡v-cosh⁡v/5v-cosh⁡vsinh⁡v/5+sinh⁡v/5v2sinh⁡v/5,b2′=-956sinh⁡2v/3sinh⁡v/5+sinh⁡v-vv2sinh⁡v/5,b2′′=-1+27/56vsinh⁡2v/3+5/48vsinh⁡v-cosh⁡vvsinh⁡v/5,b1′′=27/56cosh⁡v/5vsinh⁡2v/3+5/48cosh⁡v/5vsinh⁡v-cosh⁡v/5cosh⁡v+cosh⁡v/5-27/56vcosh⁡2v/3sinh⁡v/5-5/48vcosh⁡vsinh⁡v/5+sinh⁡vsinh⁡v/5vsinh⁡v/5.These lead to our new exponentially fitted Runge-Kutta type four-stage fifth-order explicit MRKT method denoted as EFMRKT5. The corresponding Taylor series expansion of the solution is given by(30)b1=148+12160v2+101136080000v4-5713183708000000v6-1133033981841914000000000v8-57722879134057055132000000000v10+⋯,b2=542-12160v2-2921136080000v4-136145927000000v6+5429358740920957000000000v8+6964030429670285275660000000000v10+⋯,b1′=124-17283500v4-149218700000v6-1055069341007975000000v8-230256892538959377500000000v10+⋯,b2′=2584+2412268000v4+228716123600000v6+1127781372728063800000000v8+1359935168355857106305000000000v10+⋯,b1′′=124+121600v4+16758320000v6+52838916533720000000v8+9673435845851000000000v10+⋯b2′′=125336-121600v4-2867408240000v6-302210916533720000000v8-1457821730731375000000v10+⋯,a^21=150+115000v2+111250000v4+115750000000v6+135437500000000v8+1116943750000000000v10+⋯,a^32=727+3712150v2+2876561000v4-68827516678750000v6+70974171674039150000000v8-39720321233621487034437500000000v10+⋯,a^43=935-433000v2+33231350000v4-111136579255150000000v6+26973882539344452500000000v8-72136159838800151151196250000000000v10+⋯,a32=2894860-67729000v2+261412755620000v4-467154139858075000000v6+22160081031104865839000000000v8-15664491766661484759886861250000000000v10+⋯,a43=7150-221135000v2+114463637875000v4-101047288934445250000000v6+22257681969742625996875000000v8-8431888341871633389770349418750000000000v10+⋯,γ2=1+115000v4+111250000v6+115750000000v8+135437500000000v10+⋯,γ3=1-59225v2-95691215000v4-74831911250000v6+2610442911481750000000v8-49593722925833937500000000v10+⋯γ4=1-v2+191000v4+22674050000v6-4858267255150000000v8+65060629344452500000000v10+⋯,γ^2=1+1150v2+175000v4+178750000v6+1141750000000v8+1389812500000000v10+⋯,γ^3=1-1270v2+1360750v4-2032717609375v6+1361092906317968750v8-27268073596568486328125v10+⋯,γ^4=1+1150v2+7915000v4-21619512126250000v6+78384961425250000000v8-5664933289891170503987500000000v10+⋯where γ1=1,γ^1=1.This results in the new method called EFMRKT5. Asv→0, the coefficients b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of the new method EFMRKT5 reduce to the coefficients of the original method RKT5. That is to say, b1(0),b2(0),b1′(0),b2′(0),b2′′(0),b2′′(0),a32(0),a42(0),a^21(0), a^32(0),a^43(0),γ2(0),γ3(0),γ4(0),γ^2(0),γ^3(0), and γ^4(0) are identical to b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of RKT5 method. Other than that, v→0, as EFMRKT5 method will have the same error constant as RKT5 method.


### 3.2. Trigonometrically Fitted MRKT Method

Exponentially fitted method leads to trigonometrically fitted method when replacingv=wh with iv and solving (8) to (9) to find ai,i-1, a^i,i-1,γi, and γ^i.(31)γi=cos⁡vci-12v2ci2-v3∑j=1i-1ai,jsin⁡vcj,(32)a^i,i-1=1-cos⁡vci-v2∑j=1i-2a^i,jcos⁡vcjv2cos⁡vci-1,(33)ai,i-1=-sin⁡vci+v·ci-v3∑j=1i-2ai,jcos⁡vcjv3cos⁡vci-1,(34)γ^i=sin⁡vci+v2∑j=1i-2a^i,jsin⁡vcjvci,i=2,…,s.Consider the same coefficients of fifth-order four-stage method developed by Fawzi et al.[15] as in Section 3.1. Solving the (31) to (34) and letting a^21, a^32, a^43, a32, a43, γ2, γ3, γ4, γ^2, γ^3, and γ^4 be free parameters will give(35)a^21=1-cosv/5v2,a^32=1-cos2v/3+1/27v2v2cosv/5,a^43=1-cosv-v23/10-2/35cosv/5v2cos2v/3,a32=-sin2v/3+2v/3+49/4860v3v3cosv/5,a43=v-sinv-v37/50-1/50cosv/5v3cos2v/3,γ2=cosv5-v250,γ3=cos2v3-2v29-301v34860sinv5,γ4=cosv-v22-v3-150sinv5+125sin2v3,γ^2=5vsinv5,γ^3=23vsin2v3+21-cos2v/3+v2/27sinv/53v2cosv/35,γ^4=sin⁡vv-2v35sin⁡v5+1-cos⁡v-3v2/10-2v2/5cos⁡v/5sin⁡2v/3vcos⁡2v/3.Next, solving (10) to (12), and using the above Fawzi coefficients to find b1, b2,b1′, b2′, b1′′, and b2′′,(36)b1=3cos⁡v/5sin⁡2v/3-3cos⁡2v/3sin⁡v/5112sin⁡v/5-2cos⁡v/5cos⁡v-2cos⁡v/5+cos⁡v/5v2+2sin⁡vsin⁡v/5-2vsin⁡v/52v3sin⁡v/5,b2=-3112sin⁡2v/3sin⁡v/5+2cos⁡v-2+v22v3sin⁡v/5,b1′=9cos⁡v/5sin⁡2v/3-9cos⁡2v/3sin⁡v/556sin⁡v/5+cos⁡v/5sin⁡v-vcos⁡v/5-cos⁡vsin⁡v/5+sin⁡v/5v2sin⁡v/5,b1′′=-25v/243sin⁡v+80/81cos⁡v-80/81+23v/48sin⁡2v/380v/81sin⁡v/5-v/125sin⁡2v/3,b2′′=-1+27/56vsin⁡2v/3+5/48vsin⁡v+cos⁡vvsin⁡v/5,b2′=-27/56cos⁡v/5vsin⁡2v/3+5/48cos⁡v/5vsin⁡v+cos⁡v/5cos⁡v-cos⁡v/5-27/56vcos⁡2v/3sin⁡v/5-5/48vcos⁡vsin⁡v/5+sinvsin⁡v/5vsin⁡v/5These lead to our new explicit trigonometrically fitted MRKT which is called TFMRKT5 method. The corresponding Taylor series expansion of the solution is given by(37)b1=148-12160v2+101136080000v4+5713183708000000v6-1133033981841914000000000v8+57722879134057055132000000000v10+⋯,b2=542+12160v2-2921136080000v4+136145927000000v6+5429358740920957000000000v8-6964030429670285275660000000000v10+⋯,b1′=124-17283500v4+149218700000v6-1055069341007975000000v8+230256892538959377500000000v10+⋯,b2′=2584+2412268000v4-228716123600000v6+1127781372728063800000000v8-1359935168355857106305000000000v10+⋯,b1′′=124+121600v4-16758320000v6+52838916533720000000v8-9673435845851000000000v10+⋯,b2′′=125336-121600v4+2867408240000v6-302210916533720000000v8+1457821730731375000000v10+⋯,a^21=150-115000v2+111250000v4-115750000000v6+135437500000000v8-1116943750000000000v10+⋯,a^32=727-3712150v2+2876561000v4+68827516678750000v6+70974171674039150000000v8+39720321233621487034437500000000v10+⋯,a^43=935+433000v2+33231350000v4+111136579255150000000v6+26973882539344452500000000v8+72136159838800151151196250000000000v10+⋯,a32=2894860+67729000v2+261412755620000v4+467154139858075000000v6+22160081031104865839000000000v8+15664491766661484759886861250000000000v10+⋯,a43=7150+221135000v2+114463637875000v4+101047288934445250000000v6+22257681969742625996875000000v8+8431888341871633389770349418750000000000v10+⋯,γ2=1-125v2+115000v4-111250000v6+115750000000v8-135437500000000v10+⋯,γ3=1-59225v2-95691215000v4-74831911250000v6+2610442911481750000000v8-49593722925833937500000000v10+⋯,γ4=1-v2+191000v4+22674050000v6-4858267255150000000v8+65060629344452500000000v10+⋯,γ^2=1-1150v2+175000v4-178750000v6+1141750000000v8-1389812500000000v10+⋯,γ^3=1+1270v2+1360750v4+2032717609375v6+1361092906317968750v8+27268073596568486328125v10+⋯,γ^4=1-1150v2+7915000v4+21619512126250000v6+78384961425250000000v8+5664933289891170503987500000000v10+⋯,whereγ1=1,γ^1=1.This results in the new method called TFMRKT5. Asv→0, the coefficients b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of the new method TFMRKT5 reduce to the coefficients of the original method RKT5. That is to say, b1(0),b2(0),b1′(0),b2′(0),b2′′(0),b2′′(0),a32(0),a42(0), a^21(0),a^32(0),a^43(0),γ2(0),γ3(0),γ4(0),γ^2(0),γ^3(0), and γ^4(0) are identical to b1,b2,b1′,b2′,b2′′,b2′′,a32,a42, a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of RKT5 method. Other than that, v→0, as TFMRKT5 method will have the same error constant as RKT5 method.


## 3.1. Exponentially Fitted MRKT Method

To construct the exponentially fitted Runge-Kutta type four-stage fifth-order method the functionseωx and e-ωx need to integrate exactly at each stage; therefore the following four equations are obtained:(8)e±civ=γi±civ+12ci2v2±v3∑j=1saije±cjv,(9)e±civ=1±γ^iciv±v2∑j=1sa^ije±cjv,and six more equations corresponding toy, y′, and y′′:(10)e±v=1±v+12v2±v3∑i=1sbie±civ,(11)e±v=1±v+v2∑i=1sbi′e±civ,(12)e±v=1±v∑i=1sbi′′e±civ,where v=ωh, ω∈R. The relations cosh⁡v=ev+e-v/2 and sinh⁡v=ev-e-v/2 will be used in the derivation process. The following order conditions are obtained: (13)cosh⁡vci=γi+12v2ci2+v3∑j=1i-1aijsinh⁡vcj,(14)sinh⁡vci=vci+v3∑j=1i-1aijcosh⁡vcj,(15)cosh⁡vci=1+v2∑j=1i-1a^ijcosh⁡vcj,(16)sinh⁡vci=γ^iciv+v2∑j=1i-1a^ijsinh⁡vcj,and six equations corresponding to y, y′, and y′′:(17)cosh⁡v=1+12v2+v3∑i=1sbisinh⁡vci,(18)sinh⁡v=v+v3∑i=1sbicosh⁡vci,(19)cosh⁡v=1+v2∑i=1sbi′cosh⁡vci,(20)sinh⁡v=v+v2∑i=1sbi′sinh⁡vci,(21)cosh⁡v=1+v∑i=1sbi′′sinh⁡vci,(22)sinh⁡v=v∑i=1sbi′′cosh⁡vci.Solving (13) to (16), we find ai,i-1, a^i,i-1,γi, and γ^i.(23)γi=cosh⁡vci-12v2ci2-v3∑j=1i-1ai,jsinh⁡vcj,(24)ai,i-1=sinh⁡vci-vci-v3∑j=1i-2ai,jcosh⁡vcjv3cosh⁡vci-1,(25)a^i,i-1=cosh⁡vci-1-v2∑j=1i-2a^i,jcosh⁡vcjv2cosh⁡vci-1,(26)γ^i=sinh⁡vci-v2∑j=1i-2a^i,jsinh⁡vcjvci,i=2,…,s.Referring to the following fifth-order four-stage method developed by Fawzi et al. [15]:(27)c1=0,c2=15,c3=23,c4=1,a21=0,a31=-494860,a41=750,a42=-150,a^31=-127,a^41=310,a^42=-235,b3=3112,b4=0,b3′=956,b4′=0,b3′′=756,b4′=548,we solve (23) to (26) and let a^21, a^32, a^43, a32, a43, γ2, γ3, γ4, γ^2, γ^3, and γ^4 be free parameters and yields.(28)a^21=cosh⁡v/5-1v2,a^32=cosh⁡2v/3-1+1/27v2v2cosh⁡v/5,a^43=cosh⁡v-1-v23/10-2/35cosh⁡v/5v2cosh⁡2v/3,a32=sinh⁡2v/3-2v/3+49/4860v3v3cosh⁡c2v,a43=sinh⁡v-v-7/50v3+1/50v3cosh⁡v/5v3cosh⁡2v/3,γ2=cosh⁡v5-v250,γ3=cosh⁡2v3-2v29-301v34860sinh⁡v5,γ4=cosh⁡v-v22-v3-150sinh⁡v5+125sinh⁡2v3,γ^2=5vsinh⁡v5,γ^3=23vsinh⁡2v3-2cosh⁡2v/3-1+v2/27sinh⁡v/53cosh⁡v/5v2,γ^4=1vsinh⁡v+2v35sinh⁡v5-cosh⁡v-1-3v2/10-2v2/35cosh⁡v/5sinh⁡2v/3vcosh⁡2v/3.Next, we solve (17) to (22) and use the above coefficients to find b1, b2, b1′, b2′, b1′′, and b2′′.(29)b1=3112cosh⁡v/5sinh⁡2v/3-cosh⁡2v/3sinh⁡v/5sinh⁡v/5+-2cosh⁡v/5cosh⁡v+2cosh⁡v/5+cosh⁡v/5v2+2sinh⁡vsinh⁡v/5-2vsinh⁡v/52v3sinh⁡v/5,b2=-3112sinh⁡2v/3sinh⁡v/5--2cosh⁡v+2+v22v3sinh⁡v/5,b1′=956cosh⁡v/5sinh⁡2v/3-cosh⁡2v/3sinh⁡v/5sinh⁡v/5-cosh⁡v/5sinh⁡v-cosh⁡v/5v-cosh⁡vsinh⁡v/5+sinh⁡v/5v2sinh⁡v/5,b2′=-956sinh⁡2v/3sinh⁡v/5+sinh⁡v-vv2sinh⁡v/5,b2′′=-1+27/56vsinh⁡2v/3+5/48vsinh⁡v-cosh⁡vvsinh⁡v/5,b1′′=27/56cosh⁡v/5vsinh⁡2v/3+5/48cosh⁡v/5vsinh⁡v-cosh⁡v/5cosh⁡v+cosh⁡v/5-27/56vcosh⁡2v/3sinh⁡v/5-5/48vcosh⁡vsinh⁡v/5+sinh⁡vsinh⁡v/5vsinh⁡v/5.These lead to our new exponentially fitted Runge-Kutta type four-stage fifth-order explicit MRKT method denoted as EFMRKT5. The corresponding Taylor series expansion of the solution is given by(30)b1=148+12160v2+101136080000v4-5713183708000000v6-1133033981841914000000000v8-57722879134057055132000000000v10+⋯,b2=542-12160v2-2921136080000v4-136145927000000v6+5429358740920957000000000v8+6964030429670285275660000000000v10+⋯,b1′=124-17283500v4-149218700000v6-1055069341007975000000v8-230256892538959377500000000v10+⋯,b2′=2584+2412268000v4+228716123600000v6+1127781372728063800000000v8+1359935168355857106305000000000v10+⋯,b1′′=124+121600v4+16758320000v6+52838916533720000000v8+9673435845851000000000v10+⋯b2′′=125336-121600v4-2867408240000v6-302210916533720000000v8-1457821730731375000000v10+⋯,a^21=150+115000v2+111250000v4+115750000000v6+135437500000000v8+1116943750000000000v10+⋯,a^32=727+3712150v2+2876561000v4-68827516678750000v6+70974171674039150000000v8-39720321233621487034437500000000v10+⋯,a^43=935-433000v2+33231350000v4-111136579255150000000v6+26973882539344452500000000v8-72136159838800151151196250000000000v10+⋯,a32=2894860-67729000v2+261412755620000v4-467154139858075000000v6+22160081031104865839000000000v8-15664491766661484759886861250000000000v10+⋯,a43=7150-221135000v2+114463637875000v4-101047288934445250000000v6+22257681969742625996875000000v8-8431888341871633389770349418750000000000v10+⋯,γ2=1+115000v4+111250000v6+115750000000v8+135437500000000v10+⋯,γ3=1-59225v2-95691215000v4-74831911250000v6+2610442911481750000000v8-49593722925833937500000000v10+⋯γ4=1-v2+191000v4+22674050000v6-4858267255150000000v8+65060629344452500000000v10+⋯,γ^2=1+1150v2+175000v4+178750000v6+1141750000000v8+1389812500000000v10+⋯,γ^3=1-1270v2+1360750v4-2032717609375v6+1361092906317968750v8-27268073596568486328125v10+⋯,γ^4=1+1150v2+7915000v4-21619512126250000v6+78384961425250000000v8-5664933289891170503987500000000v10+⋯where γ1=1,γ^1=1.This results in the new method called EFMRKT5. Asv→0, the coefficients b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of the new method EFMRKT5 reduce to the coefficients of the original method RKT5. That is to say, b1(0),b2(0),b1′(0),b2′(0),b2′′(0),b2′′(0),a32(0),a42(0),a^21(0), a^32(0),a^43(0),γ2(0),γ3(0),γ4(0),γ^2(0),γ^3(0), and γ^4(0) are identical to b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of RKT5 method. Other than that, v→0, as EFMRKT5 method will have the same error constant as RKT5 method.


## 3.2. Trigonometrically Fitted MRKT Method

Exponentially fitted method leads to trigonometrically fitted method when replacingv=wh with iv and solving (8) to (9) to find ai,i-1, a^i,i-1,γi, and γ^i.(31)γi=cos⁡vci-12v2ci2-v3∑j=1i-1ai,jsin⁡vcj,(32)a^i,i-1=1-cos⁡vci-v2∑j=1i-2a^i,jcos⁡vcjv2cos⁡vci-1,(33)ai,i-1=-sin⁡vci+v·ci-v3∑j=1i-2ai,jcos⁡vcjv3cos⁡vci-1,(34)γ^i=sin⁡vci+v2∑j=1i-2a^i,jsin⁡vcjvci,i=2,…,s.Consider the same coefficients of fifth-order four-stage method developed by Fawzi et al.[15] as in Section 3.1. Solving the (31) to (34) and letting a^21, a^32, a^43, a32, a43, γ2, γ3, γ4, γ^2, γ^3, and γ^4 be free parameters will give(35)a^21=1-cosv/5v2,a^32=1-cos2v/3+1/27v2v2cosv/5,a^43=1-cosv-v23/10-2/35cosv/5v2cos2v/3,a32=-sin2v/3+2v/3+49/4860v3v3cosv/5,a43=v-sinv-v37/50-1/50cosv/5v3cos2v/3,γ2=cosv5-v250,γ3=cos2v3-2v29-301v34860sinv5,γ4=cosv-v22-v3-150sinv5+125sin2v3,γ^2=5vsinv5,γ^3=23vsin2v3+21-cos2v/3+v2/27sinv/53v2cosv/35,γ^4=sin⁡vv-2v35sin⁡v5+1-cos⁡v-3v2/10-2v2/5cos⁡v/5sin⁡2v/3vcos⁡2v/3.Next, solving (10) to (12), and using the above Fawzi coefficients to find b1, b2,b1′, b2′, b1′′, and b2′′,(36)b1=3cos⁡v/5sin⁡2v/3-3cos⁡2v/3sin⁡v/5112sin⁡v/5-2cos⁡v/5cos⁡v-2cos⁡v/5+cos⁡v/5v2+2sin⁡vsin⁡v/5-2vsin⁡v/52v3sin⁡v/5,b2=-3112sin⁡2v/3sin⁡v/5+2cos⁡v-2+v22v3sin⁡v/5,b1′=9cos⁡v/5sin⁡2v/3-9cos⁡2v/3sin⁡v/556sin⁡v/5+cos⁡v/5sin⁡v-vcos⁡v/5-cos⁡vsin⁡v/5+sin⁡v/5v2sin⁡v/5,b1′′=-25v/243sin⁡v+80/81cos⁡v-80/81+23v/48sin⁡2v/380v/81sin⁡v/5-v/125sin⁡2v/3,b2′′=-1+27/56vsin⁡2v/3+5/48vsin⁡v+cos⁡vvsin⁡v/5,b2′=-27/56cos⁡v/5vsin⁡2v/3+5/48cos⁡v/5vsin⁡v+cos⁡v/5cos⁡v-cos⁡v/5-27/56vcos⁡2v/3sin⁡v/5-5/48vcos⁡vsin⁡v/5+sinvsin⁡v/5vsin⁡v/5These lead to our new explicit trigonometrically fitted MRKT which is called TFMRKT5 method. The corresponding Taylor series expansion of the solution is given by(37)b1=148-12160v2+101136080000v4+5713183708000000v6-1133033981841914000000000v8+57722879134057055132000000000v10+⋯,b2=542+12160v2-2921136080000v4+136145927000000v6+5429358740920957000000000v8-6964030429670285275660000000000v10+⋯,b1′=124-17283500v4+149218700000v6-1055069341007975000000v8+230256892538959377500000000v10+⋯,b2′=2584+2412268000v4-228716123600000v6+1127781372728063800000000v8-1359935168355857106305000000000v10+⋯,b1′′=124+121600v4-16758320000v6+52838916533720000000v8-9673435845851000000000v10+⋯,b2′′=125336-121600v4+2867408240000v6-302210916533720000000v8+1457821730731375000000v10+⋯,a^21=150-115000v2+111250000v4-115750000000v6+135437500000000v8-1116943750000000000v10+⋯,a^32=727-3712150v2+2876561000v4+68827516678750000v6+70974171674039150000000v8+39720321233621487034437500000000v10+⋯,a^43=935+433000v2+33231350000v4+111136579255150000000v6+26973882539344452500000000v8+72136159838800151151196250000000000v10+⋯,a32=2894860+67729000v2+261412755620000v4+467154139858075000000v6+22160081031104865839000000000v8+15664491766661484759886861250000000000v10+⋯,a43=7150+221135000v2+114463637875000v4+101047288934445250000000v6+22257681969742625996875000000v8+8431888341871633389770349418750000000000v10+⋯,γ2=1-125v2+115000v4-111250000v6+115750000000v8-135437500000000v10+⋯,γ3=1-59225v2-95691215000v4-74831911250000v6+2610442911481750000000v8-49593722925833937500000000v10+⋯,γ4=1-v2+191000v4+22674050000v6-4858267255150000000v8+65060629344452500000000v10+⋯,γ^2=1-1150v2+175000v4-178750000v6+1141750000000v8-1389812500000000v10+⋯,γ^3=1+1270v2+1360750v4+2032717609375v6+1361092906317968750v8+27268073596568486328125v10+⋯,γ^4=1-1150v2+7915000v4+21619512126250000v6+78384961425250000000v8+5664933289891170503987500000000v10+⋯,whereγ1=1,γ^1=1.This results in the new method called TFMRKT5. Asv→0, the coefficients b1,b2,b1′,b2′,b2′′,b2′′,a32,a42,a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of the new method TFMRKT5 reduce to the coefficients of the original method RKT5. That is to say, b1(0),b2(0),b1′(0),b2′(0),b2′′(0),b2′′(0),a32(0),a42(0), a^21(0),a^32(0),a^43(0),γ2(0),γ3(0),γ4(0),γ^2(0),γ^3(0), and γ^4(0) are identical to b1,b2,b1′,b2′,b2′′,b2′′,a32,a42, a^21,a^32,a^43,γ2,γ3,γ4,γ^2,γ^3, and γ^4 of RKT5 method. Other than that, v→0, as TFMRKT5 method will have the same error constant as RKT5 method.


## 4. Error Analysis

In this section, we will find the principal local truncation errors fory,y′, and y′′ (i.e., τn+1,τn+1′,τn+1′′) of the new exponentially fitted and trigonometrically fitted explicit modified Runge-Kutta type methods, respectively. We first find the Taylor series expansion of the actual solution y(xn+h), the first derivative of the actual solution y′(xn+h), and the second derivative of the actual solution y′′(xn+h), the approximate solution yn+1, the first derivative of the approximate solution yn+1′, and the second derivative of the approximate solution yn+1′′. The local truncation errors of y, y′, and y′′ are given as (38)τn+1=yn+1-yxn+h,τn+1′=yn+1′-y′xn+h,τn+1′′=yn+1′′-y′′xn+hTheτn+1,τn+1′, and τn+1′′ of the methods are given in the Appendix.Notes: fromτn+1,τn+1′, and τn+1′′, we can see that the order of TFMRKT5 is order 5 because all of the coefficients up to h5 vanished.


## 5. Problems Tested and Numerical Results

In this section, we will apply the new explicit exponentially fitted modified Runge-Kutta type method to somey′′′=f(x,y,y′) ODEs for problems (1)-(4) which consist of exponential solutions and the new trigonometrically fitted modified Runge-Kutta type method to some ODEs problems (5)-(8) with trigonometric functions solutions. The numerical results are compared with the results obtained when the same set of problems are reduced to a system of first-order equations and is solved using the existing Runge-Kutta of the same order.(i)
h: step sizes.(ii)
TFMRKT5:the four-stage fifth-order trigonometrically fitted RK type method derived in this paper.(iii)
EFMRKT5:the four-stage fifth-order exponentially fitted RK type method derived in this paper.(iv)
RKT5:the four-stage fifth-order RK type method given by Fawzi et al. [15].(v)
RK5B:the six-stage fifth-order RK method given in Butcher [16].(vi)
RKF5:the six-stage fifth-order RK method given in Lambert [17].(vii)
TFRK:the six-stage fifth-order trigonometrically fitted RK method given in Anastassi et al. [18].Problem 2 (homogeneous linear problem).
(39) y ′ ′ ′ x = 2 y ′ x ,  y 0 = 0 ,   y ′ 0 = 1 ,   y ′ ′ 0 = 0 ,exact solution is (40)yx=2e2x4-2e-2x4.
Estimated frequencyω=2.Problem 3 ( homogeneous linear system).
(41) y 1 ′ ′ ′ x = 8 y 3 ′ x ,  y 1 0 = 2 ,   y 1 ′ 0 = 4 ,   y 1 ′ ′ 0 = 8 , y 2 ′ ′ ′ x = 8 y 1 ′ x ,  y 2 0 = 4 ,   y 2 ′ 0 = 8 ,   y 2 ′ ′ 0 = 16 , y 3 ′ ′ ′ x = y 2 ′ x ,  y 3 0 = 1 ,   y 3 ′ 0 = 2 ,   y 3 ′ ′ 0 = 4 ,exact solutions are (42)y1x=2e2x,y2x=4e2x,y3x=e2x.
Estimated frequencyω=2.Problem 4 ( inhomogeneous linear system ).
(43) y 1 ′ ′ ′ x = y 3 ′ x + 1 ,  y 1 0 = 2 ,   y 1 ′ 0 = 3 ,   y 1 ′ ′ 0 = 5 , y 2 ′ ′ ′ x = y 1 ′ x + 2 ,  y 2 0 = 1 ,   y 2 ′ 0 = 2 ,   y 2 ′ ′ 0 = 5 , y 3 ′ ′ ′ x = y 2 ′ x + 3 ,  y 3 0 = 0 ,   y 3 ′ 0 = 4 ,   y 3 ′ ′ 0 = 5 ,exact solutions are (44)y1x=5ex-2x-3,y2x=5ex-3x-4,y3x=5ex-x-5.
Estimated frequencyω=1.Problem 5 (inhomogeneous linear problem).
(45) y ′ ′ ′ x = 5 y ′ x + sinh ⁡ x ,  y 0 = - 1 4 ,   y ′ 0 = 0 ,   y ′ ′ 0 = - 1 4 ,exact solution is (46)yx=-ex8-e-x8.
Estimated frequencyω=1.Problem 6 (homogeneous linear problem).
(47) y ′ ′ ′ x = - 25 y ′ x ,  y 0 = 0 ,   y ′ 0 = 0 ,   y ′ ′ 0 = 1 ,exact solution is (48)yx=125-125cos⁡5x.
Estimated frequencyω=5.Problem 7 (inhomogeneous linear problem).
(49) y ′ ′ ′ x = - 27 y ′ x + sin ⁡ x ,  y 0 = 1 ,   y ′ 0 = - 1 ,   y ′ ′ 0 = 0 ,exact solution is (50)yx=3702cos⁡33x-39sin⁡33x-126cos⁡x+2827.
Estimated frequencyω=33.Problem 8 (inhomogeneous linear system).
(51) y 1 ′ ′ ′ x = - 27 y 1 ′ x ,  y 1 0 = 0 ,   y 1 ′ 0 = - 1 ,   y 1 ′ 0 = 0 , y 2 ′ ′ ′ x = - 27 y 2 ′ x + cos ⁡ x ,  y 2 0 = 1 ,   y 2 ′ 0 = - 1 ,   y 2 ′ ′ 0 = 0 , y 3 ′ ′ ′ x = - 27 y 3 ′ x ,  y 3 0 = 1 ,   y 3 ′ 0 = 0 ,   y 3 ′ ′ 0 = - 1 ,exact solutions are (52)y1x=-39sin⁡33x,y2x=1-3326sin⁡33x+126sin⁡x,y3x=2627+127cos⁡33x.
Estimated frequencyω=33.Problem 9 (inhomogeneous linear system).
(53) y 1 ′ ′ ′ x = - 7 y 2 ′ x - cos ⁡ x ,  y 1 0 = 1 7 ,   y 1 ′ 0 = 0 , y 1 ′ ′ 0 = - 1 , y 2 ′ ′ ′ x = - 7 y 1 ′ x - cos ⁡ x ,  y 2 0 = 1 ,   y 2 ′ 0 = 0 ,   y 2 ′ ′ 0 = - 1 , y 3 ′ ′ ′ x = - 7 y 3 ′ x - cos ⁡ x ,  y 3 0 = 0 ,   y 3 ′ 0 = 1 ,   y 3 ′ ′ 0 = 0 ,exact solutions are (54)y1x=-16sin⁡x+742sin⁡7x+17cos⁡7x,y2x=-16sin⁡x+742sin⁡x7+67+17cos⁡7x,y3x=76sin⁡7x-16sin⁡x.
Estimated frequencyω=1.


## 6. An Application to a Problem in Thin Film Flow

Here, we will use the suggested method to a famous problem in engineering and physics based on the thin film flow of a liquid. Many researchers in the literature explain this problem more. Momoniat and Mahomed[19] constructed symmetry reduction and numerical solution of a third-order ODE from thin film flow. Tuck and Schwartz [20] discussed the movement of a thin film of viscous fluid over a solid surface and taken into account tension and gravity, as well as viscosity. The problem was evaluated and solved using third-order ODE as follows: (55)d3ydx3=fy.Many forms of the function were studied by [20] for the drainage dry surface; it has the form of f(y) which can be stated as(56)d3ydx3=-1+1y2.When the surface is prewetted by a thin film with thickness δ>0 (where δ>0 is very small), the function f is given by(57)fy=-1+1+δ+δ2y2-δ+δ2y3Problems concerning the flow of thin films of viscous fluid with a free surface in which surface tension effects play a role typically lead to third-order ODEs governing the shape of the free surface of the fluid,y=y(x). As indicated by [20], one such equation is(58)y′′′x=y-k,x≥x0with initial conditions(59)yx0=y0,y′x0=y0′,y′′x0=y0′′,where y0, y0′, and y0′′ are constants, which is of specific significance since it portrays the dynamic balance amongst surface and gooey strengths in a thin fluid layer in disregard of gravity. For compare and contrast, we utilized Runge-Kutta methods which are fifth-order (RKT5, RK5B, RKF5, and TFRKT) strategies, individually. To utilize Runge-Kutta techniques we write (1) as a system of three first-order equations. Biazar et al. [21] we can write (58) as the following system:(60)dy1dx=y2x,dy2dx=y3x,dy2dx=y1-k,where(61)y10=1,y20=1,y30=1.We have taken x0=0 and y0=y0′=y0′′=1. Unfortunately, for general k, (58) cannot be solved analytically. However, we can use these reductions to determine an efficient way to solve (1) numerically. Here, we are focusing on the cases k=2 and k=3 (see Mechee et al.[22]).


## 7. Discussion and Conclusion

In this research, we have derived exponentially fitted and trigonometrically fitted explicit modified Runge-Kutta type methods for solvingy′′′(x)=f(x,y,y′) with application to thin film flow problem. Consequently, the new four-stage fifth-order exponentially-fitted and trigonometrically-fitted methods which are denoted as EFMRKT5 and TFMRKT5, respectively, were constructed and we used in numerical comparison the criteria based on computing the maximum error in the solution max⁡ytn-yn which is equal to the maximum between absolute errors of the actual solutions and computed solutions. The numerical outcomes are plotted in Figures 1–8. Figures 1–8 demonstrate that the new TFMRKT5 and EFMRKT5 methods require less capacity assessments than the RKT5, RK5B, RKF5, and TFRK methods. The figures showed the efficiency of the new methods where the common logarithm of the maximum global error throughout the integration versus computational cost was measured by the number of function evaluations. The numerical results obtained showed clearly that the global error for a short period of integration for the new exponentially fitted method and for a large period of integration for the new trigonometrically fitted explicit modified Runge-Kutta type method is smaller than that of the other existing methods. The new EFMRKT5 and TFMRKT5 methods are much more efficient than the other existing methods when solving third-order ODEs of the form y′′′=f(x,y,y′) straightforwardly. For Tables 2 and 3 we observed that the numerical results using TFMRKT5 and EFMRKT5 methods are correct to five decimal places. Applying RK5B, RKF5, TFRK, and RKT5 to (58) for k=2 also yields five-decimal place accuracy. Tables 4 and 5 show the numerical results for the case k=3 with h=0.1 and h=0.01 since for k=3, Problem (58) cannot be solved analytically. Table 4 shows that TFMRKT5 and EFMRKT5 manage to achieve the numerical results which agree to seven decimal places when compared to RK5B, RKF5, TFRK, and RKT5 for h=0.1. In Table 5 the numerical results for TFMRKT5 and EFMRKT5 agree to nine decimal places when compared to RK5B, RKF5, TFRK, and RKT5 for h=0.01. For Table 7 we observe that RK5B, RKF5, RKT5, TFRK, TFMRKT5, and EFMRKT5 have similar order of accuracy. In Table 6 values of the error are different. Therefore it is consistent with results displayed in Tables 2 and 3. Figures 9 and 10 show that the new EFMRKT5 and TFMRKT5 methods require less function evaluations than the RK5B, RKF5, TFRK, and RKT5 methods. This is because when problem (58) is solved using RK5B, RKF5, TFRK, and RKT5 methods, it needs to be reduced to a system of first-order equations which is three times the dimension.Table 2
Numerical results for problem in Thin Film Flow (58) taking h=0.1 and k=2.
x Exact Solution RK5B RKF5 RKT5 EMFRKT5 TFRK TFMRKT5 0.0 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 0.2 1.221211030 1.2212100068 1.2212100097 1.2212100039 1.2212100052 1.2212100218 1.2212100052 0.4 1.488834893 1.4888347851 1.4888347895 1.4888347797 1.4888347885 1.4888348090 1.4888347885 0.6 1.807361404 1.8073614063 1.8073614114 1.8073613988 1.8073614237 1.8073614357 1.8073614237 0.8 2.179819234 2.1798192463 2.1798192513 2.1798192371 2.1798192873 2.1798192788 2.1798192873 1.0 2.608275822 2.6082748841 2.6082748883 2.6082748735 2.6082749587 2.6082749176 2.6082749587Table 3
Numerical results for problem in thin film flow (58) taking h=0.01 and k=2.
x Exact Solution RK5B RKF5 RKT5 EFMRKT5 TFRK TFMRKT5 0.0 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 0.2 1.221211030 1.2212100045 1.2212100045 1.2212100045 1.2212100045 1.2212100045 1.2212100045 0.4 1.488834893 1.4888347799 1.4888347799 1.4888347799 1.4888347799 1.4888347799 1.4888347799 0.6 1.807361404 1.8073613977 1.8073613977 1.8073613977 1.8073613977 1.8073613977 1.8073613977 0.8 2.179819234 2.1798192339 2.1798192339 2.1798192339 2.1798192340 2.1798192339 2.1798192340 1.0 2.608275822 2.6082748676 2.6082748676 2.6082748676 2.6082748677 2.6082748676 2.6082748677Table 4
Numerical results for problem in thin film flow (58) taking h=0.1 and k=3.
x RK5B RKF5 RKT5 EMFRKT5 TFRK TFMRKT5 0.0 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 0.2 1.2211551491 1.2211551546 1.2211551394 1.2211551412 1.2211551831 1.2211551412 0.4 1.4881052974 1.4881053065 1.4881052807 1.4881052926 1.4881053519 1.4881052926 0.6 1.8042625677 1.8042625794 1.8042625459 1.8042625786 1.8042626364 1.8042625786 0.8 2.1715228242 2.1715228376 2.1715227987 2.1715228633 2.1715229031 2.1715228633 1.0 2.5909582923 2.5909583063 2.5909582638 2.5909583715 2.5909583783 2.5909583715Table 5
Numerical results for problem in thin film flow (58) taking h=0.01 and k=3.
x RK5B RKF5 RKT5 EMFRKT5 TFRK TFMRKT5 0.0 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 1.000000000 0.2 1.2211551424 1.2211551424 1.2211551424 1.2211551424 1.2211551424 1.2211551424 0.4 1.4881052842 1.4881052842 1.4881052842 1.4881052842 1.4881052842 1.4881052842 0.6 1.8042625481 1.8042625481 1.8042625481 1.8042625482 1.8042625481 1.8042625482 0.8 2.1715227981 2.1715227981 2.1715227981 2.1715227982 2.1715227981 2.1715227982 1.0 2.5909582591 2.5909582591 2.5909582591 2.5909582592 2.5909582591 2.5909582592Table 6
Comparison of error for problem in thin film flow (58) taking h=0.1 and k=2.
x RK5B RKF5 RKT5 EFMRKT5 TFRK TFMRKT5 0.0 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.2 1.0230 (-6) 1.0200(-6) 1.2600(-6) 1.0250(-6) 1.0080(-6) 1.0250(-6) 0.4 1.0800(-7) 1.0300(-7) 1.1300(-7) 1.0500(-7) 8.4100(-7) 1.0500(-7) 0.6 2.0000(-9) 7.0000(-9) 5.0000(-8) 2.0000(-8) 3.2000(-8) 2.0000(-8) 0.8 1.2000(-8) 1.7000(-8) 3.0000(-9) 5.300(-8) 4.5000(-8) 5.3000(-8) 1.0 9.3800(-7) 9.3400(-7) 9.4800 (-7) 8.6300 (-7) 9.0400(-7) 8.6300(-7)Table 7
Comparison of error for problem in thin film flow (58) taking h=0.01 and k=2.
x RK5B RKF5 RKT5 EFMRKT5 TFRK TFMRKT5 0.0 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.2 1.0260(-6) 1.0260(-6) 1.0260(-6) 1.0260(-6) 1.0260(-6) 1.0260(-6) 0.4 6.0000(-7) 6.0000(-7) 6.0000(-7) 6.0000(-7) 6.0000(-7) 6.0000(-7) 0.6 9.0000(-9) 9.0000(-9) 9.0000(-9) 9.0000(-9) 9.0000(-9) 9.0000(-9) 0.8 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 0.0000(0) 1.0 9.5400(-7) 9.5400(-7) 9.5400(-7) 9.5400(-7) 9.5400(-7) 9.5400(-7)Figure 1
The efficiency curve for EFMRKT5, RKT5, RK5B, and RKF5 for Problem2 with xend=5 and h=0.1,0.25,0.5,0.75.Figure 2
The efficiency curve for EFMRKT5, RKT5, RK5B, and RKF5 for Problem3 with xend=2 and h=0.1,0.25,0.5,0.75.Figure 3
The efficiency curve for EFMRKT5, RKT5, RK5B, and RKF5 for Problem4 with xend=6 and h=0.1,0.25,0.5,0.75.Figure 4
The efficiency curve for EFMRKT5, RKT5, RK5B, and RKF5 Problem5 with xend=4 and h=0.1,0.25,0.5,0.75.Figure 5
The efficiency curve for TFMRKT5, RKT5, RK5B, RKF5, and TFRK for Problem6 with xend=10000 and h=0.025,0.05,0.075,0.1.Figure 6
The efficiency curve for TFMRKT5, RKT5, RK5B, RKF5, and TFRK for Problem7 with xend=10000 and h=0.025,0.05,0.075,0.1.Figure 7
The efficiency curve for TFMRKT5, RKT5, RK5B, RKF5, and TFRK for Problem8 with xend=10000 and h=0.025,0.05,0.075,0.1.Figure 8
The efficiency curve for TFMRKT5, RKT5, RK5B, RKF5, and TFRK for Problem9 with xend=10000 and h=0.025,0.05,0.075,0.1.Figure 9
Plot of graph for function evaluations against step sizeh for Problem (58) taking xend=1,h=0.1,0.01,0.001,0.0001, and k=2.Figure 10
Plot of graph for function evaluations against step sizeh for Problem (58) taking xend=1,h=0.1,0.01,0.001,0.0001 and k=2.



---

*Source: 4029371-2018-07-04.xml*