
Citation: Amit Prakash, Manish Goyal, Haci Mehmet Baskonus, Shivangi Gupta. A reliable hybrid numerical method for a time dependent vibration model of arbitrary order[J]. AIMS Mathematics, 2020, 5(2): 979-1000. doi: 10.3934/math.2020068
[1] | M. Mossa Al-Sawalha, Azzh Saad Alshehry, Kamsing Nonlaopon, Rasool Shah, Osama Y. Ababneh . Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm. AIMS Mathematics, 2022, 7(11): 19739-19757. doi: 10.3934/math.20221082 |
[2] | Shabir Ahmad, Aman Ullah, Ali Akgül, Fahd Jarad . A hybrid analytical technique for solving nonlinear fractional order PDEs of power law kernel: Application to KdV and Fornberg-Witham equations. AIMS Mathematics, 2022, 7(5): 9389-9404. doi: 10.3934/math.2022521 |
[3] | Manal Alqhtani, Khaled M. Saad, Rasool Shah, Thongchai Botmart, Waleed M. Hamanah . Evaluation of fractional-order equal width equations with the exponential-decay kernel. AIMS Mathematics, 2022, 7(9): 17236-17251. doi: 10.3934/math.2022949 |
[4] | Aslı Alkan, Halil Anaç . A new study on the Newell-Whitehead-Segel equation with Caputo-Fabrizio fractional derivative. AIMS Mathematics, 2024, 9(10): 27979-27997. doi: 10.3934/math.20241358 |
[5] | Aslı Alkan, Halil Anaç . The novel numerical solutions for time-fractional Fornberg-Whitham equation by using fractional natural transform decomposition method. AIMS Mathematics, 2024, 9(9): 25333-25359. doi: 10.3934/math.20241237 |
[6] | Maysaa Al Qurashi, Saima Rashid, Sobia Sultana, Fahd Jarad, Abdullah M. Alsharif . Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory. AIMS Mathematics, 2022, 7(7): 12587-12619. doi: 10.3934/math.2022697 |
[7] | Reny George, Sina Etemad, Ivanka Stamova, Raaid Alubady . Existence of solutions for $ [\mathtt{p},\mathtt{q}] $-difference initial value problems: application to the $ [\mathtt{p},\mathtt{q}] $-based model of vibrating eardrums. AIMS Mathematics, 2025, 10(2): 2321-2346. doi: 10.3934/math.2025108 |
[8] | Sumbal Ahsan, Rashid Nawaz, Muhammad Akbar, Saleem Abdullah, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Numerical solution of system of fuzzy fractional order Volterra integro-differential equation using optimal homotopy asymptotic method. AIMS Mathematics, 2022, 7(7): 13169-13191. doi: 10.3934/math.2022726 |
[9] | Muath Awadalla, Abdul Hamid Ganie, Dowlath Fathima, Adnan Khan, Jihan Alahmadi . A mathematical fractional model of waves on Shallow water surfaces: The Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10561-10579. doi: 10.3934/math.2024516 |
[10] | Laiq Zada, Rashid Nawaz, Sumbal Ahsan, Kottakkaran Sooppy Nisar, Dumitru Baleanu . New iterative approach for the solutions of fractional order inhomogeneous partial differential equations. AIMS Mathematics, 2021, 6(2): 1348-1365. doi: 10.3934/math.2021084 |
The vibration of membranes plays a significant role in analysis of wave mechanics in two dimensions and wave propagation, bio-engineering etc. Membranes create main components in acoustics and music such as components of microphones, speakers and related devices [1]. To investigate design of hearing aids, the knowledge of large membrane vibration [2] is vital. In bio-engineering, several human tissues are anticipated as membranes. The vibrational features of eardrum is valuable to understand hearing. The equation of vibration is used to designate vibration of membranes [3]. The integer order vibration equation [4] is,
1c2∂2 ω (r,t)∂t2=∂2 ω (r,t)∂r2+1r∂w(r,t)∂r;r≥0,t≥0, | (1) |
where ω(r,t) signifies probability density function of particle [5] at time t at a position r while c is the wave velocity of vibrations.
Fractional order derivatives offer a superb mechanism for describing hereditary and memory related properties of various real-life processes and materials [6]. The analysis of fractional differential equations (FDE) [7,8,9,10] in Mathematical Physics, vibration, oscillation, signal processing [11], visco-elasticity [12], chemical engineering [13], seismic wave propagation [14], modelling of diseases [15] etc. is a growing field of interest for researchers. In literature, there exist operational matrix method [3], decomposition method [4], homotopy perturbation scheme [5], variational iteration technique [16] etc. to solve vibration equation.
Homotopy analysis method (HAM) offers an easier way to confirm convergence of solution. HAM was introduced by Liao [17] for solving differential equations. El-Tavil and Huseen [18,19] proposed q-homotopy analysis technique (q-HAM) by generalization of HAM. But these methods have limitations like massive computation with more time consumption. So, they necessitate to be linked with a transform operator. Hybrid methods using integral transforms [20,21] are useful to treasure some solution of nonlinear FDE. Homotopy analysis Laplace transform method (HATM) is a united form of HAM and the transform of Laplace. In [1], Srivastava et al. used HATM and Laplace decomposition technique to solve vibration equation of arbitrary order. The reliability of solution procedure of a nonlinear equation is an important characteristic than modeling dimensions of equations [22,23,24].
Sumudu transform has an interesting advantage of the ‘unity’ feature over Laplace transform. It leads combinations into permutations hence it is useful in discrete systems. The function along with its Sumudu transform possess identical Taylor coefficients other than factor n. Sumudu transform is used due to its strong properties to get solution of many problems [25,26]. Watugala [27,28] proposed Sumudu transform and Asiru [29] proved its properties. Weerakoon [30,31] discussed its applications in finding solution to wave equation with variable coefficients.
Homotopy analysis Sumudu transform method (HASTM) is a graceful merger of HAM and Sumudu transform. The benefit is its power of embracing two robust computational schemes for tackling FDE. The projected approaches can reduce time and computation work as compared to existing schemes simultaneously preserving result efficiency. Singh et al. [32] used HASTM to solve fractional Drinfeld-Sokolov-Wilson equation. The q-homotopy analysis Sumudu transform method (q-HASTM) is an improvement of q∈[0,1] in HASTM to parameter q∈[0,1n],n≥1 The existence of (1n)m in solution helps in converging quickly.
Our aim is to investigate the fractional model of vibration Equation (1) and get its numerical solution by q-HASTM. The paper is presented as follows. Section 1 is introductory. In section 2, basic results of derivative in Caputo sense, Sumudu transform and its properties are provided. In section 3, mathematical model of time dependent vibration equation of fractional order is discussed along with its necessity and our motivation in finding its solution. In section 4, basic idea of q-HASTM is provided. In section 5, its implementation on fractional vibration equation is shown with convergence analysis. In section 6, we conduct numerical experiments by taking various initial conditions. In section 7, numerical results are discussed using figures and tables while in section 8, we summarize conclusion.
Definition 2.1. [33] A real function h(χ), χ>0 is said to be in spaces
a. Cζ,ζ∈R if there exists a real number q (>ζ), s.t. h(χ)=χqh1(χ),h1(χ)∈C[0,∞). Clearly, Cζ⊂Cγ if γ≤ζ.
b. Cmζ,m∈N∪{0} if h(m)∈Cζ.
Definition 2.2. [33] Fractional derivative in Caputo sense [34] of h(t),h(t)∈Cm−1,m∈N∪{0} is:
Dβth(t)={Im−βh(m)(t),m−1<β<m,m∈N,dmdtmh(t),β=m, |
a.[35] Iζth(x,t)=1Γζ∫t0(t−s)ζ−1h(x,s)ds;ζ,t>0.
b.[35] DvτV(x,τ)=Im−vτ∂mV(x,τ)∂τm,m−1<v≤m
c.[35] DζtIζth(t)=h(t),m−1<ζ≤m,m∈N.
d.[35] IζtDζth(t)=h(t)−∑m−1k=1h(k)(0+)tkk!,m−1<ζ≤m,m∈N.
e.[35] Iβtα=Γ(α+1)Γ(β+α+1)tβ+α.
Definition 2.3. [36] Sumudu transform on domain of functions:
Q={h(p)|∃N,s1,s2>0,|h(p)|<Nepsj if p∈(−1)j×[0,∞), |
is in the form of S[h(p)]=∫∞0h(wp)e−pdp,w∈(−s1,s2).
Definition 2.4. [36] Sumudu transform for arbitrary order derivative in Caputo sense is:
S[Dnβxω(x,t)]=s−nβS[ω(x,t)]−n−1∑k=0s(−nβ+k)ω(k)(0,t),n−1<nβ≤n. |
Many physical quantities are concerned with the past so to understand their physical models better by inducing the effects of memory, fractional models of such systems get more importance [3]. The FDE accomplish the systems with memory effect. The non-local property is the main benefit of working with FDE in physical models. It signifies that future system state is dependent on former states also. Thus, the models having fractional order derivatives adhere to reality. The integer order model suggested in [4] was found unable to possess memory effect in vibrational motion, so to include these effects, the integer order model is generalized to model of arbitrary order by converting derivative of integer order to fractional order in Caputo sense. The differential equation of arbitrary order [1,5,6,16] for vibration model is given as,
1c2∂αω(r,t)∂tα=∂2ω(r,t)∂r2+1r∂w(r,t)∂r,1<α≤2, | (2) |
with initial settings:
w(r,0)=φ(r),∂w(r,0)∂t=cξ(r),0≤r≤1,0≤t≤1. | (3) |
Caputo’s derivative is suitable for differentiable functions [34] and permits conditions to comprise in modelling a problem. The general response expression contains a parameter that states the arbitrary order of the derivative. It can be varied to get different responses [6]. In case α=2 Eq. (2) reduces to integer order Eq. (1).
It is also observed that the total hierarchy of moments Mk=⟨rk(t)⟩ possess similar time dependence as arbitrary Brownian motion in spite of little difference in their statistical features. In [6], power law decay of solution is found with α in disparity to exponential decay perceived in arbitrary Brownian motion. The time fractional equations depict the particle motion with memory in time. Time fractional derivative proposes inflection of memory. It is obvious that vibrational motion is affected by memory in time. It scripts suitability of fractional modeling for this system. So, the comprehensive study of Eq. (2) to find the numerical solution of this mathematical model is very important. It motivated us to solve Eq. (2) by an efficient and novel numerical scheme q-HASTM.
Consider a general fractional order non-linear partial differential equation of the form:
Pαu(x,t)=h(x,t), | (4) |
where Pα signifies the general fractional linear and nonlinear partial differential operator, u(x,t) is an unknown function of independent variables x and t Linear terms of Pα are decomposed to Dα+R where Dα is linear operator of the highest order. (Dαu)(t) is Caputo fractional derivative of u(x,t). R is remains of the linear operator. Eq. (4) may be written as,
(Dαu)(t)+Ru(x,t)+Nu(x,t)=h(x,t),(α>0,n−1<α≤n(n∈N)), | (5) |
where Nu shows non-linear terms.
Applying Sumudu transform [37] in Eq. (5), we get,
S[(Dαu)(t)]+S[Ru(x,t)]+S[Nu(x,t)]=S[h(x,t)], | (6) |
Here, S is Sumudu transform operator.
Using property [37] of Sumudu transform, we find,
S[u(x,t)]−pα∑n−1k=0u(k)(x,0)pα−k+pα(S[Ru(x,t)]+S[Nu(x,t)]−S[h(x,t)])=0, | (7) |
We state nonlinear operator as,
N[ξ(x,t;q)]=S[ξ(x,t;q)]−pα∑n−1k=0ξ(r)(x,t;q)(0)p(α−r)+pα(S[Rξ(x,t;q)]+S[Nξ(x,t;q)]−S[h(x,t)]) | (8) |
Here, q∈[0,1n];n≥1 is an embedding parameter. ξ(x,t;q) is a real valued function of x,q and t
Homotopy is constructed as:
(1−q)S[ξ(x,t;q)−u0(x,t)]=ℏqH(x,t)Pα[u(x,t)], | (9) |
Here, H≠0,ℏ≠0.H is auxiliary function and ℏ is auxiliary parameter. u0(x,t) is an initial guess of u(x,t) and ξ(x,t;q) is unknown function.
By choosing q=0 and 1n;n≥1 in Eq. (9), we get,
ξ(x,t;0)=u0(x,t) and ξ(x,t;1n)=u(x,t). | (10) |
As q surges from 0 to 1n,n≥1,ξ(x,t;q) varies from u0(x,t) to u(x,t).
Expanding ξ(x,t;q) in Taylor’s series about q, we get,
ξ(x,t;q)=u0(x,t)+∑∞m=1um(x,t)qm, | (11) |
where,
um(x,t)=1m![∂mξ(x,t;q)∂qm]q=0. | (12) |
If u0(x,t), auxiliary linear operator, ℏ and H(x,t) are suitably selected, series (11) converge at q=1n. Then, we find,
u(x,t)=u0(x,t)+∑∞m=1um(x,t)(1n)m | (13) |
Express vectors as,
→um={u0(x,t),u1(x,t),………um(x,t)}. | (14) |
Differentiating with regard to q, Eq. (9) m times then dividing by m! and taking q=0, we find deformation equation of m order :
S[um(x,t)−kmwm−1(x,t)]=ℏH(x,t)æm(→wm−1). | (15) |
Taking inverse Sumudu transform in Eq. (15),
um(x,t)=kmum−1(x,t)+ℏS−1[H(x,t)æm(→um−1)]. | (16) |
Here,
æm(→um−1)=1m−1![∂m−1N{ξ(x,t;q)}∂qm−1]q=0, | (17) |
and kr={0,r≤1,n,r>1. | (18) |
Hence, the q-HASTM solution is offered as,
u(r,t)=∑∞m=1um(x,t)(1n)m. | (19) |
We take the fractional order vibration model discussed in section 3 as,
1c2∂αω(r,t)∂tα=∂2ω(r,t)∂r2+1r∂w(r,t)∂r,1<α≤2, | (2) |
with initial settings,
w(r,0)=φ(r),∂w(r,0)∂t=cξ(r),0≤r≤1,0≤t≤1. | (3) |
Applying Sumudu transform, we attain,
S[w(r,t)]−pα∑n−1k=0w(k)(0)p(α−k)−pα[S{c2(∂2ω(r,t)∂r2+1r∂w(r,t)∂r)}]=0. | (20) |
Nonlinear operator is
N[ξ(r,t;q)]=S[ξ(r,t;q)]−pα∑n−1k=0ξ(k)(r,t;q)(0)p(α−k)−pα[S[c2(∂2ξ(r,t;q)∂r2+1r∂ξ(r,t;q)∂r)]]. | (21) |
The homotopy is,
(1−nq)S[ξ(r,t;q)−w0(r,t)]=ℏqH(r,t)N[ξ(r,t;q)], | (22) |
For q=0,ξ(r,t;0)=w0(r,t),
and, q=1n′ξ(r,t;1n)=w(r,t).
As discussed in Section 4, the q-HASTM solution will be obtained as,
w(r,t)=∞∑m=1wm(r,t)(1n)m. |
Theorem[25]. If there exists χ (a constant) as 0<χ<1 such that,
∥ζp+1(r,t)∥≤χ∥ζp(r,t)∥ for all p and
if truncated series ∑kp=0ζp(r,t) is taken as estimated solution ζ(r,t), maximum absolute truncation error is
‖ζ(r,t)−k∑p=0ζp(r,t)‖≤χk+11−χ‖ζ0(r,t)‖. |
Proof. We have,
‖ζ(r,t)−k∑p=0ζp(r,t)‖=‖∞∑p=k+1ζp(r,t)‖≤∞∑p=k+1‖ζp(r,t)‖≤∞∑p=k+1χm‖ζ0(r,t)‖ |
≤ χ p+1[1+( χ )1+( χ )2+⋯]|| ζ 0(r,t)‖≤ χ k+11− χ ‖ ζ 0(r,t)‖, |
that proves the theorem.
In this section, applicability of q-HASTM is illustrated via some test examples.
Test Example 1. Taking initial condition w(r, 0) = r2+ctr in Eq. (2) and using Sumudu transform, we get,
S [ w(r, t)]−(r2u+cru2)u-u α [S{c2(∂2 ω (r, t)∂r2+1r∂w(r,t)∂r)}]=0. | (23) |
Nonlinear operator is,
N[ φ (r,t;q)]=S[ φ (r,t;q)]−(1−kmn)(r2u+cru2)u−uα[S{c2(∂2 φ (r,t;q)∂r2+1r∂φ(r,t;q)∂r)}]. |
Deformation equation for H(r,t)=1 is,
S[wm(r,t)−kmwm−1(r,t)]=ℏæm(→wm−1), | (24) |
where, æm(→wm−1)=S[wm−1]−(1−kmn)(r2u+cru2)u−uα[S{c2(∂2wm−1∂r2+1r∂wm−1∂r)}].
Taking inverse Sumudu transform in Eq. (24), we get,
wm(r,t)=kmwm−1(r,t)+ℏS−1[Hæm(→wm−1)]. | (25) |
Simplification gives,
w0(r,t)=r2+ctr,w1(r,t)=−4c2ℏtαΓ(1+α)−c3ℏtα+1rΓ(2+α), |
w2(r,t)=−4c2ℏntαΓ(1+α)−c3ℏntα+1rΓ(2+α)−4c2ℏ2tαΓ(1+α)−c3ℏ2tα+1r(2+α)−c5ℏ2t2α+1r3Γ(2+2α)+⋯, |
and so on.
Hence, subsequent iterations wm(r,t),m≥3 can be computed using Maple software package.
The solution is
w(r,t)=w0(r,t)+∞∑m=1wm(r,t)(1n)m. |
Test Example 2. Taking w(r, 0) = r+c t r in Eq. (2) and using Sumudu transform, we get,
S[w(r,t)]−(ru+cru2)u−uα[S[c2(∂2ω(r,t)∂r2+1r∂w(r,t)∂r)]]=0. | (26) |
Also, N[ φ (r,t;q)]=S[ φ (r,t;q)]−(1−kmn)(ru+cru2)u−u α [S{c2(∂2 φ (r,t;q)∂r2+1r∂ φ (r,t;q)∂r)}].
Deformation equation for H(r,t)=1 is,
S[wm(r,t)−kmwm−1(r,t)]=ℏæm(→wm−1), | (27) |
where, æm(→wm−1)=S[wm−1]−(1−kmn)(r+ucr)−uα[S{c2(∂2wm−1∂r2+1r∂wm−1∂r)}].
Taking inverse Sumudu transform in Eq. (27), we get,
wm(r,t)=kmwm−1(r,t)+ℏS−1[æm(→wm−1)]. | (28) |
Simplification yields,
w0(r,t)=r+ctr,w1(r,t)=−c2ℏtαrΓ(1+α)−c3ℏtα+1rΓ(2+α), |
w2(r,t)=c2ℏtαr3[−(h+n)r2(1+ct+α)Γ(2+α)−c2ℏtα(1+ct+2α)Γ(2+2α)]+⋯, |
and so on.
Hence, subsequent iterations wm(r,t),m≥3 can be found.
The solution is,
w(r,t)=w0(r,t)+∞∑m=1wm(r,t)(1n)m. |
Test Example 3. Taking w(r,0)=√r+ct√r in Eq. (2) and using Sumudu transform, we get,
S[w(r,t)]−(√r+cu√r)−uα[S[c2(∂2ω(r,t)∂r2+1r∂w(r,t)∂r)]]=0. | (29) |
N[ φ (r,t;q)]=S[ φ (r,t;q)]−(1−kmn)(√r+cu√r)−u α [S{c2(∂2 φ (r,t;q)∂r2+1r∂ φ (r,t;q)∂r)}], |
Deformation equation for H(r,t)=1 is,
S[wm(r,t)−kmwm−1(r,t)]=ℏæm(→wm−1), | (30) |
where, æm(→wm−1)=S[wm−1]−(1−kmn)(√r+cu√r)−uα[S{c2(∂2wm−1∂r2+1r∂wm−1∂r)}].
Taking inverse Sumudu transform in Eq. (30), we get,
wm(r,t)=kmwm−1(r,t)+ℏS−1[æm(→wm−1)], | (31) |
Simplification yields:
w0(r,t)=√r+ct√r, |
w1(r,t)=−c2ℏtα(r+ct+rα)4r5/2Γ(2+α), |
w2(r,t)=c2ℏtα16r9/2[−4(h+n)r2(r+ct+rα)Γ(2+α)−c2ℏtα(25ct+9r(1+2α))Γ(2+2α)]+⋯, |
and so on.
Hence, subsequent iterations wm(r,t),m≥3 can be found.
The solution is, w(r,t)=w0(r,t)+∑∞m=1wm(r,t)(1n)m.
Figures 1, 2 and 3 show behavior of numerical solution w(r,t) of time fractional vibration Eq. (2) acquired at α=2 by using q-HASTM for examples 1, 2 and 3 respectively. They have been drawn for α=2 to show the nature of the unknown exact solution of the vibration model. Figure 4 shows the plots of solution w(r, t) of Eq. (2) by q-HASTM at α=1.2,1.6 and 2 taking c=5,r=5,ℏ=−1,n=1 for examples 1, 2 and 3 respectively. They reveal that w(r, t) increases with increasing t but decreases as arbitrary order α increases. This is in total agreement with the point discussed in section 3. In Figures 5, 8 and 11, behaviour of solution w(r, t) Vs. ħ at distinct α is shown. In Figures 6, 9 and 12, behaviour of w(r, t) Vs. t at diverse values of w(r, t) are carried out for ℏ for examples 1, 2, 3 similarly. Distinct ħ are taken to make residual error small and guarantee solution convergence. In Figures 7, 10 and 13, plots of w(r, t) Vs. n at different values of α for c = 5 = r, ℏ = -1, t = 0.2 are shown for examples 1, 2 and 3 respectively. We notice that by increasing n, w(r, t) increases slowly but decreases with increasing order α The legitimacy of solution in convergence region is seen through ℏ and n -curves. Figure 14 illustrates the comparison of numerical solution obtained by q-HASTM and methods in [1,5] at α=1.5 for example 2. Figure 15 depicts the absolute error between successive approximations at α=1.5 in example 2. It is clear from Figure 15 that there is a sharp decrease in the error between approximations in example 2 which confirms that the obtained solution in example 2 is convergent. Figure 16 shows the comparison of numerical solution by q-HASTM and methods in [3,4] at α=1.5 for example 3. Figure 17 illustrates the absolute error between successive approximations at α=1.5 in example 3. It is clear from Figure 17 that there is a sharp decrease in the error between approximations in example 2 which again confirms that the obtained solution in example 3 is convergent. Figure 18 shows the comparison of solution by q-HASTM and methods in [1,3,4,5,6,16] at α=1.5 for example 1. Figure 19 depicts the absolute error between successive approximations at α=1.5 in example 1 which clearly indicates that the obtained solution in example 1 is convergent. Also, the tabular comparison of results with published work are performed in Tables 1, 3 and 5 at distinct values of order α. Tables 2, 4 and 6 show that error between successive approximations is negligible and becomes zero as iterations increase. Hence, it is concluded that the q-HASTM work also for those physical models which have derivatives of arbitrary order and have no exact solution.
r | t | α | Solution by q-HASTM | Solution by operational Matrix method [3] | Solution by Methods in [1,4,5,6,16] |
0.2 | 1.5 | 0.0467182 | 0.0469 | 0.0467 | |
2 | 0.0448067 | 0.0448 | 0.0448 | ||
0.4 | 1.5 | 0.183688 | 0.1838 | 0.1837 | |
2 | 0.179227 | 0.1792 | 0.1792 | ||
0.6 | 1.5 | 0.410124 | 0.4100 | 0.4101 | |
2 | 0.40326 | 0.4032 | 0.4032 | ||
0.8 | 1.5 | 0.725746 | 0.7256 | 0.7257 | |
2 | 0.716907 | 0.7169 | 0.7169 | ||
1 | 1.5 | 1.13039 | 1.1310 | 1.1303 | |
2 | 1.12017 | 1.1203 | 1.1202 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 8.59×10-5 | 8.77×10-9 | 2.66×10-5 | 5.46×10-10 | |
0.4 | 2.43×10-4 | 1.72×10-8 | 8.97×10-5 | 1.54×10-9 | |
0.6 | 4.46×10-4 | 2.57×10-8 | 1.82×10-4 | 2.83×10-9 | |
0.8 | 6.87×10-4 | 3.42×10-8 | 3.02×10-4 | 4.35×10-9 | |
1 | 9.61×10-4 | 4.26×10-8 | 4.46×10-4 | 6.08×10-9 |
t | α | Solution by q-HASTM | Solution by methods in [1,5] |
0.2 | 1.5 | 12.397642 | 12.393637 |
2 | 12.111344 | 12.111344 | |
0.4 | 1.5 | 19.483525 | 19.444664 |
2 | 18.560015 | 18.560015 | |
0.6 | 1.5 | 27.469770 | 27.306178 |
2 | 25.526848 | 25.526871 | |
0.8 | 1.5 | 36.723342 | 36.253258 |
2 | 33.212233 | 33.212664 | |
1 | 1.5 | 47.895372 | 46.896428 |
2 | 41.851853 | 41.856135 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 1.08×10-3 | 1.78×10-5 | 3.32×10-4 | 1.23×10-6 | |
0.4 | 1.52×10-3 | 1.75×10-5 | 5.58×10-4 | 1.74×10-6 | |
0.6 | 1.86×10-3 | 1.73×10-5 | 7.57×10-4 | 2.13×10-6 | |
0.8 | 2.15×10-3 | 1.74×10-5 | 9.40×10-4 | 2.46×10-6 | |
1 | 2.40×10-3 | 1.72×10-5 | 1.11×10-3 | 2.75×10-6 |
r | t | Solution by q-HASTM | Solution by operational Matrix method [3] | Solution by method in [4] |
0.2 | 0.49301 | 0.4900 | 0.4925 | |
0.4 | 0.69733 | 0.6910 | 0.6965 | |
0.6 | 0.85405 | 0.8318 | 0.8531 | |
0.8 | 0.98618 | 0.9727 | 0.9850 | |
1 | 1.10259 | 1.1763 | 1.1013 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 2.02×10-3 | 7.80×10-5 | 6.10×10-4 | 5.24×10-6 | |
0.4 | 2×10-3 | 5.36×10-5 | 7.25×10-4 | 5.22×10-6 | |
0.6 | 1.99×10-3 | 4.33×10-5 | 8.02×10-4 | 5.21×10-6 | |
0.8 | 1.99×10-3 | 3.72×10-5 | 8.61×10-4 | 5.20×10-6 | |
1 | 1.98×10-3 | 3.31×10-5 | 9.11×10-3 | 5.21×10-6 |
In this investigation, q-HASTM is effectively applied to solve the time fractional vibration equation. The outcomes show that the derived results are trustworthy. The simulations shown confirm great accuracy of gained results. It is found that this scheme is capable of diminishing the calculation size. The q-HASTM has parameters n,ℏ that manage solution convergence. It is exciting to observe that the q-HAM works efficiently when coupled with Sumudu transform due to its ‘unity’ feature. Also, the nonlinear term can easily be handled via Sumudu transform. Hence, it is concluded that this scheme is accurate, systematic, logical, easy to use and attractive. It can be applied to study a wide variety of arbitrary order models of physical, biological, medical and social importance.
Authors declare that there is no conflict of interest in this paper.
[1] |
H. M. Srivastava, D. Kumar, J. Singh, An efficient analytical technique for fractional model of vibration equation, Appl. Math. Model., 45 (2017), 192-204. doi: 10.1016/j.apm.2016.12.008
![]() |
[2] | H. Singh, Approximate solution of fractional vibration equation using Jacobi polynomials, Appl. Math. Comput., 317 (2018), 85-100. |
[3] |
H. Singh, H. M. Srivastava, D. Kumar, A reliable numerical algorithm for the fractional vibration equation, Chaos Soliton. Fract., 103 (2017), 131-138. doi: 10.1016/j.chaos.2017.05.042
![]() |
[4] |
S. Das, A numerical solution of the vibration equation using modified decomposition method, J. Sound Vib., 320 (2009), 576-583. doi: 10.1016/j.jsv.2008.08.029
![]() |
[5] |
S. Das, P. K. Gupta, Application of homotopy perturbation method and homotopy analysis method for fractional vibration equation, Int. J. Comput. Math., 88 (2011), 430-441. doi: 10.1080/00207160903474214
![]() |
[6] | S. T. Mohyud-Din, A. Yildirim, An algorithm for solving the fractional vibration equation, Comput. Math. Model., 23 (2012), 228-237. |
[7] |
H. M. Baskonus, T. Mekkaoui, Z. Hammouch, et al. Active control of a chaotic fractional order economic system, Entropy, 17 (2015), 5771-5783. doi: 10.3390/e17085771
![]() |
[8] | H. M. Baskonus, H. Bulut, On the numerical solutions of some fractional ordinary differential equations by fractional Adams-Bashforth-Moulton method, Open Math., 13 (2015), 547-556. |
[9] | H. M. Baskonus, G. Yel, H. Bulut, Novel wave surfaces to the fractional Zakharov- Kuznetsov Benjamin-Bona-Mahony equation, AIP Conference Proceedings, 1863 (2017), 560084. |
[10] | C. Ravichandran, K. Jothimani, H. M. Baskonus, et al. New results on nondensely characterized integrodifferential equations with fractional order, Eur. Phys. J. Plus, 133 (2018), 109. |
[11] |
R. Panda, M. Dash, Fractional generalized splines and signal processing, Signal Process, 86 (2006), 2340-2350. doi: 10.1016/j.sigpro.2005.10.017
![]() |
[12] |
R. L. Bagley, P. J. Torvik, A theoretical basis for the application of fractional calculus to Viscoelasticity, J. Rheol., 27(1983), 201-210. doi: 10.1122/1.549724
![]() |
[13] | A. Prakash, M. Goyal, S. Gupta, q-homotopy analysis method for fractional Bloch model arising in nuclear magnetic resonance via the Laplace transform, Ind. J. Phys., (2019), http://doi.org/10.1007/s12648-019-01487-7. |
[14] | A. Prakash, M. Goyal, S. Gupta, Numerical simulation of space-fractional Helmholtz equation arising in Seismic wave propagation, imaging and inversion, Pramana, 93 (2019), 28. |
[15] | M. Goyal, H. M. Baskonus, A. Prakash, An efficient technique for a time fractional model of lassa hemorrhagic fever spreading in pregnant women, Eur. Phys. J. Plus, 134 (2019), 482. |
[16] | S. Das, Solution of fractional vibration equation by the variational iteration method and modified decomposition method, Int. J. Nonlin. Sci. Num., 9 (2008), 361-366. |
[17] | S. J. Liao, On the homotopy analysis method for non-linear problems, Appl. Math. Comput., 147 (2004), 499-513. |
[18] | M. A. El-Tawil, S. N. Huseen, The q-homotopy analysis method (q-HAM), Int. J. Appl. Math. Mech., 8 (2012), 51-75. |
[19] |
M. A. El-Tawil, S. N. Huseen, On convergence of the q-homotopy analysis method, Int. J. Contemp. Math. Sci., 8 (2013), 481-497. doi: 10.12988/ijcms.2013.13048
![]() |
[20] | A. Prakash, P. Veeresha, D. G. Prakasha, et al. A homotopy technique for a fractional order multi-dimensional telegraph equation via the Laplace Transform, Eur. Phys. J. Plus, 134 (2019), 19. |
[21] | A. Prakash, P. Veeresha, D. G. Prakasha, et al. A new efficient technique for solving fractional coupled Navier-Stokes equations using q-homotopy analysis transform method, Pramana, 93 (2019), 6. |
[22] |
A. Prakash, M. Goyal, S. Gupta, Fractional variational iteration method for solving time-fractional Newell-Whitehead-Segel equation, Nonlinear Eng., 8 (2019), 164-171. doi: 10.1515/nleng-2018-0001
![]() |
[23] | M. Goyal, A. Prakash, S. Gupta, Numerical simulation for time-fractional nonlinear coupled dynamical model of romantic and interpersonal relationships, Pramana, 92 (2019), 82. |
[24] | A. Prakash, M. Goyal, S. Gupta, A reliable algorithm for fractional Bloch model arising in magnetic resonance imaging, Pramana, 92 (2019), 18. |
[25] |
A. M. S. Mahdy, A. S. Mohamed, A. A. H. Mtawa, Implementation of the Homotopy perturbation Sumudu transform method for solving Klein-Gordon equation, Appl. Math., 6 (2015), 617-628. doi: 10.4236/am.2015.63056
![]() |
[26] | A. A. Elbeleze, A. Kilicman, B. M. Taib, Homotopy perturbation method for fractional Black-Scholes European option pricing equations using Sumudu transform, Math. Probl. Eng., 2013 (2013), 524852. |
[27] | G. K. Watugala, Sumudu transform: A new integral transform to solve differential equations and control engineering problems, Integr. Educ., 24 (1993), 35-43. |
[28] | G. K. Watugala, The Sumudu transform for functions of two variables, Math. Eng. Ind., 8 (2002), 293-302. |
[29] |
M. A. Asiru, Further properties of the Sumudu transform and its applications, Int. J. Math. Educ. Sci. Tech., 33 (2002), 441-449. doi: 10.1080/002073902760047940
![]() |
[30] |
S. Weerakoon, Applications of Sumudu transform to partial differential equations, Int. J. Math. Educ. Sci. Tech., 25 (1994), 277-283. doi: 10.1080/0020739940250214
![]() |
[31] | S. Weerakoon, Complex inversion formula for Sumudu transforms, Int. J. Math. Educ. Sci. Tech., 29 (1998), 618-621. |
[32] | J. Singh, D. Kumar, D. Baleanu, et al. An efficient numerical algorithm for the fractional Drinfeld-Sokolov-Wilson equation, Appl. Math. Comput., 335 (2018), 12-24. |
[33] | I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, 1999. |
[34] | M. Caputo, Elasticita e Dissipazione, Zani-Chelli, 1969. |
[35] | K. Diethelm, The Analysis of Fractional Differential Equations, Springer-Verlag, 2004. |
[36] | A. Prakash, M. Kumar, D. Baleanu, A new iterative technique for a fractional model of nonlinear Zakharov-Kuznetsov equations via Sumudu transform, Appl. Math. Comput., 334 (2004), 30-40. |
[37] | J. Choi, D. Kumar, J. Singh, et al. Analytical techniques for system of time fractional nonlinear differential equations, J. Korean Math. Soc., 54 (2017), 1209-1229. |
1. | Manish Goyal, Haci Mehmet Baskonus, Amit Prakash, Regarding new positive, bounded and convergent numerical solution of nonlinear time fractional HIV/AIDS transmission model, 2020, 139, 09600779, 110096, 10.1016/j.chaos.2020.110096 | |
2. | Majid Bagheri, Ali Khani, Soheil Salahshour, Analytical Method for Solving the Fractional Order Generalized KdV Equation by a Beta-Fractional Derivative, 2020, 2020, 1687-9139, 1, 10.1155/2020/8819183 | |
3. | Mehvish Fazal Ur Rehman, Yongyi Gu, Wenjun Yuan, Soheil Salahshour, Exact Analytical Solutions of Nonlinear Fractional Liouville Equation by Extended Complex Method, 2020, 2020, 1687-9139, 1, 10.1155/2020/8815363 | |
4. | Dumitru Baleanu, Amin Jajarmi, Hakimeh Mohammadi, Shahram Rezapour, A new study on the mathematical modelling of human liver with Caputo–Fabrizio fractional derivative, 2020, 134, 09600779, 109705, 10.1016/j.chaos.2020.109705 | |
5. | Amit Prakash, Vijay Verma, Two efficient computational technique for fractional nonlinear Hirota–Satsuma coupled KdV equations, 2020, 38, 0264-4401, 791, 10.1108/EC-02-2020-0091 | |
6. | Harendra Singh, Numerical simulation for fractional delay differential equations, 2020, 2195-268X, 10.1007/s40435-020-00671-6 | |
7. | Wei Gao, Hatıra Günerhan, Haci Mehmet Baskonus, Analytical and approximate solutions of an epidemic system of HIV/AIDS transmission, 2020, 59, 11100168, 3197, 10.1016/j.aej.2020.07.043 | |
8. | D.G. Prakasha, P. Veeresha, Analysis of Lakes pollution model with Mittag-Leffler kernel, 2020, 5, 24680133, 310, 10.1016/j.joes.2020.01.004 | |
9. | Amit Prakash, Ajay Kumar, Haci Mehmet Baskonus, Ashok Kumar, Numerical analysis of nonlinear fractional Klein–Fock–Gordon equation arising in quantum field theory via Caputo–Fabrizio fractional operator, 2021, 2008-1359, 10.1007/s40096-020-00365-2 | |
10. | Mutaz Mohammad, Carlo Cattani, A collocation method via the quasi-affine biorthogonal systems for solving weakly singular type of Volterra-Fredholm integral equations, 2020, 59, 11100168, 2181, 10.1016/j.aej.2020.01.046 | |
11. | V. Padmavathi, A. Prakash, K. Alagesan, N. Magesh, Analysis and numerical simulation of novel coronavirus (COVID‐19) model with Mittag‐Leffler Kernel, 2021, 44, 0170-4214, 1863, 10.1002/mma.6886 | |
12. | Manish Goyal, Vinod Kumar Bhardwaj, Amit Prakash, Investigating new positive, bounded, and convergent numerical solution for the nonlinear time‐dependent breast cancer dynamic competition model, 2021, 44, 0170-4214, 4636, 10.1002/mma.7057 | |
13. | Hülya Durur, Asif Yokuş, Exact solutions of (2 + 1)-Ablowitz-Kaup-Newell-Segur equation, 2020, 0, 2444-8656, 10.2478/amns.2020.2.00074 | |
14. | Uttam Ghosh, Susmita Sarkar, Bapin Mondal, Study of Stability and Bifurcation of Three Species Food Chain Model with Non-monotone Functional Response, 2021, 7, 2349-5103, 10.1007/s40819-021-01017-2 | |
15. | Shao-Wen Yao, Muhammad Farman, Maryam Amin, Mustafa Inc, Ali Akgül, Aqeel Ahmad, Fractional order COVID-19 model with transmission rout infected through environment, 2022, 7, 2473-6988, 5156, 10.3934/math.2022288 | |
16. | Amit Kumar Saraswat, Manish Goyal, Numerical simulation of time-dependent influenza model with Atangana–Baleanu non-integer order derivative in Liouville–Caputo sense, 2022, 96, 0973-7111, 10.1007/s12043-022-02335-w | |
17. | V. Padmavathi, N. Magesh, K. Alagesan, M. Ijaz Khan, Samia Elattar, Mamdooh Alwetaishi, Ahmed M. Galal, Numerical Modeling and Symmetry Analysis of a Pine Wilt Disease Model Using the Mittag–Leffler Kernel, 2022, 14, 2073-8994, 1067, 10.3390/sym14051067 | |
18. | Amit Prakash, Numerical simulation of SIR childhood diseases model with fractional Adams–Bashforth method, 2022, 0170-4214, 10.1002/mma.8785 | |
19. | Amit Prakash, Hardish Kaur, 2022, Chapter 38, 978-981-16-7663-5, 465, 10.1007/978-981-16-7664-2_38 | |
20. | Manish Goyal, Amit Prakash, Shivangi Gupta, An Efficient Perturbation Sumudu Transform Technique for the Time-Fractional Vibration Equation with a Memory Dependent Fractional Derivative in Liouville–Caputo Sense, 2021, 7, 2349-5103, 10.1007/s40819-021-01068-5 | |
21. | Amit Prakash, Hardish Kaur, Manish Goyal, 2022, Chapter 31, 978-981-16-7663-5, 383, 10.1007/978-981-16-7664-2_31 | |
22. | Manish Goyal, Haci Mehmet Baskonus, A Reliable Solution of Arbitrary Order Nonlinear Hunter–Saxton Equation with Time Dependent Derivative in Liouville–Caputo Sense, 2021, 7, 2349-5103, 10.1007/s40819-021-01065-8 | |
23. | Manish Goyal, Amit Prakash, Dumitru Baleanu, An efficient hybrid computational technique for the time dependent Lane-Emden equation of arbitrary order, 2022, 7, 24680133, 131, 10.1016/j.joes.2021.07.004 | |
24. | Rania Saadeh, Aliaa Burqan, Ahmad El-Ajou, Reliable solutions to fractional Lane-Emden equations via Laplace transform and residual error function, 2022, 61, 11100168, 10551, 10.1016/j.aej.2022.04.004 | |
25. | Viswanathan Padmavathi, Kandaswami Alagesan, Awatif Alhowaity, Muhammad Ijaz Khan, Haneen Hamam, Manivelu Angayarkanni, Vediyappan Govindan, Analysis and numerical simulation of fractional order pine wilt disease model, 2022, 36, 0217-9792, 10.1142/S0217979222502046 | |
26. | Shivangi Gupta, Manish Goyal, Amit Prakash, A Hybrid Computational Scheme with Convergence Analysis for the Dependent Rosenau-Hyman Equation of Arbitrary Order Via Caputo-Fabrizio Operator, 2021, 7, 2349-5103, 10.1007/s40819-021-01182-4 | |
27. | Hamadou Halidou, Alphonse Houwe, Souleymanou Abbagari, Mustafa Inc, Serge Y. Doka, Thomas Bouetou Bouetou, Optical and W-shaped bright solitons of the conformable derivative nonlinear differential equation, 2021, 20, 1569-8025, 1739, 10.1007/s10825-021-01758-9 | |
28. | Manish Goyal, Amit Kumar Saraswat, Amit Prakash, Numerical analysis of fractional coronavirus model with Atangana-Baleanu derivative in Liouville-Caputo sense, 2023, 97, 0973-1458, 147, 10.1007/s12648-022-02409-w | |
29. | Kazem Nouri, Hassan Ranjbar, Dumitru Baleanu, Leila Torkzadeh, Investigation on Ginzburg-Landau equation via a tested approach to benchmark stochastic Davis-Skodje system, 2021, 60, 11100168, 5521, 10.1016/j.aej.2021.04.040 | |
30. | Vinod Kumar Bhardwaj, Manish Goyal, A Reliable Solution of Nonlinear Time Dependent Fractional Model of Ebola Virus Disease with Arbitrary Order Derivative in Liouville–Caputo Sense, 2021, 7, 2349-5103, 10.1007/s40819-021-01200-5 | |
31. | A. DLAMINI, EMILE F. DOUNGMO GOUFO, M. KHUMALO, CHAOTIC BEHAVIOR OF MODIFIED STRETCH–TWIST–FOLD FLOW UNDER FRACTAL-FRACTIONAL DERIVATIVES, 2022, 30, 0218-348X, 10.1142/S0218348X22402071 | |
32. | Amit Prakash, Hardish Kaur, Numerical simulation of coupled fractional‐order Whitham‐Broer‐Kaup equations arising in shallow water with Atangana‐Baleanu derivative, 2022, 0170-4214, 10.1002/mma.8238 | |
33. | Ajay Kumar, Haci Mehmet Baskonus, Amit Prakash, 2023, Chapter 14, 978-3-031-29958-2, 232, 10.1007/978-3-031-29959-9_14 | |
34. | MUHAMMAD AZEEM, MUHAMMAD FARMAN, MARWAN ABUKHALED, KOTTAKKARAN SOOPPY NISAR, ALI AKGÜL, EPIDEMIOLOGICAL ANALYSIS OF HUMAN LIVER MODEL WITH FRACTIONAL OPERATOR, 2023, 31, 0218-348X, 10.1142/S0218348X23400479 | |
35. | Amit Prakash, Vijay Verma, Dumitru Baleanu, Two Novel Methods for Fractional Nonlinear Whitham–Broer–Kaup Equations Arising in Shallow Water, 2023, 9, 2349-5103, 10.1007/s40819-023-01497-4 | |
36. | Viswanathan Padmavathi, Kandaswami Alagesan, Samad Noeiaghdam, Unai Fernandez-Gamiz, Manivelu Angayarkanni, Vediyappan Govindan, Tobacco smoking model containing snuffing class, 2023, 9, 24058440, e20792, 10.1016/j.heliyon.2023.e20792 | |
37. | A. Dlamini, E. F. Doungmo Goufo, M. Khumalo, Strongly perturbed bondorbital attractors for generalized systems, 2025, 35, 1054-1500, 10.1063/5.0249237 |
r | t | α | Solution by q-HASTM | Solution by operational Matrix method [3] | Solution by Methods in [1,4,5,6,16] |
0.2 | 1.5 | 0.0467182 | 0.0469 | 0.0467 | |
2 | 0.0448067 | 0.0448 | 0.0448 | ||
0.4 | 1.5 | 0.183688 | 0.1838 | 0.1837 | |
2 | 0.179227 | 0.1792 | 0.1792 | ||
0.6 | 1.5 | 0.410124 | 0.4100 | 0.4101 | |
2 | 0.40326 | 0.4032 | 0.4032 | ||
0.8 | 1.5 | 0.725746 | 0.7256 | 0.7257 | |
2 | 0.716907 | 0.7169 | 0.7169 | ||
1 | 1.5 | 1.13039 | 1.1310 | 1.1303 | |
2 | 1.12017 | 1.1203 | 1.1202 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 8.59×10-5 | 8.77×10-9 | 2.66×10-5 | 5.46×10-10 | |
0.4 | 2.43×10-4 | 1.72×10-8 | 8.97×10-5 | 1.54×10-9 | |
0.6 | 4.46×10-4 | 2.57×10-8 | 1.82×10-4 | 2.83×10-9 | |
0.8 | 6.87×10-4 | 3.42×10-8 | 3.02×10-4 | 4.35×10-9 | |
1 | 9.61×10-4 | 4.26×10-8 | 4.46×10-4 | 6.08×10-9 |
t | α | Solution by q-HASTM | Solution by methods in [1,5] |
0.2 | 1.5 | 12.397642 | 12.393637 |
2 | 12.111344 | 12.111344 | |
0.4 | 1.5 | 19.483525 | 19.444664 |
2 | 18.560015 | 18.560015 | |
0.6 | 1.5 | 27.469770 | 27.306178 |
2 | 25.526848 | 25.526871 | |
0.8 | 1.5 | 36.723342 | 36.253258 |
2 | 33.212233 | 33.212664 | |
1 | 1.5 | 47.895372 | 46.896428 |
2 | 41.851853 | 41.856135 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 1.08×10-3 | 1.78×10-5 | 3.32×10-4 | 1.23×10-6 | |
0.4 | 1.52×10-3 | 1.75×10-5 | 5.58×10-4 | 1.74×10-6 | |
0.6 | 1.86×10-3 | 1.73×10-5 | 7.57×10-4 | 2.13×10-6 | |
0.8 | 2.15×10-3 | 1.74×10-5 | 9.40×10-4 | 2.46×10-6 | |
1 | 2.40×10-3 | 1.72×10-5 | 1.11×10-3 | 2.75×10-6 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 2.02×10-3 | 7.80×10-5 | 6.10×10-4 | 5.24×10-6 | |
0.4 | 2×10-3 | 5.36×10-5 | 7.25×10-4 | 5.22×10-6 | |
0.6 | 1.99×10-3 | 4.33×10-5 | 8.02×10-4 | 5.21×10-6 | |
0.8 | 1.99×10-3 | 3.72×10-5 | 8.61×10-4 | 5.20×10-6 | |
1 | 1.98×10-3 | 3.31×10-5 | 9.11×10-3 | 5.21×10-6 |
r | t | α | Solution by q-HASTM | Solution by operational Matrix method [3] | Solution by Methods in [1,4,5,6,16] |
0.2 | 1.5 | 0.0467182 | 0.0469 | 0.0467 | |
2 | 0.0448067 | 0.0448 | 0.0448 | ||
0.4 | 1.5 | 0.183688 | 0.1838 | 0.1837 | |
2 | 0.179227 | 0.1792 | 0.1792 | ||
0.6 | 1.5 | 0.410124 | 0.4100 | 0.4101 | |
2 | 0.40326 | 0.4032 | 0.4032 | ||
0.8 | 1.5 | 0.725746 | 0.7256 | 0.7257 | |
2 | 0.716907 | 0.7169 | 0.7169 | ||
1 | 1.5 | 1.13039 | 1.1310 | 1.1303 | |
2 | 1.12017 | 1.1203 | 1.1202 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 8.59×10-5 | 8.77×10-9 | 2.66×10-5 | 5.46×10-10 | |
0.4 | 2.43×10-4 | 1.72×10-8 | 8.97×10-5 | 1.54×10-9 | |
0.6 | 4.46×10-4 | 2.57×10-8 | 1.82×10-4 | 2.83×10-9 | |
0.8 | 6.87×10-4 | 3.42×10-8 | 3.02×10-4 | 4.35×10-9 | |
1 | 9.61×10-4 | 4.26×10-8 | 4.46×10-4 | 6.08×10-9 |
t | α | Solution by q-HASTM | Solution by methods in [1,5] |
0.2 | 1.5 | 12.397642 | 12.393637 |
2 | 12.111344 | 12.111344 | |
0.4 | 1.5 | 19.483525 | 19.444664 |
2 | 18.560015 | 18.560015 | |
0.6 | 1.5 | 27.469770 | 27.306178 |
2 | 25.526848 | 25.526871 | |
0.8 | 1.5 | 36.723342 | 36.253258 |
2 | 33.212233 | 33.212664 | |
1 | 1.5 | 47.895372 | 46.896428 |
2 | 41.851853 | 41.856135 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 1.08×10-3 | 1.78×10-5 | 3.32×10-4 | 1.23×10-6 | |
0.4 | 1.52×10-3 | 1.75×10-5 | 5.58×10-4 | 1.74×10-6 | |
0.6 | 1.86×10-3 | 1.73×10-5 | 7.57×10-4 | 2.13×10-6 | |
0.8 | 2.15×10-3 | 1.74×10-5 | 9.40×10-4 | 2.46×10-6 | |
1 | 2.40×10-3 | 1.72×10-5 | 1.11×10-3 | 2.75×10-6 |
r | t | Solution by q-HASTM | Solution by operational Matrix method [3] | Solution by method in [4] |
0.2 | 0.49301 | 0.4900 | 0.4925 | |
0.4 | 0.69733 | 0.6910 | 0.6965 | |
0.6 | 0.85405 | 0.8318 | 0.8531 | |
0.8 | 0.98618 | 0.9727 | 0.9850 | |
1 | 1.10259 | 1.1763 | 1.1013 |
r | t | Solution by q-HASTM | |||
α=1.5 | α=1.75 | ||||
|w2-w1| | |w3-w2| | |w2-w1| | |w3-w2| | ||
0.2 | 2.02×10-3 | 7.80×10-5 | 6.10×10-4 | 5.24×10-6 | |
0.4 | 2×10-3 | 5.36×10-5 | 7.25×10-4 | 5.22×10-6 | |
0.6 | 1.99×10-3 | 4.33×10-5 | 8.02×10-4 | 5.21×10-6 | |
0.8 | 1.99×10-3 | 3.72×10-5 | 8.61×10-4 | 5.20×10-6 | |
1 | 1.98×10-3 | 3.31×10-5 | 9.11×10-3 | 5.21×10-6 |