
Citation: Kyohei Fujita, Zhonggang Feng, Daisuke Sato, Tadashi Kosawada, Takao Nakamura, Yasuyuki Shiraishi, Mitsuo Umezu. Modulation of the mechanical properties of ventricular extracellular matrix hydrogels with a carbodiimide crosslinker and investigation of their cellular compatibility[J]. AIMS Materials Science, 2018, 5(1): 54-74. doi: 10.3934/matersci.2018.1.54
[1] | Bashir Sajo Mienda, Mohd Shahir Shamsir . Model-driven in Silico glpC Gene Knockout Predicts Increased Succinate Production from Glycerol in Escherichia Coli. AIMS Bioengineering, 2015, 2(2): 40-48. doi: 10.3934/bioeng.2015.2.40 |
[2] | Wiebke Wesseling . Beneficial biofilms in marine aquaculture? Linking points of biofilm formation mechanisms in Pseudomonas aeruginosa and Pseudoalteromonas species. AIMS Bioengineering, 2015, 2(3): 104-125. doi: 10.3934/bioeng.2015.3.104 |
[3] | Bashir Sajo Mienda, Faezah Mohd Salleh . Bio-succinic acid production: Escherichia coli strains design from genome-scale perspectives. AIMS Bioengineering, 2017, 4(4): 418-430. doi: 10.3934/bioeng.2017.4.418 |
[4] | Mingyong Xiong, Ping Yu, Jingyu Wang, Kechun Zhang . Improving Engineered Escherichia coli strains for High-level Biosynthesis of Isobutyrate. AIMS Bioengineering, 2015, 2(2): 60-74. doi: 10.3934/bioeng.2015.2.60 |
[5] | Sezer Okay, Mehmet Sezgin . Transgenic plants for the production of immunogenic proteins. AIMS Bioengineering, 2018, 5(3): 151-161. doi: 10.3934/bioeng.2018.3.151 |
[6] | Dyoni M. de Oliveira, Victor Hugo Salvador, Thatiane R. Mota, Aline Finger-Teixeira, Rodrigo F. de Almeida, Douglas A. A. Paixão, Amanda P. De Souza, Marcos S. Buckeridge, Rogério Marchiosi, Osvaldo Ferrarese-Filho, Fabio M. Squina, Wanderley D. dos Santos . Feruloyl esterase from Aspergillus clavatus improves xylan hydrolysis of sugarcane bagasse. AIMS Bioengineering, 2017, 4(1): 1-11. doi: 10.3934/bioeng.2017.1.1 |
[7] | Hadi Nazem-Bokaee, Ryan S. Senger . ToMI-FBA: A genome-scale metabolic flux based algorithm to select optimum hosts and media formulations for expressing pathways of interest. AIMS Bioengineering, 2015, 2(4): 335-374. doi: 10.3934/bioeng.2015.4.335 |
[8] | Liwei Chen, Jaslyn Lee, Wei Ning Chen . The use of metabolic engineering to produce fatty acid-derived biofuel and chemicals in Saccharomyces cerevisiae: a review. AIMS Bioengineering, 2016, 3(4): 468-492. doi: 10.3934/bioeng.2016.4.468 |
[9] | Chiara Ceresa, Maurizio Rinaldi, Letizia Fracchia . Synergistic activity of antifungal drugs and lipopeptide AC7 against Candida albicans biofilm on silicone. AIMS Bioengineering, 2017, 4(2): 318-334. doi: 10.3934/bioeng.2017.2.318 |
[10] | Bashir Sajo Mienda . Escherichia coli Genome-scale metabolic models could guide construction of proof-of-principle strains. AIMS Bioengineering, 2018, 5(2): 103-105. doi: 10.3934/bioeng.2018.2.103 |
In the past few years, fractal calculus and fractional calculus has been one of the most rapidly growing areas of mathematical analysis, which are used to model various problems of real life. Fractal calculus is considered as a fruitful field of research in science and technology and is extensively used to elucidate the phenomena of hierarchical or porous media [1]. Owing to the simplicity and effectiveness of fractal derivative it has defined an alternative approach of fractional derivatives, to explain some of the most fundamental theories with considerable ease and elegance [2,3]. The fractional derivative and fractional integral operators are the valuable tools which are used in the modeling of various physical phenomena of science and engineering. Anomalous dynamics of numerous complex nonlinear systems are elegantly modeled by the assistance of fractional differential equations, which are acknowledged as the generalization of the classical differential equations of integer order [4,5,6].
The fractional order Ricatti equation (RE) accommodates an extremely significant class of nonlinear differential equations which includes: The nth order linear homogenous ODE, one-dimensional Schrödinger equation, wave solution of nonlinear partial differential equation (PDE), etc [7]. These REs has substantial importance in classical, as well as, modern science and engineering problems because of their multifarious applications in various fields. For instance, optimal control, variational calculus, random processes, quantum mechanics, thermodynamics, robust stabilization, stochastic realization theory and diffusion problems [8,9,10,11].
The mathematical formulation for the nonlinear Ricatti differential equations with the fractional derivative defined in the Caputo sense is read as,
CDα0x(t)=p(t)x2(t)+q(t)x(t)+r(t),0⩽t⩽T | (1) |
with the initial condition defined as,
x(t0)=η, | (2) |
where, α is the order of equation, α>0 and α∈ℜ, with T,t0,η∈ℜ; p(t),q(t),r(t) are the real continuous functions and x(t)is the solution of the equation. The behavior of the aforementioned Eq. (1) depends on the parameter α, which can be varied to analyze the dynamical behavior of these equations. In order to study the nonlinear Eq. (1) various analytical and numerical techniques developed in the past were applied by many researchers and scientists to compute the solutions of these equations among which some of them are cited here: The homotopy perturbation method (HPM) [12], the fractional-order Legendre operational matrix method [13], the iterative reproducing kernel Hilbert space method [14], the optimal homotopy asymptotic method [15], the modified Laplace Adomian decomposition method [16], the B-spline operational method [17] and the Haar wavelet collocation method [18]. The HPM introduced by He has been successfully used to solve numerous linear and nonlinear problems of real life. Since, its development it has been modified by He and many other scientists to solve various ordinary and partial differential equations of integer and non-integer order [19,20,21].
The demand of global optimization technique is increasing day by day, which are utilized in the assessment of numerous nonlinear and multimodal problems of real life. The deterministic algorithm and the stochastic algorithm are the two types of optimization algorithm that are found in the literature [22,23]. The deterministic algorithm are often gradient-based whereas, the stochastic algorithm are further subcategorized as heuristic or metaheuristic algorithm. In recent trends, the popularity and demand of the nature inspired algorithms have increased extensively. Nature inspired metaheuristic algorithms which efficiently deals with the nonlinear optimization problems includes the genetic algorithm (GA) [24], simulated annealing (SA) algorithm [25], differential evolution (DE) algorithm [26], the ant colony optimization (ACO) algorithm [27], particle swarm optimization (PSO) algorithms [28], the shark smell algorithm [29] and the most powerful and demanding firefly (FA) algorithm [30,31,32].
The FA which mimics the flashing pattern and behavior of the fireflies was first developed by Yang in 2008 and since then it has been modified by various scientist and researcher to solve different types of challenging optimization problems. Some of the flashing characteristic of these unisex fireflies are idealized to develop the FA. The three idealized rule used by the firefly algorithm are as follows:
● A firefly is attracted to the other fireflies regardless of their gender.
● Attractiveness of fireflies is directly proportional to their brightness. The attractiveness and brightness of fireflies both increases as the distance between them deceases.
● Brightness of a firefly is determined by the objective function.
The literature of this modern, self-adaptive, highly efficient and truly intelligent algorithm has expanded dramatically [30,31,32,33,34].
The fundamental aim of this paper is to provide a numerical technique for the assessment of nonlinear fractional differential model defined as (1) and (2). The concept of classical homotopy perturbation method is merged with the modern metaheuristic optimization technique for the development of an expedite homptopy perturbation method (EHPM). The developed method-EHPM transforms the fractional model into a system of algebraic equations leading to a fitness function determined by a particular fragment of the weighted series solution, which are trained by the using the powerful and reliable optimization technique-FA. The optimal values achieved by the FA are utilized to attain the accurate, convergent and reliable solutions. Comparative study is conducted by the comparing the EHPM computed results with the available exact solution and the solutions obtained by the modified homotopy perturbation (MHPM) [35], the residual power series method (RPSM) [36] and the Adam bashforth method (ABMA) [37]. Furthermore, the accuracy and competency of the EHPM is also ratified by computing results by the proposed design methodology in combination with accelerated particle swarm optimization (APSO), i.e., the fitness function determined by EHPM is optimized by using APSO. Various error measures are also carried out to validate the correctness and accuracy of the suggested scheme.
Some basic definitions of the fractional calculus which are going to be utilized in the further discussion are stated below:
Definition 2.1. Let f(x) be a differential function with β∈(0,1],then the Caputo order fractional derivative CDβ0f(x), is defined as [36],
CDβ0f(x)={1Γ(1−β)x∫0(x−ϖ)−βf′(ϖ)dϖ,0<β<1,df(x)dx,β=1, | (3) |
with CDβ0ω=0, for some constant ω∈R.
Definition 2.2. For any function f(x) with x:[0,∞)→ℜ, the Riemann fractional integral of order β is given as [35],
RIβ0f(x)=1Γ(β)x∫0(x−ϖ)β−1f(ϖ)dϖ, | (4) |
where β>0 and β∈ℜ.
This section comprehensively describes the procedure to determine the approximate solutions of the nonlinear fractional differential equations using the classical idea of HPM with the modern optimizing tool. A brief review of the learning solver FA, which is utilized for the development of the presented scheme (EHPM) is also demonstrated here.
To exemplify the basic ideas of the proposed scheme consider the nonlinear FDE of the form,
ψ|(x,t,g)=0,Ω=[0,T]∈ℜ | (5) |
with the initial condition defined as
x(t0)=η,η,t0∈ℜ, | (6) |
where g and xare function of t and Ω is the boundary of the domain. Then ψ can be expressed as
L(x,t)+N(x,t)+A(g(t))=0, | (7) |
Where L is the linear operator, N is the nonlinear operator and A is the known analytical function. By the homotopy procedure we rewrite Eq. (7) as
H(X(t),δ)≡L(X(t),t)−L(x0(t),t)+δL(x0(t),t)+δ(N(X(t),t)+A(g(t)))=0 | (8) |
where x0(t) is an initial approximation, and δ is an embedding parameter defined in some closed interval [0,1]. Let the homotopy solution of equation Eq. (8) can be written as
X(t)=∞∑k=0δkXk(t) | (9) |
where Xk(t) are attained by the kth-order homotopy derivative defined as
Xk(t)=1k!∂kH∂δk|δ=0 | (10) |
Now, we let the initial approximation of Eq. (8) of the form
x0(t)=n∑i=0tiαξiΓ(iα+1) | (11) |
where ξ0,ξ1,ξ2,…,ξn are the unknown weights to be determined. By utilizing the Eq. (11), Eq. (9) and the basic ideas of MHPM [35] with X1(t)=0 and L−1 as the inverse operator of L we obtain an algebraic system expressed as
X0=L−1(n∑i=0tiαξiΓ(iα+1)) | (12) |
X1(t)=L−1(L(x0(t),t)+N(X(t),t)+A(g(t))) | (13) |
and X2(t)=X3(t)=…=0. Setting X1(t)=0 by MHPM, automatically leads the solution of Eq. (8) written as:
x(t)=X(t)=X0=L−1(n∑i=0tiαξiΓ(iα+1)) | (14) |
The aforementioned approximate solution (14) comprises of the unknown weights ξ0,ξ1,ξ2,…,ξn which are determined by constructing a fitness function given as
E(ξi)=minm∑j=1|X1(tj)| | (15) |
The weights in Eq. (15) are learned by using a modern and powerful optimization technique FA. The graphical abstract of the above presented scheme EHPM is portrayed in Figure 1.
Lemma 3.1. In problem (5) if x(t)=X0 is the solutions of Eq. (9), which comprises of the unknown weights ξi;i=0,1,2,...,n, then the minimization problem will be simplified as:
E(ξi)=minm∑j=1|(tj,X0(tj,ξi),DαtX0(tj,ξi))| | (16) |
Theorem 3.2. If the solution of Eq. (8) satisfies X1(t)=0, then Eq. (10) results X2(t)=X3(t)=…=0 and x(t)=X0 as the solution of Eq. (8).
Here we mention that if g(r(t)) and X0 are analytic at t=t0, then their power series is defined as
X0=∞∑n=0ξn(t−t0)nandg(t)=∞∑n=0ξn∗(t−t0)n | (17) |
where ξ0∗,ξ1∗,ξ2∗,…are known coefficients and ξ0,ξ1,ξ2,…are unknown ones which are to be determined. In order to obtain ξ0,ξ1,ξ2,… we adapt a nature inspired strategy (FA) to find the approximate solution, by using the above Lemma 3.1. for constructing the fitness function formulated as:
E(ξi)=minm∑j=1|L−1(A(g(tj)))+L−1(n∑i=0tjiαΓ(iα+1)ξi)+L−1(N(X0(tj),tj))| | (18) |
Theorem 3.3. Let x(t) be an integrable and continuous function defined in some domain [a,T] where a∈ℜ, tiα is bounded and continuous function defined in the same interval [a,T] for some, ε∈ℜ+ and real bounded weights, ξi∈ℜ then by substituting the real bounded weights attained by the FA, the obtained series solution converges as n approaches to infinity.
By simple calculations and assuming |ξi|⩽ˆM; ˆM∈ℜ+ then Eq. (14) can be written as:
|x(t)|=|∫Tan∑i=0tiαξiΓ(iα+1)dt|⩽ˆMε(T−a)Γ((iα+1)+1) | (19) |
As the number of terms in the initial approximation, n→∞ it leads|Γ((nα+1)+1)|>1.
The FA is based on the flashing pattern and behavior of the fireflies. These flies are unisex, so they are attracted to the other fireflies regardless of their gender. The attractiveness of these tropical fireflies is proportional to their brightness. The main purpose of this attraction is to enable an algorithm to converge quickly, by allowing the swarming agents to interact and move towards the true global optimality. In FA the attractiveness between the fireflies at βi and βj are determined as
βk+1i=βki+θ(βi,βj)(βkj−βki)+˜α(μ−12) | (20) |
where the middle term is due to the attraction θ(βi,βj) of the fireflies, which varies with the distance rij between them and can be modeled as:
θ(βi,βj)=θ0e−γr2ij | (21) |
where
rij=‖βi−βj‖ | (22) |
The attractiveness at distance zero is represented by θ0, ˜α∈[0,1] is a randomization parameter and μ is a random number generator that is uniformly distributed in the interval [0,1]. The light absorption coefficient is denoted by γ∈(0,∞], which is assumed to be very large but in practice it is determined by the characteristic distance Γ of the system, over which the attractiveness varies from θ0 to θ0e−1. The convergence of the algorithm can be further improved by varying the randomization parameter ˜α that is decreased gradually as the optima is approached.
The nonlinear updating equation used by the FA yields a richer behavior and higher convergence than the other optimization algorithms with linear updating equating. The middle term of the updating equation becomes negligible by setting a very large value of γ, leading to the standard SA algorithm. For ˜α=0 and γ→0 the exponential term in Eq. (21) tends to one, leading Eq. (20) to be a variant of the DE algorithm. Moreover, the APSO algorithm and the harmonic search (HS) algorithm are also the special cases of the FA. Therefore, we can essentially say that the FA is an amalgamation of all these four algorithm (SA, DE, APSO, HS) to a certain extent and can easily outperform other modern and powerful optimization algorithms [33,34].
The step by step procedure of the proposed algorithm for Eq. (1) is as follows:
Step1: Construct homotopy X(r(t),δ):Ω×[0,T]→ℜ by using Eqs. (9) and (11).
Step2: Fix the number of terms in the initial approximation x0(t).
Step3: Set X1(t)=0 and the equidistant points for in the interval
Step4:Configure the error function E(ξi) specified in Eq. (15).
Step5: Utilize the FA with γ→∞, to optimize the derived error function and attain the minimum fitness value for an effective value of each unknown.
Step6: Substitute the optimal values of ξi in Eq. (14) to acquire the approximate solution.
Step7: Calculate error norms to validate to correctness and efficiency of the suggested scheme.
In this section, the accuracy and competency of the proposed scheme is illustrated by considering several examples of the form (1). Comparison is made with the available exact solution and the results obtained by some former techniques such as, ABMA, APSO, MHPM and RPSM. Also, the accuracy of the algorithm are assessed in terms of the error norms mathematically formulated as
Labs=|xrfj−x∗j|, | (23) |
L∞=maxLabs, | (24) |
Lrms=√1ˆnˆn∑j=1|xrfj−x∗j|2, | (25) |
where ˆn stands for the number of input grid points in the computational domain [t0,T], for T∈ℜ+ and xrfj and x∗j represents the reference solution and the approximate solution, respectively.
Test problem 1
Consider a nonlinear fractional differential equation
cDα0x(t)=−x(t)+x2(t),0<α⩽1, | (26) |
with initial condition x(0)=0.5. The exact solution of above problem for α=1 is given as
x(t)=e−te−t+1. | (27) |
Test problem 2
Consider a nonlinear factional Riccati differential equation with the fractional derivative defined in the Caputo sense expressed as
CDα0x(t)=−x2(t)+2x(t)+1,0<α⩽1 | (28) |
and the initial condition defined as x(0)=0. The exact solution of the above problem for α=1 is given as
x(t)=1+√2tanh(√2t+12log(√2−1√2+1)). | (29) |
Test problem 3
Now we consider another fractional differential equation
CDα0x(t)=−x2(t)+1,0<α⩽1, | (30) |
with initial condition x(0)=0 and the exact solution of Eq. (27) for α=1 is
x(t)=e2t−1e2t+1. | (31) |
Test problem 4
Consider another fractional differential equation
CDα0x(t)=x(t)t+x2(t)t+2t2−αΓ(3−α)−t3(1+t2),0<α⩽1, | (32) |
with initial condition x(0)=0. By the use of definition (2.1) the exact solution of Eq. (32) is given as
x(t)=t2. | (33) |
The proposed methodology EHPM accompanied by an efficient and powerful optimization technique is applied to the nonlinear problems defined above. The fitness function derived in the form (18) for the nonlinear problems (26), (28), (30) and (32) are optimized using the FA. The optimal weights achieved by the implementation of the developed technique at α=1, n=5 and x∈[0,1] yielding a fitness value 3.177187×10−17, 4.316209×10−12, and 1.43145×10−17 for the test problems 1–4 are displayed in Figure 2, respectively. The numerical results accomplished by the utilization of these optimal weights for the respective test problems are presented in Tables 1–4. Comparative study is conducted by comparing the EHPM computed results with the respective exact solution and the results attained by ABMA with the step size ˉh=0.001 and ˉK=1000, 10 term approximate solution by MHPM, the APSO computed results (for 800 iterations and swarm size equal to 40) and the seventh term approximate solutions obtained by the RPSM. To illustrate the reliability and accuracy of the suggested scheme EHPM the values of absolute error attained by the considering the exact solution as the reference solution of the respective problems 1–4 are also presented in Tables 1–4, respectively. It is observed that generally, the EHPM computed results for all the problems coincides the exact solution within five to ten decimal places of accuracy. The error norms Lrms and L∞ acquired for the above problems by consuming the developed scheme EHPM for the computational domain x∈[0,1] and at distinct values n are tabulated in Table 5. These error estimates presented in Table 5 ensures about the convergence and accuracy of the deliberated technique which is seen to be increased as n increases.
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.474993 | 2.46866×10-9 |
0.2 | 0.450166 | 0.450166 | 0.450166 | 0.450166 | 0.450165 | 0.449984 | 8.11591×10-9 |
0.3 | 0.425557 | 0.425557 | 0.425557 | 0.425557 | 0.425552 | 0.425056 | 2.62871×10-9 |
0.4 | 0.401312 | 0.401312 | 0.401312 | 0.401312 | 0.401291 | 0.400332 | 6.47150×10-9 |
0.5 | 0.377541 | 0.377541 | 0.377541 | 0.377541 | 0.377476 | 0.375992 | 1.62494×10-9 |
0.6 | 0.354344 | 0.354344 | 0.354344 | 0.354344 | 0.354182 | 0.352276 | 6.16224×10-9 |
0.7 | 0.331812 | 0.331812 | 0.331812 | 0.331813 | 0.331462 | 0.329493 | 6.81599×10-10 |
0.8 | 0.310026 | 0.310026 | 0.310026 | 0.310028 | 0.309343 | 0.308033 | 6.05932×10-9 |
0.9 | 0.289050 | 0.289051 | 0.289051 | 0.289058 | 0.287820 | 0.288369 | 3.75236×10-9 |
1.0 | 0.268941 | 0.268941 | 0.268941 | 0.268961 | 0.266858 | 0.271073 | 4.53833×10-9 |
t | yEXACT | yEHPM | yABMA | yRPSM | yNHPM | yAPSO | Pabs |
0.1 | 0.110295 | 0.110323 | 0.110295 | 0.110295 | 0.110294 | 0.132598 | 2.79559×10-5 |
0.2 | 0.241976 | 0.241942 | 0.241976 | 0.241976 | 0.241965 | 0.275045 | 3.40483×10-5 |
0.3 | 0.395104 | 0.395107 | 0.395104 | 0.395089 | 0.395106 | 0.427310 | 2.56995×10-6 |
0.4 | 0.567812 | 0.567860 | 0.567811 | 0.56766 | 0.568115 | 0.589222 | 4.80912×10-5 |
0.5 | 0.756014 | 0.756031 | 0.756014 | 0.755134 | 0.757564 | 0.760469 | 1.68103×10-5 |
0.6 | 0.953566 | 0.953530 | 0.953565 | 0.949964 | 0.958259 | 0.940590 | 3.57596×10-5 |
0.7 | 1.152949 | 1.152943 | 1.152949 | 1.141423 | 1.163459 | 1.128970 | 5.97971×10-6 |
0.8 | 1.346364 | 1.346424 | 1.346363 | 1.315723 | 1.365240 | 1.324840 | 6.04658×10-5 |
0.9 | 1.526911 | 1.526896 | 1.526911 | 1.456545 | 1.554960 | 1.527270 | 1.51114×10-5 |
1.0 | 1.689498 | 1.689546 | 1.689498 | 1.546030 | 1.723810 | 1.735140 | 4.79364×10-5 |
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.099668 | 0.099666 | 0.099668 | 0.099668 | 0.099668 | 0.103891 | 1.23113×10-6 |
0.2 | 0.197375 | 0.197378 | 0.197375 | 0.197375 | 0.197375 | 0.202649 | 2.55077×10-6 |
0.3 | 0.291313 | 0.291313 | 0.291313 | 0.291312 | 0.291312 | 0.295810 | 1.45452×10-7 |
0.4 | 0.379949 | 0.379946 | 0.379949 | 0.379944 | 0.379944 | 0.382931 | 2.41309×10-6 |
0.5 | 0.462117 | 0.462117 | 0.462117 | 0.462078 | 0.462078 | 0.463587 | 1.66903×10-7 |
0.6 | 0.537050 | 0.537052 | 0.537049 | 0.536857 | 0.536857 | 0.537380 | 2.46396×10-6 |
0.7 | 0.604368 | 0.604368 | 0.604368 | 0.603631 | 0.603631 | 0.603935 | 4.64338×10-7 |
0.8 | 0.664037 | 0.664034 | 0.664037 | 0.661706 | 0.661706 | 0.662905 | 2.81496×10-6 |
0.9 | 0.716298 | 0.716299 | 0.716298 | 0.709919 | 0.709919 | 0.713971 | 1.30183×10-6 |
1.0 | 0.761594 | 0.761592 | 0.761594 | 0.746032 | 0.746032 | 0.756838 | 2.22537×10-6 |
t | yEXACT | yEHPM | yABMA | yMHPM | yAPSO | Pabs |
0.1 | 0.01 | 0.01 | 0.019966 | 0.001 | 0.012254 | 2.91364×10-13 |
0.2 | 0.04 | 0.04 | 0.079449 | 0.008 | 0.044342 | 1.58054×10-13 |
0.3 | 0.09 | 0.09 | 0.177089 | 0.027 | 0.096103 | 2.13080×10-13 |
0.4 | 0.16 | 0.16 | 0.310209 | 0.064 | 0.167389 | 4.90163×10-14 |
0.5 | 0.25 | 0.25 | 0.473957 | 0.125 | 0.258064 | 2.18325×10-13 |
0.6 | 0.36 | 0.36 | 0.659526 | 0.216 | 0.368016 | 8.49321×10-14 |
0.7 | 0.49 | 0.49 | 0.850389 | 0.343 | 0.497157 | 1.76026×10-13 |
0.8 | 0.64 | 0.64 | 1.014271 | 0.512 | 0.645426 | 1.66533×10-13 |
0.9 | 0.81 | 0.81 | 1.086366 | 0.729 | 0.812795 | 2.24820×10-13 |
1.0 | 1.00 | 1.00 | 0.938449 | 1.000 | 0.999274 | 3.39728×10-14 |
n | Test Problem 1 | Test Problem 2 | Test Problem 3 | Test Problem 4 | ||||
Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | |
5 | 4.82×10-9 | 8.11×10-9 | 3.46×10-5 | 6.04×10-5 | 1.86×10-6 | 2.81×10-6 | 1.79×10-13 | 2.91×10-13 |
4 | 1.98×10-7 | 2.63×10-7 | 1.67×10-4 | 2.63×10-4 | 1.78×10-5 | 2.36×10-5 | 1.73×10-14 | 3.25×10-14 |
3 | 1.86×10-6 | 3.04×10-6 | 1.12×10-3 | 2.01×10-3 | 3.33×10-5 | 6.10×10-5 | 6.49×10-9 | 1.13×10-8 |
2 | 2.94×10-5 | 5.27×10-5 | 3.51×10-3 | 6.39×10-3 | 1.03×10-3 | 1.76×10-3 | 5.17×10-9 | 8.62×10-9 |
Numerical solutions of the above problems are also constructed by the discussed technique for the non-integer order. The solutions of the above nonlinear problems are attained by the learning of unknown weights in the derived fitness function at two distinct values of α i.e., 0.95 and 0.85. The optimal weights achieved by the proposed scheme for the considered fractional values with n=2 and x∈[0,1] for problem 1, with n=2 and x∈[0,4] for problem 2, with n=5 and x∈[0,5] for problem 3 and with n=5 and x∈[0,1] for problem 4 are showcased in Figure 3, respectively. The results obtained by utilizing these real valued and bounded optimal weights for the respective problems 1–4 are presented in Tables 6–9, respectively. The results obtained by EHPM are compared with the available exact solution, MHPM computed results, seventh term approximate solutions obtained by the RPSM and the results attained by ABMA for the step ˉh=0.001 and the values ˉK equal to 4000 and 5000 for the test problems 1–3, respectively. In problem 1 by taking smaller value n for a smaller domain [0,1] the results achieved by the EHPM and the RPSM shows a constructive agreement with the results obtained by ABMA. While, by taking a smaller value n for a larger domain [0,4] in problem 2, the RPSM is seen to diverge whereas, the solutions obtained by the EHPM remains convergent. By taking a larger value n for a larger domain [0,5] in problem 3 the RPSM computed results diverges whereas, the results accomplished by the proposed scheme shows a constructive agreement with the results attained by ABMA. The values of absolute error achieved by EHPM for the problems 1–3 with the results attained by ABMA as the reference solution are also depicted in Tables 6–8. In problem 4 by taking a larger value of n for a smaller domain [0,1] the EHPM computed results are compared with available exact solution and the results obtained by MHPM, with the values of absolute error tabulated in Table 9. Furthermore, the correctness and reliability of the suggested scheme is validated by the minimize fitness function and the related graphical solution attained at distinct values of α = 1, 0.95, 0.85 for a lager span, that are portrayed in Figures 4a–7a and Figures 4b–7b for the problems 1–4, respectively. One can infer that the accuracy of the presented scheme can be further enhanced but at the cost of more computation.
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.1 | 0.458464 | 0.462742 | 0.462742 | 4.27763×10-3 | 0.470459 | 0.471407 | 0.471407 | 9.4861×10-4 |
0.2 | 0.425630 | 0.433189 | 0.433189 | 7.55906×10-3 | 0.443171 | 0.444939 | 0.444939 | 1.76834×10-3 |
0.3 | 0.395869 | 0.406386 | 0.406387 | 1.05170×10-2 | 0.416934 | 0.419474 | 0.419474 | 2.54007×10-3 |
0.4 | 0.368265 | 0.381546 | 0.381546 | 1.32809×10-2 | 0.391583 | 0.394880 | 0.394880 | 3.29609×10-3 |
0.5 | 0.342440 | 0.358340 | 0.358340 | 1.58994×10-2 | 0.367092 | 0.371144 | 0.371145 | 4.05246×10-3 |
0.6 | 0.318209 | 0.336593 | 0.336596 | 1.83840×10-2 | 0.343479 | 0.348290 | 0.348291 | 4.81164×10-3 |
0.7 | 0.295471 | 0.316194 | 0.316201 | 2.07229×10-2 | 0.320785 | 0.326348 | 0.326350 | 5.56302×10-3 |
0.8 | 0.274172 | 0.297059 | 0.297079 | 2.28877×10-2 | 0.299065 | 0.305348 | 0.305353 | 6.28306×10-3 |
0.9 | 0.254283 | 0.279120 | 0.279167 | 2.48372×10-2 | 0.278377 | 0.285313 | 0.285327 | 6.93573×10-3 |
1.0 | 0.235793 | 0.262314 | 0.262417 | 2.65208×10-2 | 0.258786 | 0.266259 | 0.266294 | 7.47296×10-3 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.4 | 1.128780 | 0.772076 | 0.701892 | 3.56704×10-1 | 0.951216 | 0.629409 | 0.609030 | 3.21808×10-1 |
0.8 | 1.732609 | 1.534233 | 1.292070 | 1.98376×10-1 | 1.582695 | 1.414490 | 1.339680 | 1.68205×10-1 |
1.2 | 2.095352 | 1.965731 | -0.689702 | 1.29622×10-1 | 1.999238 | 1.962805 | 1.116260 | 3.64331×10-2 |
1.6 | 2.300624 | 2.159265 | -12.26610 | 1.41359×10-1 | 2.252742 | 2.217871 | -1.984430 | 3.48714×10-2 |
2.0 | 2.399864 | 2.248003 | -48.77540 | 1.51861×10-1 | 2.384753 | 2.319886 | -9.224490 | 6.48668×10-2 |
2.4 | 2.430901 | 2.293382 | -139.1670 | 1.37519×10-1 | 2.431978 | 2.361076 | -17.66350 | 7.09018×10-2 |
2.8 | 2.424109 | 2.319480 | -332.8130 | 1.04629×10-1 | 2.428249 | 2.379249 | -13.99210 | 4.89996×10-2 |
3.2 | 2.405166 | 2.336069 | -707.8470 | 6.90973×10-2 | 2.405453 | 2.388335 | 34.79380 | 1.71183×10-2 |
3.6 | 2.396527 | 2.347472 | -1381.040 | 4.90542×10-2 | 2.394053 | 2.393514 | 194.3750 | 5.39429×10-4 |
4.0 | 2.418283 | 2.355791 | -2519.200 | 6.24925×10-2 | 2.423411 | 2.396826 | 580.2210 | 2.65850×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.5 | 0.571511 | 0.515900 | 0.481043 | 5.56103×10-2 | 0.497222 | 0.481137 | 0.515392 | 1.60848×10-2 |
1.0 | 0.824167 | 0.749435 | 0.731907 | 7.47319×10-2 | 0.778543 | 0.758713 | 0.675115 | 1.98303×10-2 |
1.5 | 0.933117 | 0.855319 | 0.263158 | 7.77975×10-2 | 0.915829 | 0.888118 | -0.334334 | 2.77111×10-2 |
2.0 | 0.972671 | 0.906026 | -4.400130 | 6.66455×10-2 | 0.968915 | 0.943982 | -7.027070 | 2.49329×10-2 |
2.5 | 0.986094 | 0.932703 | -25.95750 | 5.33910×10-2 | 0.985180 | 0.968440 | -32.40280 | 1.67394×10-2 |
3.0 | 0.993557 | 0.948165 | -97.30430 | 4.53923×10-2 | 0.992469 | 0.979829 | -104.6460 | 1.26404×10-2 |
3.5 | 0.999402 | 0.957932 | -288.1960 | 4.14707×10-2 | 0.999263 | 0.985618 | -275.6610 | 1.36449×10-2 |
4.0 | 1.000703 | 0.964556 | -728.4980 | 3.61731×10-2 | 1.001750 | 0.988867 | -631.4190 | 1.28808×10-2 |
4.5 | 0.997296 | 0.969315 | -1638.840 | 2.79813×10-2 | 0.997197 | 0.990876 | -1304.100 | 6.32101×10-3 |
5.0 | 1.002409 | 0.972892 | -3369.480 | 2.95975×10-2 | 1.003140 | 0.992229 | -2486.000 | 1.09120×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yEXACT | yMHPM | Pabs | yEHPM | yEXACT | yMHPM | Pabs | |
0.1 | 0.009336 | 0.01 | 0.001527 | 6.63272×10-4 | 0.009782 | 0.01 | 0.001150 | 2.17468×10-4 |
0.2 | 0.037202 | 0.04 | 0.011010 | 2.79761×10-3 | 0.039123 | 0.04 | 0.008892 | 8.76525×10-4 |
0.3 | 0.083691 | 0.09 | 0.034966 | 6.30802×10-3 | 0.087988 | 0.09 | 0.029411 | 2.01164×10-3 |
0.4 | 0.148502 | 0.16 | 0.079383 | 1.14971×10-2 | 0.156347 | 0.16 | 0.068718 | 3.65298×10-3 |
0.5 | 0.231266 | 0.25 | 0.149941 | 1.87333×10-2 | 0.244079 | 0.25 | 0.132726 | 5.92024×10-3 |
0.6 | 0.331497 | 0.36 | 0.251096 | 2.85027×10-2 | 0.351006 | 0.36 | 0.227269 | 8.99314×10-3 |
0.7 | 0.448290 | 0.49 | 0.391189 | 4.17099×10-2 | 0.476835 | 0.49 | 0.358125 | 1.31646×10-2 |
0.8 | 0.579872 | 0.64 | 0.572354 | 6.01272×10-2 | 0.621014 | 0.64 | 0.531020 | 1.89856×10-2 |
0.9 | 0.723061 | 0.81 | 0.800662 | 8.69390×10-2 | 0.782501 | 0.81 | 0.751641 | 2.74983×10-2 |
1.0 | 0.872643 | 1.00 | 1.081081 | 1.27356×10-1 | 0.959450 | 1.00 | 1.025641 | 4.05497×10-2 |
In this work, the nonlinear fractional differential equations with the fractional derivative and integral operators defined in the Caputo sense are successfully solved by a simple, accurate and reliable numerical technique. The designed methodology EHPM is an amalgamation of the classical homotopy perturbation technique with the modern bio-inspired metaheuristic technique. The FA which is based on the flashing behavior of fireflies fast track the procedure to determine the approximate solution of the considered nonlinear FDE. Some instructive examples were undertaken, with varying span and different number of terms in the assumed initial approximation to expound the potential ability of the proposed technique. The validity, applicability and computational efficiency of the EHPM is exposed by comparing the fast track outcomes achieved with the available exact solution and the solutions obtained by ABMA, APSO, MHPM and the RPSM. High accuracy and consistent convergence were ascertained by the accomplishment of the optimal values of the error norms. Thus, the contribution of the discussed scheme is abridged as follows:
● The expedite design methodology is well suited to solve the nonlinear FDE.
● The classical approach transformed the nonlinear FDE into an algebraic system which leads to a weighted series solution.
● Fitness function determined by a trivial fragment of the weighted series solution was optimized using the FA.
● Fast track outcomes were achieved by the modern optimization technique (FA).
● Numerical experiments ratified the correctness and reliability of the EHPM.
● Converging solutions were attained even for a larger computational domain.
● High accuracy and consistent convergence achieved by EHPM can be further enhanced by increasing the number of term in the initial approximation.
As the proposed design methodology exhibits a merger of classical idea with the modern optimizing metaheuristic technique, which offers promising converging results. One may utilize the recommended fusion to solve numerous other mathematical models or may even attempt to determine a new combination of modern optimization technique with the former traditional schemes.
The authors declare no conflict of interest.
[1] | Vunjak-Novakovic G, Tandon N, Godier A, et al. (2009) Challenges in cardiac tissue engineering. Tissue Eng Part B 16: 169–187. |
[2] |
Leone M, Magadum A, Engel FB (2015) Cardiomyocyte proliferation in cardiac development and regeneration: a guide to methodologies and interpretations. Am J Physiol-Heart C 309: H1237–H1250. doi: 10.1152/ajpheart.00559.2015
![]() |
[3] |
Mummery CL, Zhang J, Ng ES, et al. (2012) Differentiation of human embryonic stem cells and induced pluripotent stem cells to cardiomyocytes: a method overview. Circ Res 111: 344–358. doi: 10.1161/CIRCRESAHA.110.227512
![]() |
[4] |
Schwach V, Passier R (2016) Generation and purification of human stem cell-derived cardiomyocytes. Differentiation 91: 126–138. doi: 10.1016/j.diff.2016.01.001
![]() |
[5] |
Kempf H, Andree B, Zweigerdt R (2016) Large-scale production of human pluripotent stem cell derived cardiomyocytes. Adv Drug Deliver Rev 96: 18–30. doi: 10.1016/j.addr.2015.11.016
![]() |
[6] |
Duan Y, Liu Z, O'Neill J, et al. (2011) Hybrid gel composed of native heart matrix and collagen induces cardiac differentiation of human embryonic stem cells without supplemental growth factors. J Cardiovasc Transl 4: 605–615. doi: 10.1007/s12265-011-9304-0
![]() |
[7] |
Higuchi S, Lin Q, Wang J, et al. (2013) Heart extracellular matrix supports cardiomyocyte differentiation of mouse embryonic stem cells. J Biosci Bioeng 115: 320–325. doi: 10.1016/j.jbiosc.2012.10.004
![]() |
[8] |
Engler AJ, Sen S, Sweeney HL, et al. (2006) Matrix elasticity directs stem cell lineage specification. Cell 126: 677–689. doi: 10.1016/j.cell.2006.06.044
![]() |
[9] |
Chaudhuri O, Gu L, Darnell M, et al. (2015) Substrate stress relaxation regulates cell spreading. Nat Commun 6: 6364. doi: 10.1038/ncomms7364
![]() |
[10] |
Brandl F, Sommer F, Goepferich A (2007) Rational design of hydrogels for tissue engineering: impact of physical factors on cell behavior. Biomaterials 28: 134–146. doi: 10.1016/j.biomaterials.2006.09.017
![]() |
[11] |
Feng Z, Wagatsuma Y, Kikuchi M, et al. (2014) The mechanisms of fibroblast-mediated compaction of collagen gels and the mechanical niche around individual fibroblasts. Biomaterials 35: 8078–8091. doi: 10.1016/j.biomaterials.2014.05.072
![]() |
[12] |
Feng Z, Ishiguro Y, Fujita K, et al. (2015) A fibril-based structural constitutive theory reveals the dominant role of network characteristics on the mechanical behavior of fibroblast-compacted collagen gels. Biomaterials 67: 365–381. doi: 10.1016/j.biomaterials.2015.07.038
![]() |
[13] |
Johnson TD, Lin SY, Christman KL (2011) Tailoring material properties of a nanofibrous extracellular matrix derived hydrogel. Nanotechnology 22: 494015. doi: 10.1088/0957-4484/22/49/494015
![]() |
[14] |
Freytes DO, Martin J, Velankar SS, et al. (2008) Preparation and rheological characterization of a gel form of the porcine urinary bladder matrix. Biomaterials 29: 1630–1637. doi: 10.1016/j.biomaterials.2007.12.014
![]() |
[15] |
Williams C, Budina E, Stoppel WL, et al. (2015) Cardiac Extracellular Matrix-Fibrin Hybrid Scaffolds with Tunable Properties for Cardiovascular Tissue Engineering. Acta Biomater 14: 84–95. doi: 10.1016/j.actbio.2014.11.035
![]() |
[16] |
Grover GN, Rao N, Christman KL (2014) Myocardial Matrix-Polyethylene Glycol Hybrid Hydrogels for Tissue Engineering. Nanotechnology 25: 014011. doi: 10.1088/0957-4484/25/1/014011
![]() |
[17] |
Singelyn JM, Christman KL (2011) Modulation of material properties of a decellularized myocardial matrix scaffold. Macromol Biosci 11: 731–738. doi: 10.1002/mabi.201000423
![]() |
[18] | Everaerts F, Torrianni M, Hendriks M, et al. (2008) Biomechanical properties of carbodiimide crosslinked collagen: Influence of the formation of ester crosslinks. J Biomed Mater Res A 85: 547–555. |
[19] | Fujita K, Tuchida Y, Seki H, et al. (2015) Characterizing and modulating the mechanical properties of hydrogels from ventricular extracellular matrix. 2015 10th Asian Control Conference (ASCC) 1–5. |
[20] | Kikuchi M, Feng Z, Kosawada T, et al. (2016) Stress relaxation and stress-strain characteristics of porcine amniotic membrane. Bio-med Mater Eng 27: 603–611. |
[21] |
Johnson TD, DeQuach JA, Gaetani R, et al. (2014) Human versus porcine tissue sourcing for an injectable myocardial matrix hydrogel. Biomater Sci 2: 735–744. doi: 10.1039/C3BM60283D
![]() |
[22] |
Ishihama Y, Oda Y, Tabata T, et al. (2005) Exponentially Modified Protein Abundance Index (emPAI) for Estimation of Absolute Protein Amount in Proteomics by the Number of Sequenced Peptides per Protein. Mol Cell Proteomics 4: 1265–1272. doi: 10.1074/mcp.M500061-MCP200
![]() |
[23] | Fung YC (1993) Biomechanics: Mechanical properties of living tissues, 2 Eds, Springer. |
[24] |
Fancey KS (2005) A mechanical model for creep, recovery and stress relaxation in polymeric materials. J Mater Sci 40: 4827–4831. doi: 10.1007/s10853-005-2020-x
![]() |
[25] |
Feng Z, Matsumoto T, Nakamura T, et al. (2003) Measurements of the mechanical properties of contracted collagen gels populated with rat fibroblasts or cardiomyocytes. J Artif Organs 6: 192–196. doi: 10.1007/s10047-003-0230-z
![]() |
[26] |
Ott HC, Matthiesen TS, Goh SK, et al. (2008) Perfusion-decellularized matrix: using nature's platform to engineer a bioartificial heart. Nat Med 14: 213–221. doi: 10.1038/nm1684
![]() |
[27] |
Eiselt P, Lee KY, Mooney DJ (1999) Rigidity of Two-Component Hydrogels Prepared from Alginate and Poly(ethylene glycol)-Diamines. Macromolecules 32: 5561–5566. doi: 10.1021/ma990514m
![]() |
[28] |
Lee KY, Rowley JA, Eiselt P, et al. (2000) Controlling Mechanical and Swelling Properties of Alginate Hydrogels Independently by Cross-Linker Type and Cross-Linking Density. Macromolecules 33: 4291–4294. doi: 10.1021/ma9921347
![]() |
[29] |
Wang Y, Wang G, Luo X, et al. (2012) Substrate stiffness regulated the proliferation, migration, and differentiation of epidermal cells. Burns 38: 414–420. doi: 10.1016/j.burns.2011.09.002
![]() |
[30] |
Ungerleider JL, Johnson TD, Rao N, et al. (2015) Fabrication and Characterization of Injectable Hydrogels Derived from Decellularized Skeletal and Cardiac Muscle. Methods 84: 53–59. doi: 10.1016/j.ymeth.2015.03.024
![]() |
[31] |
Johnson TD, Dequach JA, Gaetani R, et al. (2014) Human versus porcine tissue sourcing for an injectable myocardial matrix hydrogel. Biomater Sci 2: 735–744. doi: 10.1039/C3BM60283D
![]() |
[32] |
Seif-Naraghi SB, Horn D, Schup-Magoffin PA, et al. (2011) Patient-to-patient variability in autologous pericardial matrix scaffolds for cardiac repair. J Cardiovasc Transl 4: 545–556. doi: 10.1007/s12265-011-9293-z
![]() |
[33] |
Seif-Naraghi SB, Horn D, Schup-Magoffin PJ, et al. (2012) Injectable extracellular matrix derived hydrogel provides a platform for enhanced retention and delivery of a heparin-binding growth factor. Acta Biomater 8: 3695–3703. doi: 10.1016/j.actbio.2012.06.030
![]() |
[34] |
Huang-Lee LL, Cheung DT, Nimni ME (1990) Biochemical changes and cytotoxicity associated with the degradation of polymeric glutaraldehyde derived crosslinks. J Biomed Mater Res A 24: 1185–1201. doi: 10.1002/jbm.820240905
![]() |
[35] | Cornwell KG, Lei P, Andreadis ST, et al. (2007) Crosslinking of discrete self-assembled collagen threads: Effects on mechanical strength and cell-matrix interactions. J Biomed Mater Res A 80: 362–371. |
[36] |
Yunoki S, Mori K, Suzuki T, et al. (2007) Novel elastic material from collagen for tissue engineering. J Mater Sci-Mater M 18: 1369–1375. doi: 10.1007/s10856-007-0121-6
![]() |
[37] |
Park SN, Park JC, Kim HO, et al. (2002) Characterization of porous collagen/hyaluronic acid scaffold modified by 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide cross-linking. Biomaterials 23: 1205–1212. doi: 10.1016/S0142-9612(01)00235-6
![]() |
[38] |
Ma DH, Lai JY, Cheng HY, et al. (2010) Carbodiimide cross-linked amniotic membranes for cultivation of limbal epithelial cells. Biomaterials 31: 6647–6658. doi: 10.1016/j.biomaterials.2010.05.034
![]() |
[39] |
Piao Y, Chen B (2017) Synthesis and mechanical properties of double cross-linked gelatin-graphene oxide hydrogels. Int J Biol Macromol 101: 791–798. doi: 10.1016/j.ijbiomac.2017.03.155
![]() |
[40] |
Tanpichaia S, Oksmana K (2016) Cross-linked nanocomposite hydrogels based on cellulose nanocrystals and PVA: Mechanical properties and creep recovery. Compos Part A-Appl S 88: 226–233. doi: 10.1016/j.compositesa.2016.06.002
![]() |
[41] | DeQuach J, Lin J, Cam C, et al. (2013) Injectable skeletal muscle matrix hydrogel promotes neovascularization and muscle cell infiltration in a hindlimb ischemia model. Eur Cells Mater 23: 400–412. |
[42] |
Massensini AR, Ghuman H, Saldin LT, et al. (2015) Concentration-dependent rheological properties of ECM hydrogel for intracerebral delivery to a stroke cavity. Acta Biomater 27: 116–130. doi: 10.1016/j.actbio.2015.08.040
![]() |
[43] |
Wolf MT, Daly KA, Brennan-Pierce EP, et al. (2012) A hydrogel derived from decellularized dermal extracellular matrix. Biomaterials 33: 7028–7038. doi: 10.1016/j.biomaterials.2012.06.051
![]() |
[44] |
Lee JS, Shin J, Park HM, et al. (2014) Liver extracellular matrix providing dual functions of two-dimensional substrate coating and three-dimensional injectable hydrogel platform for liver tissue engineering. Biomacromolecules 15: 206–218. doi: 10.1021/bm4015039
![]() |
[45] |
Pouliot RA, Link PA, Mikhaiel NS, et al. (2016) Development and characterization of a naturally derived lung extracellular matrix hydrogel. J Biomed Mater Res A 104: 1922–1935. doi: 10.1002/jbm.a.35726
![]() |
[46] | Ma L, Gao C, Mao Z, et al. (2004) Biodegradability and cell-mediated contraction of porous collagen scaffolds: the effect of lysine as a novel crosslinking bridge. J Biomed Mater Res A 71: 334–342. |
[47] |
Wang HB, Dembo M, Wang YL (2000) Substrate flexibility regulate growth and apoptosis of normal but not transformed cells. Am J Physiol-Cell Ph 279: C1345–C1350. doi: 10.1152/ajpcell.2000.279.5.C1345
![]() |
[48] |
Engler AJ, Carag-Kriegr C, Johnson CP, et al. (2008) Embryonic cardiomyocytes beat best on a matrix with heart-like elasticity: scar-like rigidity inhibits beating. J Cell Sci 121: 3794–3802. doi: 10.1242/jcs.029678
![]() |
1. | Najeeb Alam Khan, Samreen Ahmed, Oyoon Abdul Razzaq, Ahmad Kamil Mahmood, Exploring fractional order 2‐D Helmholtz equation using finite difference scheme through the bat optimization algorithm, 2021, 0170-4214, 10.1002/mma.7271 | |
2. | Najeeb Alam Khan, Samreen Ahmed, Finite Difference Method with Metaheuristic Orientation for Exploration of Time Fractional Partial Differential Equations, 2021, 7, 2349-5103, 10.1007/s40819-021-01061-y | |
3. | Moa’ath N. Oqielat, Tareq Eriqat, Zeyad Al-Zhour, Osama Ogilat, Ahmad El-Ajou, Ishak Hashim, Construction of fractional series solutions to nonlinear fractional reaction–diffusion for bacteria growth model via Laplace residual power series method, 2022, 2195-268X, 10.1007/s40435-022-01001-8 | |
4. | Mohamed A. Abd El Salam, Mohamed A. Ramadan, Mahmoud A. Nassar, Praveen Agarwal, Yu-Ming Chu, Matrix computational collocation approach based on rational Chebyshev functions for nonlinear differential equations, 2021, 2021, 1687-1847, 10.1186/s13662-021-03481-y | |
5. | Ahmed B. Khoshaim, Muhammad Naeem, Ali Akgul, Nejib Ghanmi, Shamsullah Zaland, Lakhdar Ragoub, Novel Analysis of Fractional-Order Fifth-Order Korteweg–de Vries Equations, 2022, 2022, 2314-4785, 1, 10.1155/2022/1883268 |
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.474993 | 2.46866×10-9 |
0.2 | 0.450166 | 0.450166 | 0.450166 | 0.450166 | 0.450165 | 0.449984 | 8.11591×10-9 |
0.3 | 0.425557 | 0.425557 | 0.425557 | 0.425557 | 0.425552 | 0.425056 | 2.62871×10-9 |
0.4 | 0.401312 | 0.401312 | 0.401312 | 0.401312 | 0.401291 | 0.400332 | 6.47150×10-9 |
0.5 | 0.377541 | 0.377541 | 0.377541 | 0.377541 | 0.377476 | 0.375992 | 1.62494×10-9 |
0.6 | 0.354344 | 0.354344 | 0.354344 | 0.354344 | 0.354182 | 0.352276 | 6.16224×10-9 |
0.7 | 0.331812 | 0.331812 | 0.331812 | 0.331813 | 0.331462 | 0.329493 | 6.81599×10-10 |
0.8 | 0.310026 | 0.310026 | 0.310026 | 0.310028 | 0.309343 | 0.308033 | 6.05932×10-9 |
0.9 | 0.289050 | 0.289051 | 0.289051 | 0.289058 | 0.287820 | 0.288369 | 3.75236×10-9 |
1.0 | 0.268941 | 0.268941 | 0.268941 | 0.268961 | 0.266858 | 0.271073 | 4.53833×10-9 |
t | yEXACT | yEHPM | yABMA | yRPSM | yNHPM | yAPSO | Pabs |
0.1 | 0.110295 | 0.110323 | 0.110295 | 0.110295 | 0.110294 | 0.132598 | 2.79559×10-5 |
0.2 | 0.241976 | 0.241942 | 0.241976 | 0.241976 | 0.241965 | 0.275045 | 3.40483×10-5 |
0.3 | 0.395104 | 0.395107 | 0.395104 | 0.395089 | 0.395106 | 0.427310 | 2.56995×10-6 |
0.4 | 0.567812 | 0.567860 | 0.567811 | 0.56766 | 0.568115 | 0.589222 | 4.80912×10-5 |
0.5 | 0.756014 | 0.756031 | 0.756014 | 0.755134 | 0.757564 | 0.760469 | 1.68103×10-5 |
0.6 | 0.953566 | 0.953530 | 0.953565 | 0.949964 | 0.958259 | 0.940590 | 3.57596×10-5 |
0.7 | 1.152949 | 1.152943 | 1.152949 | 1.141423 | 1.163459 | 1.128970 | 5.97971×10-6 |
0.8 | 1.346364 | 1.346424 | 1.346363 | 1.315723 | 1.365240 | 1.324840 | 6.04658×10-5 |
0.9 | 1.526911 | 1.526896 | 1.526911 | 1.456545 | 1.554960 | 1.527270 | 1.51114×10-5 |
1.0 | 1.689498 | 1.689546 | 1.689498 | 1.546030 | 1.723810 | 1.735140 | 4.79364×10-5 |
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.099668 | 0.099666 | 0.099668 | 0.099668 | 0.099668 | 0.103891 | 1.23113×10-6 |
0.2 | 0.197375 | 0.197378 | 0.197375 | 0.197375 | 0.197375 | 0.202649 | 2.55077×10-6 |
0.3 | 0.291313 | 0.291313 | 0.291313 | 0.291312 | 0.291312 | 0.295810 | 1.45452×10-7 |
0.4 | 0.379949 | 0.379946 | 0.379949 | 0.379944 | 0.379944 | 0.382931 | 2.41309×10-6 |
0.5 | 0.462117 | 0.462117 | 0.462117 | 0.462078 | 0.462078 | 0.463587 | 1.66903×10-7 |
0.6 | 0.537050 | 0.537052 | 0.537049 | 0.536857 | 0.536857 | 0.537380 | 2.46396×10-6 |
0.7 | 0.604368 | 0.604368 | 0.604368 | 0.603631 | 0.603631 | 0.603935 | 4.64338×10-7 |
0.8 | 0.664037 | 0.664034 | 0.664037 | 0.661706 | 0.661706 | 0.662905 | 2.81496×10-6 |
0.9 | 0.716298 | 0.716299 | 0.716298 | 0.709919 | 0.709919 | 0.713971 | 1.30183×10-6 |
1.0 | 0.761594 | 0.761592 | 0.761594 | 0.746032 | 0.746032 | 0.756838 | 2.22537×10-6 |
t | yEXACT | yEHPM | yABMA | yMHPM | yAPSO | Pabs |
0.1 | 0.01 | 0.01 | 0.019966 | 0.001 | 0.012254 | 2.91364×10-13 |
0.2 | 0.04 | 0.04 | 0.079449 | 0.008 | 0.044342 | 1.58054×10-13 |
0.3 | 0.09 | 0.09 | 0.177089 | 0.027 | 0.096103 | 2.13080×10-13 |
0.4 | 0.16 | 0.16 | 0.310209 | 0.064 | 0.167389 | 4.90163×10-14 |
0.5 | 0.25 | 0.25 | 0.473957 | 0.125 | 0.258064 | 2.18325×10-13 |
0.6 | 0.36 | 0.36 | 0.659526 | 0.216 | 0.368016 | 8.49321×10-14 |
0.7 | 0.49 | 0.49 | 0.850389 | 0.343 | 0.497157 | 1.76026×10-13 |
0.8 | 0.64 | 0.64 | 1.014271 | 0.512 | 0.645426 | 1.66533×10-13 |
0.9 | 0.81 | 0.81 | 1.086366 | 0.729 | 0.812795 | 2.24820×10-13 |
1.0 | 1.00 | 1.00 | 0.938449 | 1.000 | 0.999274 | 3.39728×10-14 |
n | Test Problem 1 | Test Problem 2 | Test Problem 3 | Test Problem 4 | ||||
Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | |
5 | 4.82×10-9 | 8.11×10-9 | 3.46×10-5 | 6.04×10-5 | 1.86×10-6 | 2.81×10-6 | 1.79×10-13 | 2.91×10-13 |
4 | 1.98×10-7 | 2.63×10-7 | 1.67×10-4 | 2.63×10-4 | 1.78×10-5 | 2.36×10-5 | 1.73×10-14 | 3.25×10-14 |
3 | 1.86×10-6 | 3.04×10-6 | 1.12×10-3 | 2.01×10-3 | 3.33×10-5 | 6.10×10-5 | 6.49×10-9 | 1.13×10-8 |
2 | 2.94×10-5 | 5.27×10-5 | 3.51×10-3 | 6.39×10-3 | 1.03×10-3 | 1.76×10-3 | 5.17×10-9 | 8.62×10-9 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.1 | 0.458464 | 0.462742 | 0.462742 | 4.27763×10-3 | 0.470459 | 0.471407 | 0.471407 | 9.4861×10-4 |
0.2 | 0.425630 | 0.433189 | 0.433189 | 7.55906×10-3 | 0.443171 | 0.444939 | 0.444939 | 1.76834×10-3 |
0.3 | 0.395869 | 0.406386 | 0.406387 | 1.05170×10-2 | 0.416934 | 0.419474 | 0.419474 | 2.54007×10-3 |
0.4 | 0.368265 | 0.381546 | 0.381546 | 1.32809×10-2 | 0.391583 | 0.394880 | 0.394880 | 3.29609×10-3 |
0.5 | 0.342440 | 0.358340 | 0.358340 | 1.58994×10-2 | 0.367092 | 0.371144 | 0.371145 | 4.05246×10-3 |
0.6 | 0.318209 | 0.336593 | 0.336596 | 1.83840×10-2 | 0.343479 | 0.348290 | 0.348291 | 4.81164×10-3 |
0.7 | 0.295471 | 0.316194 | 0.316201 | 2.07229×10-2 | 0.320785 | 0.326348 | 0.326350 | 5.56302×10-3 |
0.8 | 0.274172 | 0.297059 | 0.297079 | 2.28877×10-2 | 0.299065 | 0.305348 | 0.305353 | 6.28306×10-3 |
0.9 | 0.254283 | 0.279120 | 0.279167 | 2.48372×10-2 | 0.278377 | 0.285313 | 0.285327 | 6.93573×10-3 |
1.0 | 0.235793 | 0.262314 | 0.262417 | 2.65208×10-2 | 0.258786 | 0.266259 | 0.266294 | 7.47296×10-3 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.4 | 1.128780 | 0.772076 | 0.701892 | 3.56704×10-1 | 0.951216 | 0.629409 | 0.609030 | 3.21808×10-1 |
0.8 | 1.732609 | 1.534233 | 1.292070 | 1.98376×10-1 | 1.582695 | 1.414490 | 1.339680 | 1.68205×10-1 |
1.2 | 2.095352 | 1.965731 | -0.689702 | 1.29622×10-1 | 1.999238 | 1.962805 | 1.116260 | 3.64331×10-2 |
1.6 | 2.300624 | 2.159265 | -12.26610 | 1.41359×10-1 | 2.252742 | 2.217871 | -1.984430 | 3.48714×10-2 |
2.0 | 2.399864 | 2.248003 | -48.77540 | 1.51861×10-1 | 2.384753 | 2.319886 | -9.224490 | 6.48668×10-2 |
2.4 | 2.430901 | 2.293382 | -139.1670 | 1.37519×10-1 | 2.431978 | 2.361076 | -17.66350 | 7.09018×10-2 |
2.8 | 2.424109 | 2.319480 | -332.8130 | 1.04629×10-1 | 2.428249 | 2.379249 | -13.99210 | 4.89996×10-2 |
3.2 | 2.405166 | 2.336069 | -707.8470 | 6.90973×10-2 | 2.405453 | 2.388335 | 34.79380 | 1.71183×10-2 |
3.6 | 2.396527 | 2.347472 | -1381.040 | 4.90542×10-2 | 2.394053 | 2.393514 | 194.3750 | 5.39429×10-4 |
4.0 | 2.418283 | 2.355791 | -2519.200 | 6.24925×10-2 | 2.423411 | 2.396826 | 580.2210 | 2.65850×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.5 | 0.571511 | 0.515900 | 0.481043 | 5.56103×10-2 | 0.497222 | 0.481137 | 0.515392 | 1.60848×10-2 |
1.0 | 0.824167 | 0.749435 | 0.731907 | 7.47319×10-2 | 0.778543 | 0.758713 | 0.675115 | 1.98303×10-2 |
1.5 | 0.933117 | 0.855319 | 0.263158 | 7.77975×10-2 | 0.915829 | 0.888118 | -0.334334 | 2.77111×10-2 |
2.0 | 0.972671 | 0.906026 | -4.400130 | 6.66455×10-2 | 0.968915 | 0.943982 | -7.027070 | 2.49329×10-2 |
2.5 | 0.986094 | 0.932703 | -25.95750 | 5.33910×10-2 | 0.985180 | 0.968440 | -32.40280 | 1.67394×10-2 |
3.0 | 0.993557 | 0.948165 | -97.30430 | 4.53923×10-2 | 0.992469 | 0.979829 | -104.6460 | 1.26404×10-2 |
3.5 | 0.999402 | 0.957932 | -288.1960 | 4.14707×10-2 | 0.999263 | 0.985618 | -275.6610 | 1.36449×10-2 |
4.0 | 1.000703 | 0.964556 | -728.4980 | 3.61731×10-2 | 1.001750 | 0.988867 | -631.4190 | 1.28808×10-2 |
4.5 | 0.997296 | 0.969315 | -1638.840 | 2.79813×10-2 | 0.997197 | 0.990876 | -1304.100 | 6.32101×10-3 |
5.0 | 1.002409 | 0.972892 | -3369.480 | 2.95975×10-2 | 1.003140 | 0.992229 | -2486.000 | 1.09120×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yEXACT | yMHPM | Pabs | yEHPM | yEXACT | yMHPM | Pabs | |
0.1 | 0.009336 | 0.01 | 0.001527 | 6.63272×10-4 | 0.009782 | 0.01 | 0.001150 | 2.17468×10-4 |
0.2 | 0.037202 | 0.04 | 0.011010 | 2.79761×10-3 | 0.039123 | 0.04 | 0.008892 | 8.76525×10-4 |
0.3 | 0.083691 | 0.09 | 0.034966 | 6.30802×10-3 | 0.087988 | 0.09 | 0.029411 | 2.01164×10-3 |
0.4 | 0.148502 | 0.16 | 0.079383 | 1.14971×10-2 | 0.156347 | 0.16 | 0.068718 | 3.65298×10-3 |
0.5 | 0.231266 | 0.25 | 0.149941 | 1.87333×10-2 | 0.244079 | 0.25 | 0.132726 | 5.92024×10-3 |
0.6 | 0.331497 | 0.36 | 0.251096 | 2.85027×10-2 | 0.351006 | 0.36 | 0.227269 | 8.99314×10-3 |
0.7 | 0.448290 | 0.49 | 0.391189 | 4.17099×10-2 | 0.476835 | 0.49 | 0.358125 | 1.31646×10-2 |
0.8 | 0.579872 | 0.64 | 0.572354 | 6.01272×10-2 | 0.621014 | 0.64 | 0.531020 | 1.89856×10-2 |
0.9 | 0.723061 | 0.81 | 0.800662 | 8.69390×10-2 | 0.782501 | 0.81 | 0.751641 | 2.74983×10-2 |
1.0 | 0.872643 | 1.00 | 1.081081 | 1.27356×10-1 | 0.959450 | 1.00 | 1.025641 | 4.05497×10-2 |
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.475021 | 0.474993 | 2.46866×10-9 |
0.2 | 0.450166 | 0.450166 | 0.450166 | 0.450166 | 0.450165 | 0.449984 | 8.11591×10-9 |
0.3 | 0.425557 | 0.425557 | 0.425557 | 0.425557 | 0.425552 | 0.425056 | 2.62871×10-9 |
0.4 | 0.401312 | 0.401312 | 0.401312 | 0.401312 | 0.401291 | 0.400332 | 6.47150×10-9 |
0.5 | 0.377541 | 0.377541 | 0.377541 | 0.377541 | 0.377476 | 0.375992 | 1.62494×10-9 |
0.6 | 0.354344 | 0.354344 | 0.354344 | 0.354344 | 0.354182 | 0.352276 | 6.16224×10-9 |
0.7 | 0.331812 | 0.331812 | 0.331812 | 0.331813 | 0.331462 | 0.329493 | 6.81599×10-10 |
0.8 | 0.310026 | 0.310026 | 0.310026 | 0.310028 | 0.309343 | 0.308033 | 6.05932×10-9 |
0.9 | 0.289050 | 0.289051 | 0.289051 | 0.289058 | 0.287820 | 0.288369 | 3.75236×10-9 |
1.0 | 0.268941 | 0.268941 | 0.268941 | 0.268961 | 0.266858 | 0.271073 | 4.53833×10-9 |
t | yEXACT | yEHPM | yABMA | yRPSM | yNHPM | yAPSO | Pabs |
0.1 | 0.110295 | 0.110323 | 0.110295 | 0.110295 | 0.110294 | 0.132598 | 2.79559×10-5 |
0.2 | 0.241976 | 0.241942 | 0.241976 | 0.241976 | 0.241965 | 0.275045 | 3.40483×10-5 |
0.3 | 0.395104 | 0.395107 | 0.395104 | 0.395089 | 0.395106 | 0.427310 | 2.56995×10-6 |
0.4 | 0.567812 | 0.567860 | 0.567811 | 0.56766 | 0.568115 | 0.589222 | 4.80912×10-5 |
0.5 | 0.756014 | 0.756031 | 0.756014 | 0.755134 | 0.757564 | 0.760469 | 1.68103×10-5 |
0.6 | 0.953566 | 0.953530 | 0.953565 | 0.949964 | 0.958259 | 0.940590 | 3.57596×10-5 |
0.7 | 1.152949 | 1.152943 | 1.152949 | 1.141423 | 1.163459 | 1.128970 | 5.97971×10-6 |
0.8 | 1.346364 | 1.346424 | 1.346363 | 1.315723 | 1.365240 | 1.324840 | 6.04658×10-5 |
0.9 | 1.526911 | 1.526896 | 1.526911 | 1.456545 | 1.554960 | 1.527270 | 1.51114×10-5 |
1.0 | 1.689498 | 1.689546 | 1.689498 | 1.546030 | 1.723810 | 1.735140 | 4.79364×10-5 |
t | yEXACT | yEHPM | yABMA | yRPSM | yMHPM | yAPSO | Pabs |
0.1 | 0.099668 | 0.099666 | 0.099668 | 0.099668 | 0.099668 | 0.103891 | 1.23113×10-6 |
0.2 | 0.197375 | 0.197378 | 0.197375 | 0.197375 | 0.197375 | 0.202649 | 2.55077×10-6 |
0.3 | 0.291313 | 0.291313 | 0.291313 | 0.291312 | 0.291312 | 0.295810 | 1.45452×10-7 |
0.4 | 0.379949 | 0.379946 | 0.379949 | 0.379944 | 0.379944 | 0.382931 | 2.41309×10-6 |
0.5 | 0.462117 | 0.462117 | 0.462117 | 0.462078 | 0.462078 | 0.463587 | 1.66903×10-7 |
0.6 | 0.537050 | 0.537052 | 0.537049 | 0.536857 | 0.536857 | 0.537380 | 2.46396×10-6 |
0.7 | 0.604368 | 0.604368 | 0.604368 | 0.603631 | 0.603631 | 0.603935 | 4.64338×10-7 |
0.8 | 0.664037 | 0.664034 | 0.664037 | 0.661706 | 0.661706 | 0.662905 | 2.81496×10-6 |
0.9 | 0.716298 | 0.716299 | 0.716298 | 0.709919 | 0.709919 | 0.713971 | 1.30183×10-6 |
1.0 | 0.761594 | 0.761592 | 0.761594 | 0.746032 | 0.746032 | 0.756838 | 2.22537×10-6 |
t | yEXACT | yEHPM | yABMA | yMHPM | yAPSO | Pabs |
0.1 | 0.01 | 0.01 | 0.019966 | 0.001 | 0.012254 | 2.91364×10-13 |
0.2 | 0.04 | 0.04 | 0.079449 | 0.008 | 0.044342 | 1.58054×10-13 |
0.3 | 0.09 | 0.09 | 0.177089 | 0.027 | 0.096103 | 2.13080×10-13 |
0.4 | 0.16 | 0.16 | 0.310209 | 0.064 | 0.167389 | 4.90163×10-14 |
0.5 | 0.25 | 0.25 | 0.473957 | 0.125 | 0.258064 | 2.18325×10-13 |
0.6 | 0.36 | 0.36 | 0.659526 | 0.216 | 0.368016 | 8.49321×10-14 |
0.7 | 0.49 | 0.49 | 0.850389 | 0.343 | 0.497157 | 1.76026×10-13 |
0.8 | 0.64 | 0.64 | 1.014271 | 0.512 | 0.645426 | 1.66533×10-13 |
0.9 | 0.81 | 0.81 | 1.086366 | 0.729 | 0.812795 | 2.24820×10-13 |
1.0 | 1.00 | 1.00 | 0.938449 | 1.000 | 0.999274 | 3.39728×10-14 |
n | Test Problem 1 | Test Problem 2 | Test Problem 3 | Test Problem 4 | ||||
Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | Lrms | L∞ | |
5 | 4.82×10-9 | 8.11×10-9 | 3.46×10-5 | 6.04×10-5 | 1.86×10-6 | 2.81×10-6 | 1.79×10-13 | 2.91×10-13 |
4 | 1.98×10-7 | 2.63×10-7 | 1.67×10-4 | 2.63×10-4 | 1.78×10-5 | 2.36×10-5 | 1.73×10-14 | 3.25×10-14 |
3 | 1.86×10-6 | 3.04×10-6 | 1.12×10-3 | 2.01×10-3 | 3.33×10-5 | 6.10×10-5 | 6.49×10-9 | 1.13×10-8 |
2 | 2.94×10-5 | 5.27×10-5 | 3.51×10-3 | 6.39×10-3 | 1.03×10-3 | 1.76×10-3 | 5.17×10-9 | 8.62×10-9 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.1 | 0.458464 | 0.462742 | 0.462742 | 4.27763×10-3 | 0.470459 | 0.471407 | 0.471407 | 9.4861×10-4 |
0.2 | 0.425630 | 0.433189 | 0.433189 | 7.55906×10-3 | 0.443171 | 0.444939 | 0.444939 | 1.76834×10-3 |
0.3 | 0.395869 | 0.406386 | 0.406387 | 1.05170×10-2 | 0.416934 | 0.419474 | 0.419474 | 2.54007×10-3 |
0.4 | 0.368265 | 0.381546 | 0.381546 | 1.32809×10-2 | 0.391583 | 0.394880 | 0.394880 | 3.29609×10-3 |
0.5 | 0.342440 | 0.358340 | 0.358340 | 1.58994×10-2 | 0.367092 | 0.371144 | 0.371145 | 4.05246×10-3 |
0.6 | 0.318209 | 0.336593 | 0.336596 | 1.83840×10-2 | 0.343479 | 0.348290 | 0.348291 | 4.81164×10-3 |
0.7 | 0.295471 | 0.316194 | 0.316201 | 2.07229×10-2 | 0.320785 | 0.326348 | 0.326350 | 5.56302×10-3 |
0.8 | 0.274172 | 0.297059 | 0.297079 | 2.28877×10-2 | 0.299065 | 0.305348 | 0.305353 | 6.28306×10-3 |
0.9 | 0.254283 | 0.279120 | 0.279167 | 2.48372×10-2 | 0.278377 | 0.285313 | 0.285327 | 6.93573×10-3 |
1.0 | 0.235793 | 0.262314 | 0.262417 | 2.65208×10-2 | 0.258786 | 0.266259 | 0.266294 | 7.47296×10-3 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.4 | 1.128780 | 0.772076 | 0.701892 | 3.56704×10-1 | 0.951216 | 0.629409 | 0.609030 | 3.21808×10-1 |
0.8 | 1.732609 | 1.534233 | 1.292070 | 1.98376×10-1 | 1.582695 | 1.414490 | 1.339680 | 1.68205×10-1 |
1.2 | 2.095352 | 1.965731 | -0.689702 | 1.29622×10-1 | 1.999238 | 1.962805 | 1.116260 | 3.64331×10-2 |
1.6 | 2.300624 | 2.159265 | -12.26610 | 1.41359×10-1 | 2.252742 | 2.217871 | -1.984430 | 3.48714×10-2 |
2.0 | 2.399864 | 2.248003 | -48.77540 | 1.51861×10-1 | 2.384753 | 2.319886 | -9.224490 | 6.48668×10-2 |
2.4 | 2.430901 | 2.293382 | -139.1670 | 1.37519×10-1 | 2.431978 | 2.361076 | -17.66350 | 7.09018×10-2 |
2.8 | 2.424109 | 2.319480 | -332.8130 | 1.04629×10-1 | 2.428249 | 2.379249 | -13.99210 | 4.89996×10-2 |
3.2 | 2.405166 | 2.336069 | -707.8470 | 6.90973×10-2 | 2.405453 | 2.388335 | 34.79380 | 1.71183×10-2 |
3.6 | 2.396527 | 2.347472 | -1381.040 | 4.90542×10-2 | 2.394053 | 2.393514 | 194.3750 | 5.39429×10-4 |
4.0 | 2.418283 | 2.355791 | -2519.200 | 6.24925×10-2 | 2.423411 | 2.396826 | 580.2210 | 2.65850×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yABMA | yRPSM | Pabs | yEHPM | yABMA | yRPSM | Pabs | |
0.5 | 0.571511 | 0.515900 | 0.481043 | 5.56103×10-2 | 0.497222 | 0.481137 | 0.515392 | 1.60848×10-2 |
1.0 | 0.824167 | 0.749435 | 0.731907 | 7.47319×10-2 | 0.778543 | 0.758713 | 0.675115 | 1.98303×10-2 |
1.5 | 0.933117 | 0.855319 | 0.263158 | 7.77975×10-2 | 0.915829 | 0.888118 | -0.334334 | 2.77111×10-2 |
2.0 | 0.972671 | 0.906026 | -4.400130 | 6.66455×10-2 | 0.968915 | 0.943982 | -7.027070 | 2.49329×10-2 |
2.5 | 0.986094 | 0.932703 | -25.95750 | 5.33910×10-2 | 0.985180 | 0.968440 | -32.40280 | 1.67394×10-2 |
3.0 | 0.993557 | 0.948165 | -97.30430 | 4.53923×10-2 | 0.992469 | 0.979829 | -104.6460 | 1.26404×10-2 |
3.5 | 0.999402 | 0.957932 | -288.1960 | 4.14707×10-2 | 0.999263 | 0.985618 | -275.6610 | 1.36449×10-2 |
4.0 | 1.000703 | 0.964556 | -728.4980 | 3.61731×10-2 | 1.001750 | 0.988867 | -631.4190 | 1.28808×10-2 |
4.5 | 0.997296 | 0.969315 | -1638.840 | 2.79813×10-2 | 0.997197 | 0.990876 | -1304.100 | 6.32101×10-3 |
5.0 | 1.002409 | 0.972892 | -3369.480 | 2.95975×10-2 | 1.003140 | 0.992229 | -2486.000 | 1.09120×10-2 |
t | α=0.85 | α=0.95 | ||||||
yEHPM | yEXACT | yMHPM | Pabs | yEHPM | yEXACT | yMHPM | Pabs | |
0.1 | 0.009336 | 0.01 | 0.001527 | 6.63272×10-4 | 0.009782 | 0.01 | 0.001150 | 2.17468×10-4 |
0.2 | 0.037202 | 0.04 | 0.011010 | 2.79761×10-3 | 0.039123 | 0.04 | 0.008892 | 8.76525×10-4 |
0.3 | 0.083691 | 0.09 | 0.034966 | 6.30802×10-3 | 0.087988 | 0.09 | 0.029411 | 2.01164×10-3 |
0.4 | 0.148502 | 0.16 | 0.079383 | 1.14971×10-2 | 0.156347 | 0.16 | 0.068718 | 3.65298×10-3 |
0.5 | 0.231266 | 0.25 | 0.149941 | 1.87333×10-2 | 0.244079 | 0.25 | 0.132726 | 5.92024×10-3 |
0.6 | 0.331497 | 0.36 | 0.251096 | 2.85027×10-2 | 0.351006 | 0.36 | 0.227269 | 8.99314×10-3 |
0.7 | 0.448290 | 0.49 | 0.391189 | 4.17099×10-2 | 0.476835 | 0.49 | 0.358125 | 1.31646×10-2 |
0.8 | 0.579872 | 0.64 | 0.572354 | 6.01272×10-2 | 0.621014 | 0.64 | 0.531020 | 1.89856×10-2 |
0.9 | 0.723061 | 0.81 | 0.800662 | 8.69390×10-2 | 0.782501 | 0.81 | 0.751641 | 2.74983×10-2 |
1.0 | 0.872643 | 1.00 | 1.081081 | 1.27356×10-1 | 0.959450 | 1.00 | 1.025641 | 4.05497×10-2 |