
The present study proposes a hybrid numerical technique to discuss the solution of non-linear reaction-diffusion equations with variable coefficients. The perturbation parameter was assumed to be time-dependent. The spatial domain was discretized using the cubic Hermite splines collocation method. These splines are smooth enough to interpolate the function as well as its tangent at the node points. The temporal domain was discretized using the Crank-Nicolson scheme, commonly known as the CN scheme. The cubic Hermite splines are convergent of order h4, and the CN scheme is convergent of order Δt2. The technique is found to be convergent of order O(h2(γ2εjΔt+γ0(1+ˉα)h2)+Δt2). The step size in the space direction is taken to be h, and the step size in the time direction is Δt. Stability of the proposed scheme was studied using the L2 and L∞ norms. The proposed scheme has been applied to different sets of problems and is found to be more efficient than existing schemes.
Citation: Abdul-Majeed Ayebire, Inderpreet Kaur, Dereje Alemu Alemar, Mukhdeep Singh Manshahia, Shelly Arora. A robust technique of cubic Hermite splines to study the non-linear reaction-diffusion equation with variable coefficients[J]. AIMS Mathematics, 2024, 9(4): 8192-8213. doi: 10.3934/math.2024398
[1] | Abdul-Majeed Ayebire, Saroj Sahani, Priyanka, Shelly Arora . Numerical study of soliton behavior of generalised Kuramoto-Sivashinsky type equations with Hermite splines. AIMS Mathematics, 2025, 10(2): 2098-2130. doi: 10.3934/math.2025099 |
[2] | Mohammad Tamsir, Neeraj Dhiman, Deependra Nigam, Anand Chauhan . Approximation of Caputo time-fractional diffusion equation using redefined cubic exponential B-spline collocation technique. AIMS Mathematics, 2021, 6(4): 3805-3820. doi: 10.3934/math.2021226 |
[3] | Ajmal Ali, Tayyaba Akram, Azhar Iqbal, Poom Kumam, Thana Sutthibutpong . A numerical approach for 2D time-fractional diffusion damped wave model. AIMS Mathematics, 2023, 8(4): 8249-8273. doi: 10.3934/math.2023416 |
[4] | Azhar Iqbal, Abdullah M. Alsharif, Sahar Albosaily . Numerical study of non-linear waves for one-dimensional planar, cylindrical and spherical flow using B-spline finite element method. AIMS Mathematics, 2022, 7(8): 15417-15435. doi: 10.3934/math.2022844 |
[5] | Shanshan Wang . Split-step quintic B-spline collocation methods for nonlinear Schrödinger equations. AIMS Mathematics, 2023, 8(8): 19794-19815. doi: 10.3934/math.20231009 |
[6] | Haoran Sun, Siyu Huang, Mingyang Zhou, Yilun Li, Zhifeng Weng . A numerical investigation of nonlinear Schrödinger equation using barycentric interpolation collocation method. AIMS Mathematics, 2023, 8(1): 361-381. doi: 10.3934/math.2023017 |
[7] | Hasim Khan, Mohammad Tamsir, Manoj Singh, Ahmed Hussein Msmali, Mutum Zico Meetei . Numerical approximation of the time-fractional regularized long-wave equation emerging in ion acoustic waves in plasma. AIMS Mathematics, 2025, 10(3): 5651-5670. doi: 10.3934/math.2025261 |
[8] | Asifa Tassaddiq, Muhammad Yaseen, Aatika Yousaf, Rekha Srivastava . Computational study of the convection-diffusion equation using new cubic B-spline approximations. AIMS Mathematics, 2021, 6(5): 4370-4393. doi: 10.3934/math.2021259 |
[9] | Shafeeq Rahman Thottoli, Mohammad Tamsir, Mutum Zico Meetei, Ahmed H. Msmali . Numerical investigation of nonlinear extended Fisher-Kolmogorov equation via quintic trigonometric B-spline collocation technique. AIMS Mathematics, 2024, 9(7): 17339-17358. doi: 10.3934/math.2024843 |
[10] | Jin Xie, Xiaoyan Liu, Lei Zhu, Yuqing Ma, Ke Zhang . The C3 parametric eighth-degree interpolation spline function. AIMS Mathematics, 2023, 8(6): 14623-14632. doi: 10.3934/math.2023748 |
The present study proposes a hybrid numerical technique to discuss the solution of non-linear reaction-diffusion equations with variable coefficients. The perturbation parameter was assumed to be time-dependent. The spatial domain was discretized using the cubic Hermite splines collocation method. These splines are smooth enough to interpolate the function as well as its tangent at the node points. The temporal domain was discretized using the Crank-Nicolson scheme, commonly known as the CN scheme. The cubic Hermite splines are convergent of order h4, and the CN scheme is convergent of order Δt2. The technique is found to be convergent of order O(h2(γ2εjΔt+γ0(1+ˉα)h2)+Δt2). The step size in the space direction is taken to be h, and the step size in the time direction is Δt. Stability of the proposed scheme was studied using the L2 and L∞ norms. The proposed scheme has been applied to different sets of problems and is found to be more efficient than existing schemes.
Non-linear partial differential equations are widely used by scientists and engineers to interpret various physical phenomena, such as shallow water wave behavior, viscoelastic behavior of fluids, and shock wave behavior in fluids. In different branches of science and engineering, various physical processes are described theoretically through ordinary or partial differential equations. The Lane-Emden equation [1,2], Emden Fowler equation [3,4,5], and Riccati equation [6] are examples of ordinary differential equations which are used to describe various physical, chemical, and biological phenomena, including the problems in fluid flow and engineering. Burger's equation is used to study the physical processes in aerospace engineering and fluid dynamics [7]. The non-local diffusion equation is used to describe neutron transport in a nuclear reactor [8]. The Burger-Huxley equation, Burger-Fisher equation [7,9,10], Hodgkin-Huxley equation [11,12], and Fitzough-Nagumo equation [13,14] have numerous applications in the field of physics, biology, fluid dynamics, engineering, optics, and plasma physics, among others.
Consider the one-dimensional non-linear reaction-diffusion equation:
∂u∂τ=ε(τ)∂2u∂ζ2+α(τ)f(u),Ω=I×(0,T), | (1.1) |
where f(u) is the non-linear source function, ε(τ) is the function of time whose value lies between 0 and 1, and I=(a,b). The initial condition is taken to be continuous whereas boundary conditions are assumed to be of the Dirichlet type. The initial condition is defined as:
u(ζ,0)=g(ζ). | (1.2) |
Both the homogeneous Dirichlet type and non-homogeneous Dirichlet type boundary conditions have been considered and are defined as:
u(a,τ)=K1(τ)andu(b,τ)=K2(τ), | (1.3) |
where, K1 and K2 are continuous functions of τ.
The simplified form of Eq (1.1) for ε(τ)=1 was first proposed by Hodgkin and Huxley to explain the propagation of action potentials in the squid giant axiom. The generalized Burger's, generalized Burger's-Huxley, and generalized Burger's-Fisher equations are the special cases of Eq (1.1) with non-linear advection term and source term. These equations have numerous applications in reaction mechanics and diffusion transport in fluids and semi-solid particles.
A variety of numerical techniques have been developed by different investigators to study the solution behavior of non-linear partial differential equations [1,5,6,15,16,17,18]. Lagurre polynomials were utilized in [1] to study the behavior of the Lane-Emden equation whereas [2] has used the Bessel collocation approach. The Emden-Fowler equation was solved using artificial neural networks by [5]. Orthogonal collocation on finite elements was used by [19] to discretize the non-linear heat conduction equation. The coupled Burger and coupled Burger-Huxley equations were discretized using the Bessel collocation method in [7]. The two phase non-linear reaction diffusion model was discretized using cubic Hermite splines in [20], whereas singularly perturbed reaction diffusion equations were solved by orthogonal cubic splines of order 3 and order 5 in [21]. The non-local Nagumo type equation was solved by using a finite volume scheme in [22]. The non-local fourth-order sub-diffusion equation was solved using orthogonal spline collocation in [18], and the non-local heat model was discretized using the sinc collocation method in [23].
In the present study, cubic Hermite splines with the Crank-Nicolson scheme is used to discretize the non-linear reaction-diffusion equation. The cubic Hermite splines are third-order Hermite splines that interpolate the function as well as its first-order derivative at node points. These splines have an advantage over the cubic B-splines as these splines do not involve any fictitious node points and thus are easily adapted to a computer program. The combination of cubic Hermite splines with the CN scheme in the time direction gives an advantage to discretize the equation in both the the space and, the time directions. It reduces the stiffness not only in space direction but also in time direction.
The present paper has been weaved into different sections and sub-sections, from the introduction to the conclusions. The CN scheme is discussed in Section 2 and the cubic Hermite splines are explained in Section 3. The stability and convergence of the proposed technique is discussed in Section 4, whereas Section 5 is composed of the numerical applications. The conclusions of the study are presented in Section 6.
A variety of numerical techniques for solving two point boundary value problems in the time direction have been developed. The Crank-Nicolson (CN) and Alternate Direction Implicit scheme are two such examples. The Alternate Direction Implicit scheme, commonly known as ADI scheme, was introduced by Peaceman and Rachford to study second-order parabolic boundary value problems over the rectangular domain. ADI methods are capable of solving multidimensional problems into collection of one-dimensional problems [24]. Other than the simple ADI scheme or the Peaceman-Rachford ADI scheme, [25] discussed certain other classifications of ADI schemes, such as the Douglas scheme, Craig-Sneyd scheme, modified Craig-Sneyd scheme, Hundsdorfer-Verwer scheme, etc. Crank-Nicolson scheme, or CN scheme, is also an efficient time integration scheme to solve one-dimensional reaction-diffusion problems. In the present study, the CN scheme as proposed by [26] operates on the mean of a function over the interval [τjτj+1] with a uniform distribution of points [16,27].
Define the partition πτ:0=τ0<τ1<…<τM=T with Δt=τj+1−τj. Consider the general non-linear reaction-diffusion equation defined by Eq (1.1)
u(ζ,τj+1)−u(ζ,τj)Δt=ε(τj)(uj+1ζζ+ujζζ2)+α(τj)f(uj+1)+f(uj)2,u(ζ,τ0)=g(ζ). | (2.1) |
For convenience, we write u(ζ,τj+1)=uj+1.
Lemma 1. [16,28] If |∂iu∂τi|≤C, i=0,1, and (ζ,τ)∈ˉΩ, then the local truncation error e is given by
e=CΔt3, |
and the global truncation error E is given by
E=CΔt2. |
Theorem 1. [16,28] Let U(ζ, τ) be a function such that £U = 0 and U≤C. Then, Uτ is also bounded by some constant C, where C is an arbitrary constant.
For the implementation of the technique on non-linear singularly perturbed equations, let uj+1 be the approximating polynomial function at time step τj+1. Also, let Δt=T/M, M>0. Then, τj=jΔt for j=0,1,2,…,M. The semi-discretized form of Eq (2.1) follows as:
uj+1−ujΔt=εjuj+12ζζ+αjfj+12, | (2.2) |
with corresponding boundary conditions are given by Eqs (1.2) and (1.3), respectively. For any function F,
Fj+12=Fj+1+Fj2. | (2.3) |
Throughout this paper, C denotes the generic constant.
After the implementation of the CN scheme, Eq (1.1) reduces to the semi-discretized form as follows:
u(ζ,τj+1)−u(ζ,τj)Δt=εj(uj+1ζζ+ujζζ2)+αjf(uj+1)+f(uj)2,j=0,1,2,…,M. | (2.4) |
The non-linear function term f(uj+1) is quasi-linearized using the formula given by [29]:
f(uj+1)=f(uj)+(uj+1−uj)dfduj+O(Δt2), | (2.5) |
where, dfduj represents the derivative of f(u) at τj. After substituting Eq (2.5) in Eq (2.4) and rearranging the terms, one gets the iterated scheme
uj+1−εjΔt2uj+1ζζ−αjΔtfuj2uj+1=uj+εjΔt2ujζζ−αjΔtfuj2uj+αjΔtf(uj), | (2.6) |
where fju represents dfduj for convenience. The iterative form of Eq (2.6) can be expressed in operator form at the (j+1)th time step in the following manner:
L(uj+1)=Gτj, | (2.7) |
where, L(uj+1)=uj+1−εjΔt2uj+1ζζ−αjΔtfuj2uj+1 and Gτj=uj+εjΔt2ujζζ−αjΔtfuj2uj+αjΔtf(uj). Since there is no singular term on the right-hand side, Eq (2.7) is convergent, which leads to the solution of Eq (1.1).
Orthogonal splines represent the piecewise orthogonal polynomials used to interpolate the function at node points. Hermite and Lagrange's interpolating polynomials are such orthogonal splines that are often used to interpolate functions. Hermite interpolating polynomials of order 'k' are considered as an extension of kth-order Lagrangian interpolating polynomials [7,9,10,16,20,21,27]. A continuity condition is imposed at node points in the Lagrangian interpolating polynomials as these polynomials can only interpolate the function at node points. Hermite interpolating polynomials of order '2k+1' interpolate the function as well as its kth order derivative at node points. This feature of Hermite interpolating polynomials makes it superior over Lagrangian interpolating polynomials. Hence, th kth-order Hermite polynomial in ζ is a polynomial of order 2k+1, and, therefore, cubic Hermite interpolating polynomials are a particular case of general Hermite interpolating polynomials for k=1. This consists of two node points and two tangents at these node points. The explanation of cubic Hermite splines is given hereunder.
Consider an interval I=(a,b), and let π be the partition of (a,b) such that:
πζ:a=ζ0≤ζ1≤ζ2≤…≤ζn=b. |
Let P3 be the set of all polynomials of degree less than or equal to 3 defined on [ζi−1,ζi]. Let v be a continuously differentiable function defined on ˉI such that v is the cubic Hermite approximation of u on I. Let Mμ be the space of all continuously differentiable functions defined on [ζi−1,ζi] such that
Mμ={vϵC1[a,b]|v∈P3on[ζi−1,ζi],i=1,2,…,n},M0μ={vϵMμ|v(a)=K1 and v(b)=K2}. | (3.1) |
If boundary conditions are of the homogeneous type, then the dimension of Mμ=2n, and if the boundary conditions are of the non-homogeneous type, then the dimension of Mμ=2n+2 [24].
The cubic Hermite polynomials Pμ(ζ) and Qμ(ζ) for k=1 are defined as
Pμ(ζ)={3(ζ−ζμ−1ζμ−ζμ−1)2−2(ζ−ζμ−1ζμ−ζμ−1)3,ζμ−1≤ζ≤ζμ,3(ζμ+1−ζζμ+1−ζμ)2−2(ζμ+1−ζζμ+1−ζμ)3,ζμ≤ζ≤ζμ+1,0,elsewhere, | (3.2) |
Qμ(ζ)={−(ζ−ζμ−1)2ζμ−ζμ−1+(ζ−ζμ−1)3(ζμ−ζμ−1)2,ζμ−1≤ζ≤ζμ,(ζμ+1−ζ)2ζμ+1−ζμ−(ζμ+1−ζ)3(ζμ+1−ζμ)2,ζμ≤ζ≤ζμ+1,0,elsewhere. | (3.3) |
These piecewise cubics are designed such that the following identities hold: Pμ(ζi)=δμi; P′μ(ζi)=0; Qμ(ζi)=0; Q′μ(ζi)=δμi.
To apply collocation in the interval [ζi−1,ζi], a new variable ξ is introduced in such a way that, as ζ varies from ζi−1 to ζi, the variable ξ varies from 0 to 1. It converts the polynomials defined in Eqs (3.2) and (3.3) into the following form:
H1(ξ)=1−3ξ2+2ξ3,H2(ξ)=hξ2(3−2ξ),H3(ξ)=ξ(ξ−1)2,H4(ξ)=hξ2(ξ−1). |
These four cubic Hermite interpolating polynomials form the basis of Mμ for k=1. These interpolating polynomials are bounded and vary from −0.5 to 1. The graphical representation of these polynomials is shown in Figure 1. In Table 1, the values of these polynomials and their respective derivatives are shown at boundary points 0 and 1.
ξ=0 | ξ=1 | ξ=0 | ξ=1 | ||
H1 | 1 | 0 | H′1 | 0 | 0 |
H2 | 0 | 1 | H′2 | 0 | 0 |
H3 | 0 | 0 | H′3 | 1 | 0 |
H4 | 0 | 0 | H′4 | 0 | 1 |
Let G be the space of all those points which are zeros of the shifted Legendre polynomials of order n. Let ηj be the node points of the Gauss-Legendre quadrature formula defined on [0,1], where j = 1,2, …, n−1. The collocation points are defined by
ξ(μ−1)(n−1)=ζμ−1+hμηj,μ=1,2,…,mandj=1,2,…,n−1, |
where hμ=ζμ−ζμ−1. Runge's divergence formula also states that non-uniform collocation points give less error compared to uniform collocation points [15,19].
Theorem 2. [30] The necessary and sufficient condition for a simple set Bn(ξ) of real polynomials to be orthogonal w.r.t. the weight function ω(ξ) on a≤ξ≤b is
∫baω(ξ)ξmPn(ξ)dξ=0,m=0,1,2,3,...,(n−1), |
where, Pn(ξ)∈Bn(ξ). The roots of these polynomials are always real and distinct and lie in the interval a≤ξ≤b.
The collocation points are chosen to be the zeros of the Legendre polynomials as these polynomials have the tendency to reduce the error at corners as well as on averages [19]. The zeros of the Legendre polynomials are calculated using the following recurrence relation given in [31]:
Pi(ξ)=(ξ−0.5)Pi−1(ξ)−(i−1)24(2i−3)(2i−1)Pi−2(ξ),i=1,2,…,n+1. |
The details of these polynomials are given elsewhere [19,31].
An application of orthogonal collocation using cubic Hermite splines as basis polynomials is applied with consideration of the approximate solution of a semi-discrete scheme within each sub-element of the partition πζ. A linear combination of four stencils with time dependent coefficients is proposed as follows:
u(ξ,τj)=4∑i=1Hi(ξ)σji, | (3.4) |
where σji are the unknown coefficients to be determined. After substituting Eq (3.4) into Eq (2.6), the following system of equations is obtained:
4∑i=1(Hi(ξ)−εjΔt2h2H″i(ξ)−αjΔtfuj2Hi(ξ))σj+1=4∑i=1(Hi(ξ)+εjΔt2h2H″i(ξ)−αjΔtfuj2Hi(ξ))σj+αjΔtf(4∑i=1Hi(ξ)σj). | (3.5) |
At the kth collocation point, Eq (3.5) reduces to
4∑i=1(Hki−εjΔt2h2H″ki−αjΔtfuj2Hki)σj+1=4∑i=1(Hki+εjΔt2h2H″ki−αjΔtfuj2Hki)σj+αjΔtfjk, | (3.6) |
where Hki and H″ki are the Hermite spline and the second-order derivative of the Hermite spline at the kth collocation point, respectively, and fjk=f(∑4i=1Hkiσj). This system of equations reduces to the matrix form
Qρj+1=Hρj+Fj, | (3.7) |
where Q and H are almost quad-diagonal dominant matrices of order 2n×2n, where n denotes the number of sub-elements in πζ. ρn=[σn1,σn2,…,σn2n]′; Fn=[fn(ξ1),fn(ξ2),…,fn(ξ2n)]′. The resulting system of matrix equations can be solved by any iterative technique.
In combined form, Eq (3.7) can be written as
Qρ=H, | (3.8) |
where H=Hρj+Fj.
The implementation of any numerical technique depends upon its stability. A method is said to be stable if its jth iteration is bounded [9,10]. The function f(u) is bounded over the domain I×[0,T] and satisfies the Lipschitz condition. dfdu is bounded over the desired domain such that ∣f(u)∣≤K1 and ∣dfdu∣≤K.
Define the inner product <u,v>=∫Ωu(ζ)v(ζ)dζ, where <u,u>=∣∣u∣∣2.
Lemma 2. Let f(u) be defined on the domain I×[0,T] such that f(u) satisfies Lagrange's theorem. If the derivative dfdu is bounded over the domain I×[0,T] s.t. ∣dfdu∣≤K, then f(u) satisfies the Lipschitz condition
∣f(u1)−f(u2)∣≤K∣u1−u2∣,∀u1,u2 ϵ I×[0,T]. | (4.1) |
Theorem 3. Cauchy-Schwarz inequality: Let u and v be any two functions defined on R. Then,
∣<u,v>∣≤∣∣u∣∣∣∣v∣∣. |
Lemma 3. [28] Let u(ζ,0) be the positive continuous initial solution of the non-linear reaction-diffusion equation with 0≤u(ζ,0)≤1. Then, the bounds for solution u(ζ,t) remains constant with time.
Throughout this paper, C denotes the generic constant.
Lemma 4. [28] Let u(ζ,τ) be the solution of Eq (1.1) in ˉΩ. Then, the bounds on u(ζ,τ) is given by
|u(ζ,τ)|≤C, ∀(ζ,τ)∈ˉΩ. | (4.2) |
Lemma 5. [28] Maximum Principle: Let u(ζ,τ) be the solution of Eq (1.1) in ˉΩ, u(ζ,τ)≥0 on ∂Ω, and Lu(ζ,τ)≥0 on Ω. Then, u(ζ,τ)≥0 on ˉΩ.
Lemma 6. [28] Let u(ζ,τ) be the solution of Eq (1.1) in ˉΩ. Then, the bound on the derivative of u w.r.t. τ is given by
|∂u∂τ|≤C, ∀(ζ,τ)∈ˉΩ. | (4.3) |
Lemma 7. [28] Let u(ζ,τ) be the solution of Eq (1.1) in ˉΩ. Then, the bound on the derivative of u w.r.t. ζ is given by
|∂iu∂ζi|≤C, ∀(ζ,τ)∈ˉΩ, i=0,1,2. | (4.4) |
Lemma 8. [28] Let u(ζ,τ) be the solution of Eq (1.1) in ˉΩ. Then, the bound on the derivatives of u are given by
|∂i+ju∂ζi∂τj|≤C, ∀(ζ,τ)∈ˉΩ, 0≤i+j≤3. | (4.5) |
Theorem 4. Let uj be the approximate solution of Eq (1.1) at τj. Then, the following inequality holds:
∣∣uj+1∣∣≤∣∣uj∣∣+η2η1∣∣ujζζ∣∣+αjΔtη1∣∣fj∣∣. | (4.6) |
Proof. Now, using the bound for dfdu and redefining Eq (2.6), the following equation is obtained:
(1−αjΔtK2)uj+1−εjΔt2uj+1ζζ=(1−αjΔtK2)uj+εjΔt2ujζζ+αjΔtfj, | (4.7) |
where η1=(1−αjΔtK2) and η2=εjΔt2.
Taking the inner product of both sides of uj+1
η1<uj+1,uj+1>−η2<uj+1ζζ,uj+1> ≤ η1<uj,uj+1>+η2<ujζζ,uj+1>+αjΔt<fj,uj+1>. | (4.8) |
Now, using the fact that −<uj+1ζζ,uj+1> ≥ 0, and from Cauchy-Schwarz inequality that ∣<uj+1,uj>∣ ≤ ∣∣uj+1∣∣ ∣∣uj∣∣, Eq (4.8) takes the following form:
η1∣∣uj+1∣∣2≤η1∣∣uj+1∣∣∣∣uj∣∣+η2∣∣ujζζ∣∣∣∣uj+1∣∣+αjΔt∣∣fj∣∣∣∣uj+1∣∣. | (4.9) |
Dividing both sides by η1∣∣uj+1∣∣,
∣∣uj+1∣∣≤∣∣uj∣∣+η2η1∣∣ujζζ∣∣+αjΔtη1∣∣fj∣∣. | (4.10) |
Using Lemmas 2, 3, and 8 and the principle of induction, it is clear that right-hand side is bounded. This implies ∣∣uj+1∣∣ is bounded, and hence the result.
Lemma 9. [32] If any Hermite interpolation polynomial based on at most two distinct points has a unique solution in Ω, then any Hermite interpolation polynomial has a unique solution in Ω. This implies that any Hermite interpolation polynomial has a unique solution in the polynomial space Pn of degree n on I.
Lemma 10. [32] The univariate Hermite interpolating polynomials are regular for any set of node points. Moreover, Hermite interpolating polynomials are regular for any nodal set of points as well as for any choice of derivatives to be interpolated.
Lemma 11. Let ˉH be the space of all Hermite interpolating polynomials of order 3 defined on [0 1]. Then, the Hermite splines of order 3 are bounded with upper bound unity. Moreover, ∑4i=1∣Hi(ξ)∣≤4 for all 0≤ξ≤1.
Lemma 12. [27,33,34] Let H ∈H4Δξ(ξ) be the piecewise Hermite spline of degree 3 over the subinterval [ξi,ξi+1] approximating U ∈C4[a, b]. Then, ∣∣H(r)−U(r)∣∣∞≤Cγrh4−r,r=0,1,2,3, where H(r) represents the rth-order derivative of H. The values of γr are given in [27,33,34].
Theorem 5. Let u(ζ) be the exact solution of Eq (1.1) and ˉu(ζ) be the approximate solution of Eq (1.1) in the space H3 of cubic Hermite interpolating polynomials of order 3 such that ˉu(ζ)∈C4[a,b]. Then, the uniform error estimate is given by
∣∣u(ζ)−ˉu(ζ)∣∣∞≤Ch2(γ2εjΔt+γ0(1+ˉα)h2). | (4.11) |
Proof. Let H(ζ) be the unique cubic Hermite spline interpolate of u(ζ), and πζ=a=ζ0≤ζ1≤,...,≤ζn=b be the equi-spaced partition of [a, b] with uniform step-size h.
Consider
||L(u(ζk))−L(H(ζk))||∞=||εjΔt2(u″(ζk)−H″(ζk))+(1−αjfujΔt2)(u(ζk)−H(ζk))||∞≤|εjΔt2| ||(u″(ζk)−H″(ζk))||∞+ˉα2||(u(ζk)−H(ζk))||∞, | (4.12) |
where ˉα/2 = |1−αjfujΔt2|.
Using Lemma 12, the following inequality is obtained:
||L(u(ζk))−L(H(ζk))||∞≤Ch2(γ2εjΔt+γ0ˉαh2), | (4.13) |
where C is the generic constant and the values of γ0 and γ2 are given in [27,33,34].
Now, consider
||L(ˉu(ζk))−L(H(ζk))||∞=||Gτj(ζk)−L(H(ζk))||∞≤||L(u(ζk))−L(H(ζk))||∞≤Ch2(γ2εjΔt+γ0ˉαh2). | (4.14) |
Now, we use the fact that L(ˉu(ζk))=ˉGτj(ζk) and Qρ=H.
This implies that Q(ρ−ˉρ)=H−ˉH, and using Eq (4.14) we get
||H−ˉH||∞≤max||Gτj(ζk)−ˉGτj(ζk)||∞≤max||L(ˉu(ζk))−L(H(ζk))||≤Ch2(γ2εjΔt+γ0ˉαh2). | (4.15) |
Also, Q, the matrix of collocation coefficients, is bounded. Therefore, ||Q||≤C and
⇒||H−ˉH||∞≤Ch2(γ2εjΔt+γ0(1+αjη2Δt)h2). | (4.16) |
Now, ˉu(ζ)−H(ζ)=∑4i=1Hi(ζ)(σi−ˉσi), and we get
⇒||ˉu(ζ)−H(ζ)||∞=||4∑i=1Hi(ζ)(σi−ˉσi)||∞≤4∑i=1|Hi(ζ)|||(σi−ˉσi)||∞≤Ch2(γ2εjΔt+γ0(1+αjηΔt)h2). | (4.17) |
According to [33],
||u(ζ)−H(ζ)||∞=Cγ0h4. | (4.18) |
Now, using the triangle inequality,
||u(ζ)−ˉu(ζ)||∞≤||u(ζ)−H(ζ)||∞+||ˉu(ζ)−H(ζ)||∞≤Cγ0h4+Ch2(γ2εjΔt+γ0ˉαh2)≤Ch2(γ2εjΔt+γ0(1+ˉα)h2). | (4.19) |
Theorem 6. Let u(ζ,τ) represent the exact solution and ˉu(ζ,τ) be the approximate solution. Then, the uniform error estimate is given by
||u−ˉu||∞≤C(h2(γ2εjΔt+γ0(1+ˉα)h2)+Δt2). |
Proof. The convergence of the proposed method in the time direction is O(Δt2) as given in Lemma 1, and from Theorem 5 the error estimate in the space direction is Ch2(γ2εjΔt+γ0(1+ˉα)h2). Therefore, combining the results from Lemma 1 and Theorem 5, the following error estimate is obtained:
||u−ˉu||∞≤C(h2(γ2εjΔt+γ0(1+ˉα)h2)+Δt2). | (4.20) |
The error analysis of any technique is incomplete if it does not consider the Eucledian and supremum norms, also known as the L2-norm and L∞-norm, respectively. These norms are defined as
∣∣u−ˉu∣∣2=√n∑i=0u(ζi,τ)2−ˉu(ζi,τ)2. |
∣∣u−ˉu∣∣∞=max.∣u(ζi,τ)−ˉu(ζi,τ)∣, i=0,1,2,…,n. |
Problem 1. Consider the Fitzough-Nagumo equation
uτ=uζζ+u(1−u)(u−μ),∀(ζ,τ)∈(−10,10)×(0,6]. | (5.1) |
The boundary as well as initial conditions are given below:
u(−10,τ)=0.5+0.5×tanh(12√2(−10−2μ−1√2τ)), |
u(10,τ)=0.5+0.5×tanh(12√2(10−2μ−1√2τ)), |
u(ζ,0)=0.5+0.5×tanh(ζ2√2). |
Problem 1 was solved using cubic Hermite splines with a=−10 and b=10. The numerical findings of the L2-norm and the L∞-norm have been calculated for μ=0.25 and μ=−0.5, which is presented in Table 2. The norms show the stability and the efficiency of the technique. In Tables 3 and 4, the absolute error is shown for μ=0.25 and μ=−0.5 for varying number of node points. It is clear from these tables that the absolute error varies from order 10−6 to 10−4. In Table 5, the numerical values obtained from present technique are compared with the values given in [35]. The L∞-norm is found to be better than the values given in the literature, and the L2-norm is found to be on par with the L2-norm given in the literature, and at some points even better. The graphical representation of the numerical solutions is visualized in Figures 2–5 for varying values of μ. It is observed from these figures that the numerical values are bounded and lie between 0 and 1. From these tables and figures, it is observed that these quantitative results show the efficiency and accuracy of the method.
μ=0.25 | μ=−0.5 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 6.53005×10−6 | 6.10156×10−5 | 2.45277×10−5 | 2.42908×10−4 |
0.6 | 1.46464×10−5 | 1.91197×10−4 | 2.39772×10−5 | 2.23646×10−4 |
0.8 | 1.88315×10−5 | 2.46652×10−4 | 2.40393×10−5 | 2.18534×10−4 |
1 | 2.10045×10−5 | 3.66614×10−4 | 2.37522×10−5 | 2.12365×10−4 |
3 | 3.12686×10−5 | 4.85899×10−4 | 3.86837×10−5 | 5.53539×10−4 |
6 | 3.51878×10−5 | 5.44749×10−4 | 5.13096×10−4 | 9.43100×10−3 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 2.420467×10−2 | 3.032328×10−3 | 4.816493×10−5 | 9.623057×10−5 | 1.519802×10−5 |
−6 | 3.805404×10−4 | 1.942463×10−5 | 4.888632×10−6 | 4.283628×10−6 | 4.157908×10−6 |
−4 | 7.190073×10−5 | 1.778285×10−5 | 1.505554×10−5 | 1.498958×10−5 | 1.498550×10−5 |
−2 | 3.583102×10−4 | 1.537630×10−5 | 4.005334×10−5 | 4.066483×10−5 | 4.070273×10−5 |
2 | 3.518331×10−4 | 5.601424×10−5 | 3.817221×10−5 | 3.772792×10−5 | 3.770037×10−5 |
4 | 2.669094×10−5 | 9.245578×10−6 | 1.320194×10−5 | 1.330002×10−5 | 1.330608×10−5 |
6 | 7.101890×10−4 | 3.316729×10−5 | 5.893619×10−6 | 4.357652×10−6 | 3.636630×10−6 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.46995×10−2 | 4.31885×10−3 | 1.18035×10−4 | 1.37984×10−4 | 2.41337×10−5 |
−6 | 5.42797×10−4 | 4.09318×10−5 | 2.07135×10−5 | 1.98629×10−5 | 1.96841×10−5 |
−4 | 6.67514×10−5 | 6.80546×10−5 | 6.95134×10−5 | 6.95535×10−5 | 6.95561×10−5 |
−2 | 2.85835×10−4 | 1.45990×10−4 | 1.76459×10−4 | 1.77204×10−4 | 1.77250×10−4 |
2 | 2.95911×10−4 | 1.38367×10−4 | 1.30731×10−4 | 1.30538×10−4 | 1.30526×10−4 |
4 | 1.37342×10−5 | 3.88842×10−5 | 4.33501×10−5 | 4.34606×10−5 | 4.34674×10−5 |
6 | 6.26818×10−4 | 3.35472×10−5 | 1.32485×10−5 | 1.21421×10−5 | 1.16359×10−5 |
8 | 3.16242×10−2 | 5.11491×10−3 | 2.48175×10−4 | 8.23161×10−5 | 4.64385×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞−norm | L∞-norm |
0.2 | 2.3012×10−6 | 1.0808×10−6 | 4.7416×10−5 | 1.9189×10−5 |
0.5 | 5.5695×10−6 | 3.5274×10−6 | 1.2312×10−4 | 5.6503×10−5 |
1 | 1.1864×10−5 | 1.0948×10−5 | 2.6261×10−4 | 1.9484×10−4 |
1.5 | 1.9400×10−5 | 3.6190×10−5 | 4.2096×10−4 | 5.2431×10−4 |
2 | 2.8162×10−5 | 2.7306×10−5 | 5.9999×10−4 | 5.3048×10−4 |
3 | 4.9735×10−5 | 6.9690×10−5 | 1.0324×10−3 | 9.9714×10−4 |
5 | 1.1395×10−4 | 1.1158×10−4 | 2.3020×10−3 | 1.9205×10−3 |
Problem 2. Consider the variable order Fitzough-Nagumo equation
uτ=cosτuζζ−cosτuζ+2cosτ(u(1−u)(u−μ)),∀(ζ,τ)∈(−10,10)×(0,1],u(−10,τ)=0.5μ+0.5μtanh(0.5μ(−10−(3−μ)sinτ)),u(10,τ)=0.5μ+0.5μtanh(0.5μ(10−(3−μ)sinτ)). | (5.2) |
Initially,
u(ζ,0)=0.5μ+0.5μtanh(0.5μζ). |
Problem 2 was solved using cubic Hermite splines as basis functions. To check the stability and the efficiency of the technique, the numerical findings of the L2-norm and L∞-norm were calculated for μ=0.25 and μ=0.05. The numerical values of the L2-norm and the L∞-norm are presented in Table 6. It is observed that the L2-norm varies from 10−6 to 10−5, whereas the L∞-norm varies from 10−5 to 10−3. In Tables 7 and 8, the absolute error is shown for μ=0.25 and μ=0.05 for varying number of node points. It is observed that the absolute error stabilizes after 100 node points for μ=0.25 and μ=0.05. In Table 9, the L2-norm and L∞-norm are compared to the values given in [35]. Both the L2-norm and L∞-norm are found to be better than the values given in the literature. It is also observed that the absolute error is reduced for small values of μ. The graphical representation of the numerical solution is presented in Figures 6–9 for different values of μ. It is clear from these figures that u is bounded and lies between 0 and 0.25.
μ=0.25 | μ=0.05 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 7.99141×10−5 | 5.58252×10−4 | 6.09782×10−6 | 7.95415×10−5 |
0.4 | 7.80824×10−5 | 7.43978×10−4 | 5.60292×10−6 | 8.00388×10−5 |
0.6 | 7.89700×10−5 | 7.70075×10−4 | 1.19943×10−5 | 1.69741×10−4 |
0.8 | 8.83713×10−5 | 7.05602×10−4 | 1.71748×10−5 | 3.17626×10−4 |
1 | 9.50768×10−5 | 1.42997×10−3 | 1.92533×10−5 | 2.93372×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.56960×10−1 | 8.40339×10−2 | 1.84735×10−2 | 1.01369×10−3 | 4.19861×10−4 |
−6 | 2.71322×10−2 | 1.76316×10−2 | 6.73960×10−3 | 4.46380×10−4 | 3.73789×10−4 |
−4 | 2.11950×10−2 | 1.10186×10−3 | 1.70975×10−3 | 2.62175×10−4 | 2.59010×10−4 |
−2 | 2.46704×10−2 | 2.77489×10−3 | 2.77545×10−3 | 2.23101×10−4 | 2.23101×10−4 |
2 | 4.97110×10−3 | 4.97179×10−3 | 1.06710×10−3 | 1.06710×10−3 | 4.38632×10−4 |
4 | 5.34227×10−3 | 5.34505×10−3 | 1.77438×10−4 | 1.77435×10−4 | 1.77435×10−4 |
6 | 4.94204×10−3 | 4.96493×10−3 | 4.16657×10−4 | 4.16660×10−4 | 4.16660×10−4 |
8 | 3.40494×10−3 | 4.05882×10−3 | 6.38229×10−4 | 6.40161×10−4 | 1.80930×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.71171×10−1 | 1.18063×10−2 | 2.34982×10−3 | 1.69686×10−3 | 9.16820×10−5 |
−6 | 1.15766×10−2 | 1.73286×10−3 | 1.04057×10−3 | 4.97002×10−4 | 2.65137×10−5 |
−4 | 1.49387×10−3 | 5.60174×10−4 | 4.02871×10−5 | 4.88264×10−5 | 1.61169×10−5 |
−2 | 1.16173×10−3 | 7.66383×10−5 | 7.72591×10−5 | 7.72621×10−5 | 1.53091×10−5 |
2 | 8.10909×10−5 | 8.08199×10−5 | 8.08200×10−5 | 8.08200×10−5 | 1.84049×10−5 |
4 | 7.97483×10−5 | 8.20047×10−5 | 8.20047×10−5 | 4.23380×10−5 | 4.23380×10−5 |
6 | 3.48375×10−5 | 8.27213×10−5 | 8.27595×10−5 | 4.04893×10−5 | 2.11177×10−5 |
8 | 1.33879×10−3 | 5.81959×10−5 | 8.22777×10−5 | 3.90138×10−5 | 2.23853×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞-norm | L∞-norm |
0.2 | 1.3122×10−6 | 2.1142×10−6 | 1.2350×10−5 | 3.0019×10−5 |
0.5 | 6.0995×10−6 | 4.8016×10−6 | 5.1986×10−5 | 8.6760×10−5 |
1 | 2.1213×10−5 | 2.1057×10−5 | 6.3283×10−4 | 3.7204×10−4 |
1.5 | 3.2340×10−5 | 2.0231×10−5 | 8.5383×10−4 | 3.3954×10−4 |
Problem 3. Consider the Newell-Whitehead-Segel equation
uτ=εuζζ+u−u3,∀(ζ,τ)∈(0,1)×(0,2]. | (5.3) |
The boundary conditions can be taken as
u(0,τ)=−0.5+0.5tanh(−0.75τ), |
u(1,τ)=−0.5+0.5tanh(0.3536−0.75τ). |
The initial condition is
u(ζ,0)=−0.5+0.5tanh(0.3536ζ). |
The above problem is solved using cubic Hermite splines. The absolute error obtained from the cubic Hermite splines is compared with values given in [36] at different time periods in Tables 10–14. It is observed that the values obtained from the cubic Hermite splines are better than those given in the literature. In Table 15, the absolute error is shown for ε=1 and τ=0.5 for different node points. It is observed that the absolute error stabilizes after 200 node points. The graphical representation of the numerical solution is shown in Figures 10 and 11 for different values of τ. The 2D and 3D graphs show that the absolute value of u is less than 1.
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 1.16440×10−11 | 2.92000×10−9 | 9.47831×10−10 |
−15 | 1.37206×10−8 | 3.43787×10−6 | 3.07688×10−9 |
25 | 4.83920×10−10 | 3.38887×10−9 | 8.42926×10−10 |
30 | 1.40960×10−11 | 9.72520×10−11 | 4.76281×10−9 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 8.82590×10−11 | 5.43320×10−9 | 5.3765×10−10 |
−15 | 1.03995×10−7 | 6.39620×10−6 | 2.8219×10−9 |
25 | 1.98915×10−9 | 7.32706×10−9 | 9.7918×10−10 |
30 | 5.79430×10−11 | 2.10500×10−10 | 1.2574×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 2.86520×10−10 | 5.43180×10−9 | 6.26624×10−10 |
−15 | 3.37614×10−7 | 6.39128×10−6 | 2.31068×10−9 |
25 | 4.60456×10−9 | 7.32168×10−9 | 1.04420×10−9 |
30 | 1.34120×10−10 | 2.08880×10−10 | 2.55231×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 6.55230×10−10 | 5.43041×10−9 | 3.48943×10−9 |
−15 | 7.72060×10−7 | 6.38636×10−6 | 2.05489×10−9 |
25 | 8.43340×10−9 | 7.31631×10−9 | 9.50172×10−9 |
30 | 2.45660×10−10 | 2.07260×10−10 | 4.66049×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 1.45683×10−6 | 5.42902×10−9 | 8.09276×10−9 |
−15 | 1.23638×10−9 | 6.38144×10−6 | 1.57301×10−9 |
25 | 1.35958×10−8 | 7.31090×10−10 | 3.57346×10−8 |
30 | 3.96040×10−10 | 2.05640×10−10 | 7.88424×10−8 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−25 | 4.124052×10−2 | 3.071154×10−6 | 3.087868×10−8 | 9.088583×10−9 | 8.092759×10−9 |
−15 | 1.623276×10−3 | 8.594511×10−6 | 6.143898×10−6 | 2.501633×10−8 | 1.573009×10−9 |
15 | 6.569221×10−4 | 3.970126×10−4 | 5.117331×10−5 | 6.449049×10−8 | 6.246470×10−8 |
25 | 5.053591×10−4 | 4.519585×10−6 | 4.490637×10−7 | 5.625901×10−8 | 3.573463×10−8 |
Orthogonal collocation with cubic Hermite splines was applied on non-linear reaction-diffusion type equations. Two equations of Fitzough-Nagumo type and one of Newell-Whitehead-Segel type were solved numerically using cubic Hermite splines with the CN scheme. It is observed that the proposed algorithm is applicable to non-linear reaction-diffusion equations with both variable and constant coefficients. The absolute error is found to be of order 10−4, which proves the efficiency of the proposed algorithm. The proposed technique can also be applied on higher-dimensional linear as well as non-linear problems. The proposed technique can also be applied to fractional-order partial differential equations.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Dr. Shelly Arora is thankful to DST-SERB for providing grant via SPG/2022/001269 and Mr. Abdul-Majeed Ayebire is thankful to ICCR for providing financial assistance via reference number LY7930408566174.
Authors are thankful to the worthy reviewers for their valuable comments to improve this manuscript.
Authors declare that they have no conflict of interest.
[1] |
B. Gurbuz, M. Sezer, Laguerre polynomial approach for solving Lane-Emden type functional differential equations, Appl. Math. Comput., 242 (2014), 255–264. http://dx.doi.org/10.1016/j.amc.2014.05.058 doi: 10.1016/j.amc.2014.05.058
![]() |
[2] |
S. Yuzbasi, M. Sezer, An improved Bessel collocation method with a residual error function to solve a class of Lane-Emden differential equations, Math. Comput. Model., 57 (2013), 1298–1311. http://dx.doi.org/10.1016/j.mcm.2012.10.032 doi: 10.1016/j.mcm.2012.10.032
![]() |
[3] |
B. Mehta, R. Aris, A note on a form of the Emden-Fowler equation, J. Math. Anal. Appl., 36 (1971), 611–621. http://dx.doi.org/10.1016/0022-247X(71)90043-6 doi: 10.1016/0022-247X(71)90043-6
![]() |
[4] |
J. Wong, On the generalized Emden-Fowler equation, SIAM Rev., 17 (1975), 339–360. http://dx.doi.org/10.1137/1017036 doi: 10.1137/1017036
![]() |
[5] |
A. Verma, M. Kumar, Numerical solution of third-order Emden-Fowler type equations using artificial neural network technique, Eur. Phys. J. Plus, 135 (2020), 751. http://dx.doi.org/10.1140/epjp/s13360-020-00780-3 doi: 10.1140/epjp/s13360-020-00780-3
![]() |
[6] |
G. File, T. Aga, Numerical solution of quadratic Riccati differential equations, Egyptian Journal of Basic and Applied Sciences, 3 (2016), 392–397. http://dx.doi.org/10.1016/j.ejbas.2016.08.006 doi: 10.1016/j.ejbas.2016.08.006
![]() |
[7] | S. Arora, I. Bala, Numerical study of the coupled Burger and Burger Huxley equations using Bessel collocation scheme, MESA, 14 (2023), 323. |
[8] |
W. Wang, H. Zhang, X. Jiang, X. Yang, A high-order and efficient numerical technique for the nonlocal neutron diffusion equation representing neutron transport in a nuclear reactor, Ann. Nucl. Energy, 195 (2024), 110163. http://dx.doi.org/10.1016/j.anucene.2023.110163 doi: 10.1016/j.anucene.2023.110163
![]() |
[9] |
I. Kaur, S. Arora, I. Bala, An improvised technique of quintic Hermite splines to discretize generalized Burger Huxley type equations, Iranian Journal of Numerical Analysis and Optimization, 13 (2023), 59–79. http://dx.doi.org/10.22067/ijnao.2022.75871.1120 doi: 10.22067/ijnao.2022.75871.1120
![]() |
[10] | S. Arora, R. Jain, V. Kukreja, A robust Hermite spline collocation technique to study generalized Burgers-Huxley equation, generalized Burgers-Fisher equation and Modified Burgers' equation, J. Ocean. Eng. Sci., in press. http://dx.doi.org/10.1016/j.joes.2022.05.016 |
[11] |
M. Hausser, The Hodgkin-Huxley theory of the action potential, Nat. Neurosci., 3 (2000), 1165. http://dx.doi.org/10.1038/81426 doi: 10.1038/81426
![]() |
[12] |
K. Petousakis, A. Apostolopoulou, P. Poirazi, The impact of Hodgkin-Huxley models on dendritic research, J. Physiol., 601 (2023), 3091–3102. http://dx.doi.org/10.1113/JP282756 doi: 10.1113/JP282756
![]() |
[13] |
J. Bisquert, A frequency domain analysis of the excitability and bifurcations of the FitzHugh-Nagumo neuron model, J. Phys. Chem. Lett., 12 (2021), 11005–11013. http://dx.doi.org/10.1021/acs.jpclett.1c03406 doi: 10.1021/acs.jpclett.1c03406
![]() |
[14] |
A. Cevikel, A. Bekir, O. Arqub, M. Abukhaled, Solitary wave solutions of Fitzhugh-Nagumo-type equations with conformable derivatives, Front. Phys., 10 (2022), 1028668. http://dx.doi.org/10.3389/fphy.2022.1028668 doi: 10.3389/fphy.2022.1028668
![]() |
[15] |
J. Villadsen, W. Stewart, Solution of boundary value problem by orthogonal collocation, Chem. Eng. Sci., 20 (1995), 3981–3996. http://dx.doi.org/10.1016/0009-2509(96)81831-8 doi: 10.1016/0009-2509(96)81831-8
![]() |
[16] |
S. Arora, I. Kaur, Applications of quintic Hermite collocation with time discretization to singularly perturbed problems, Appl. Math. Comput., 316 (2018), 409–421. http://dx.doi.org/10.1016/j.amc.2017.08.040 doi: 10.1016/j.amc.2017.08.040
![]() |
[17] |
M. Noor, M. Waseem, Some iterative method for solving a system of nonlinear equations, Comput. Math. Appl., 57 (2009), 101–106. http://dx.doi.org/10.1016/j.camwa.2008.10.067 doi: 10.1016/j.camwa.2008.10.067
![]() |
[18] |
H. Zhang, X. Yang, Q. Tang, D. Xu, A robust error analysis of the OSC method for a multi-term fourth-order sub-diffusion equation, Comput. Math. Appl., 109 (2022), 180–190. http://dx.doi.org/10.1016/j.camwa.2022.01.007 doi: 10.1016/j.camwa.2022.01.007
![]() |
[19] |
S. Arora, S. Dhaliwal, V. Kukreja, Solution of two point boundary value problems using orthogonal collocation on finite elements, Appl. Math. Comput., 171 (2005), 358–370. http://dx.doi.org/10.1016/j.amc.2005.01.049 doi: 10.1016/j.amc.2005.01.049
![]() |
[20] |
S. Arora, I. Kaur, H. Kumar, V. Kukreja, A robust technique of cubic Hermite collocation for solution of two phase non linear model, Journal of King Saud University-Engineering Sciences, 29 (2017), 159–165. http://dx.doi.org/10.1016/j.jksues.2015.06.003 doi: 10.1016/j.jksues.2015.06.003
![]() |
[21] | P. Mishra, K. Sharma, A. Pani, G. Fairweather, Orthogonal spline collocation for singularly perturbed reaction diffusion problems in one dimension, Int. J. Numer. Anal. Mod., 16 (2019), 647–667. |
[22] |
X. Yang, Z. Zhang, On conservative, positivity preserving, nonlinear FV scheme on distorted meshes for the multi-term nonlocal Nagumo-type equations, Appl. Math. Lett., 150 (2024), 108972. http://dx.doi.org/10.1016/j.aml.2023.108972 doi: 10.1016/j.aml.2023.108972
![]() |
[23] |
X. Yang, L. Wu, H. Zhang, A space-time spectral order sinc-collocation method for the fourth-order nonlocal heat model arising in viscoelasticity, Appl. Math. Comput., 457 (2023), 128192. http://dx.doi.org/10.1016/j.amc.2023.128192 doi: 10.1016/j.amc.2023.128192
![]() |
[24] |
B. Bialecki, R. Fernandes, An alternating-direction implicit orthogonal spline collocation scheme for nonlinear parabolic problems on rectangular polygons, SIAM J. Sci. Comput., 28 (2006), 1054–1077. http://dx.doi.org/10.1137/050627885 doi: 10.1137/050627885
![]() |
[25] |
C. Hendricks, M. Ehrhardt, M. Gunther, High-order ADI schemes for diffusion equations with mixed derivatives in the combination technique, Appl. Numer. Math., 101 (2016), 36–52. http://dx.doi.org/10.1016/j.apnum.2015.11.003 doi: 10.1016/j.apnum.2015.11.003
![]() |
[26] |
M. Kadalbajoo, A. Awasthi, A numerical method based on Crank-Nicolson scheme for Burgers' equation, Appl. Math. Comput., 182 (2006), 1430–1442. http://dx.doi.org/10.1016/j.amc.2006.05.030 doi: 10.1016/j.amc.2006.05.030
![]() |
[27] |
Priyanka, S. Arora, F. Mebrek-Oudina, S. Sahani, Super convergence analysis of fully discrete Hermite splines to simulate wave behaviour of Kuramoto-Sivashinsky equation, Wave Motion, 121 (2023), 103187. http://dx.doi.org/10.1016/j.wavemoti.2023.103187 doi: 10.1016/j.wavemoti.2023.103187
![]() |
[28] |
D. Kumar, M. Kadalbajoo, A parameter-uniform numerical method for time-dependent singularly perturbed differential difference equations, Appl. Math. Model., 35 (2011), 2805–2819. http://dx.doi.org/10.1016/j.apm.2010.11.074 doi: 10.1016/j.apm.2010.11.074
![]() |
[29] | S. Rubin, R. Graves, A cubic spline approximation for problems in fluid dynamics, NASA Technical Report, 1975, 19750025272. |
[30] | E. Rainville, Special functions, New York: The Macmillan Company, 1960. |
[31] | I. Sneddon, Special function of mathematical physics and chemistry, 3 Eds, London: Longman Mathematical Texts, 1980. |
[32] |
M. Mazure, On the Hermite interpolation, CR Math. 340 (2005), 177–180. http://dx.doi.org/10.1016/j.crma.2004.11.004 doi: 10.1016/j.crma.2004.11.004
![]() |
[33] |
C. Hall, On error bounds for spline interpolation, J. Approx. Theory, 1 (1968), 209–218. http://dx.doi.org/10.1016/0021-9045(68)90025-7 doi: 10.1016/0021-9045(68)90025-7
![]() |
[34] | P. Prenter, Splines and variational methods, New York: Wiley interscience publication, 1975. |
[35] |
R. Jiwari, R. Gupta, V. Kumar, Polynomial differential quadrature method for numerical solutions of the generalized Fitzhugh-Nagumo equation with time-dependent coefficients, Ain Shams Eng. J., 5 (2014), 1343–1350. http://dx.doi.org/10.1016/j.asej.2014.06.005 doi: 10.1016/j.asej.2014.06.005
![]() |
[36] |
R. Ezzati, K. Shakibi, Using Adomian's decomposition and multiquadric quasi-interpolation methods for solving Newell-Whitehead equation, Procedia Computer Science, 3 (2011), 1043–1048. http://dx.doi.org/10.1016/j.procs.2010.12.171 doi: 10.1016/j.procs.2010.12.171
![]() |
1. | Abdul-Majeed Ayebire, Atul Pasrija, Mukhdeep Singh Manshahia, Shelly Arora, A Novel Hybrid Computational Technique to Study Conformable Burgers’ Equation, 2024, 29, 2297-8747, 114, 10.3390/mca29060114 | |
2. | Jeong-Wook Lee, Gang-Ju Yoo, Jae-Hoon Kim, Seong-Dae Lee, Enhancing navigation safety for small vessels with an electronic chart GPS plotter and autopilot, 2024, 48, 2234-7925, 494, 10.5916/jamet.2024.48.6.494 |
ξ=0 | ξ=1 | ξ=0 | ξ=1 | ||
H1 | 1 | 0 | H′1 | 0 | 0 |
H2 | 0 | 1 | H′2 | 0 | 0 |
H3 | 0 | 0 | H′3 | 1 | 0 |
H4 | 0 | 0 | H′4 | 0 | 1 |
μ=0.25 | μ=−0.5 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 6.53005×10−6 | 6.10156×10−5 | 2.45277×10−5 | 2.42908×10−4 |
0.6 | 1.46464×10−5 | 1.91197×10−4 | 2.39772×10−5 | 2.23646×10−4 |
0.8 | 1.88315×10−5 | 2.46652×10−4 | 2.40393×10−5 | 2.18534×10−4 |
1 | 2.10045×10−5 | 3.66614×10−4 | 2.37522×10−5 | 2.12365×10−4 |
3 | 3.12686×10−5 | 4.85899×10−4 | 3.86837×10−5 | 5.53539×10−4 |
6 | 3.51878×10−5 | 5.44749×10−4 | 5.13096×10−4 | 9.43100×10−3 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 2.420467×10−2 | 3.032328×10−3 | 4.816493×10−5 | 9.623057×10−5 | 1.519802×10−5 |
−6 | 3.805404×10−4 | 1.942463×10−5 | 4.888632×10−6 | 4.283628×10−6 | 4.157908×10−6 |
−4 | 7.190073×10−5 | 1.778285×10−5 | 1.505554×10−5 | 1.498958×10−5 | 1.498550×10−5 |
−2 | 3.583102×10−4 | 1.537630×10−5 | 4.005334×10−5 | 4.066483×10−5 | 4.070273×10−5 |
2 | 3.518331×10−4 | 5.601424×10−5 | 3.817221×10−5 | 3.772792×10−5 | 3.770037×10−5 |
4 | 2.669094×10−5 | 9.245578×10−6 | 1.320194×10−5 | 1.330002×10−5 | 1.330608×10−5 |
6 | 7.101890×10−4 | 3.316729×10−5 | 5.893619×10−6 | 4.357652×10−6 | 3.636630×10−6 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.46995×10−2 | 4.31885×10−3 | 1.18035×10−4 | 1.37984×10−4 | 2.41337×10−5 |
−6 | 5.42797×10−4 | 4.09318×10−5 | 2.07135×10−5 | 1.98629×10−5 | 1.96841×10−5 |
−4 | 6.67514×10−5 | 6.80546×10−5 | 6.95134×10−5 | 6.95535×10−5 | 6.95561×10−5 |
−2 | 2.85835×10−4 | 1.45990×10−4 | 1.76459×10−4 | 1.77204×10−4 | 1.77250×10−4 |
2 | 2.95911×10−4 | 1.38367×10−4 | 1.30731×10−4 | 1.30538×10−4 | 1.30526×10−4 |
4 | 1.37342×10−5 | 3.88842×10−5 | 4.33501×10−5 | 4.34606×10−5 | 4.34674×10−5 |
6 | 6.26818×10−4 | 3.35472×10−5 | 1.32485×10−5 | 1.21421×10−5 | 1.16359×10−5 |
8 | 3.16242×10−2 | 5.11491×10−3 | 2.48175×10−4 | 8.23161×10−5 | 4.64385×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞−norm | L∞-norm |
0.2 | 2.3012×10−6 | 1.0808×10−6 | 4.7416×10−5 | 1.9189×10−5 |
0.5 | 5.5695×10−6 | 3.5274×10−6 | 1.2312×10−4 | 5.6503×10−5 |
1 | 1.1864×10−5 | 1.0948×10−5 | 2.6261×10−4 | 1.9484×10−4 |
1.5 | 1.9400×10−5 | 3.6190×10−5 | 4.2096×10−4 | 5.2431×10−4 |
2 | 2.8162×10−5 | 2.7306×10−5 | 5.9999×10−4 | 5.3048×10−4 |
3 | 4.9735×10−5 | 6.9690×10−5 | 1.0324×10−3 | 9.9714×10−4 |
5 | 1.1395×10−4 | 1.1158×10−4 | 2.3020×10−3 | 1.9205×10−3 |
μ=0.25 | μ=0.05 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 7.99141×10−5 | 5.58252×10−4 | 6.09782×10−6 | 7.95415×10−5 |
0.4 | 7.80824×10−5 | 7.43978×10−4 | 5.60292×10−6 | 8.00388×10−5 |
0.6 | 7.89700×10−5 | 7.70075×10−4 | 1.19943×10−5 | 1.69741×10−4 |
0.8 | 8.83713×10−5 | 7.05602×10−4 | 1.71748×10−5 | 3.17626×10−4 |
1 | 9.50768×10−5 | 1.42997×10−3 | 1.92533×10−5 | 2.93372×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.56960×10−1 | 8.40339×10−2 | 1.84735×10−2 | 1.01369×10−3 | 4.19861×10−4 |
−6 | 2.71322×10−2 | 1.76316×10−2 | 6.73960×10−3 | 4.46380×10−4 | 3.73789×10−4 |
−4 | 2.11950×10−2 | 1.10186×10−3 | 1.70975×10−3 | 2.62175×10−4 | 2.59010×10−4 |
−2 | 2.46704×10−2 | 2.77489×10−3 | 2.77545×10−3 | 2.23101×10−4 | 2.23101×10−4 |
2 | 4.97110×10−3 | 4.97179×10−3 | 1.06710×10−3 | 1.06710×10−3 | 4.38632×10−4 |
4 | 5.34227×10−3 | 5.34505×10−3 | 1.77438×10−4 | 1.77435×10−4 | 1.77435×10−4 |
6 | 4.94204×10−3 | 4.96493×10−3 | 4.16657×10−4 | 4.16660×10−4 | 4.16660×10−4 |
8 | 3.40494×10−3 | 4.05882×10−3 | 6.38229×10−4 | 6.40161×10−4 | 1.80930×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.71171×10−1 | 1.18063×10−2 | 2.34982×10−3 | 1.69686×10−3 | 9.16820×10−5 |
−6 | 1.15766×10−2 | 1.73286×10−3 | 1.04057×10−3 | 4.97002×10−4 | 2.65137×10−5 |
−4 | 1.49387×10−3 | 5.60174×10−4 | 4.02871×10−5 | 4.88264×10−5 | 1.61169×10−5 |
−2 | 1.16173×10−3 | 7.66383×10−5 | 7.72591×10−5 | 7.72621×10−5 | 1.53091×10−5 |
2 | 8.10909×10−5 | 8.08199×10−5 | 8.08200×10−5 | 8.08200×10−5 | 1.84049×10−5 |
4 | 7.97483×10−5 | 8.20047×10−5 | 8.20047×10−5 | 4.23380×10−5 | 4.23380×10−5 |
6 | 3.48375×10−5 | 8.27213×10−5 | 8.27595×10−5 | 4.04893×10−5 | 2.11177×10−5 |
8 | 1.33879×10−3 | 5.81959×10−5 | 8.22777×10−5 | 3.90138×10−5 | 2.23853×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞-norm | L∞-norm |
0.2 | 1.3122×10−6 | 2.1142×10−6 | 1.2350×10−5 | 3.0019×10−5 |
0.5 | 6.0995×10−6 | 4.8016×10−6 | 5.1986×10−5 | 8.6760×10−5 |
1 | 2.1213×10−5 | 2.1057×10−5 | 6.3283×10−4 | 3.7204×10−4 |
1.5 | 3.2340×10−5 | 2.0231×10−5 | 8.5383×10−4 | 3.3954×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−25 | 4.124052×10−2 | 3.071154×10−6 | 3.087868×10−8 | 9.088583×10−9 | 8.092759×10−9 |
−15 | 1.623276×10−3 | 8.594511×10−6 | 6.143898×10−6 | 2.501633×10−8 | 1.573009×10−9 |
15 | 6.569221×10−4 | 3.970126×10−4 | 5.117331×10−5 | 6.449049×10−8 | 6.246470×10−8 |
25 | 5.053591×10−4 | 4.519585×10−6 | 4.490637×10−7 | 5.625901×10−8 | 3.573463×10−8 |
ξ=0 | ξ=1 | ξ=0 | ξ=1 | ||
H1 | 1 | 0 | H′1 | 0 | 0 |
H2 | 0 | 1 | H′2 | 0 | 0 |
H3 | 0 | 0 | H′3 | 1 | 0 |
H4 | 0 | 0 | H′4 | 0 | 1 |
μ=0.25 | μ=−0.5 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 6.53005×10−6 | 6.10156×10−5 | 2.45277×10−5 | 2.42908×10−4 |
0.6 | 1.46464×10−5 | 1.91197×10−4 | 2.39772×10−5 | 2.23646×10−4 |
0.8 | 1.88315×10−5 | 2.46652×10−4 | 2.40393×10−5 | 2.18534×10−4 |
1 | 2.10045×10−5 | 3.66614×10−4 | 2.37522×10−5 | 2.12365×10−4 |
3 | 3.12686×10−5 | 4.85899×10−4 | 3.86837×10−5 | 5.53539×10−4 |
6 | 3.51878×10−5 | 5.44749×10−4 | 5.13096×10−4 | 9.43100×10−3 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 2.420467×10−2 | 3.032328×10−3 | 4.816493×10−5 | 9.623057×10−5 | 1.519802×10−5 |
−6 | 3.805404×10−4 | 1.942463×10−5 | 4.888632×10−6 | 4.283628×10−6 | 4.157908×10−6 |
−4 | 7.190073×10−5 | 1.778285×10−5 | 1.505554×10−5 | 1.498958×10−5 | 1.498550×10−5 |
−2 | 3.583102×10−4 | 1.537630×10−5 | 4.005334×10−5 | 4.066483×10−5 | 4.070273×10−5 |
2 | 3.518331×10−4 | 5.601424×10−5 | 3.817221×10−5 | 3.772792×10−5 | 3.770037×10−5 |
4 | 2.669094×10−5 | 9.245578×10−6 | 1.320194×10−5 | 1.330002×10−5 | 1.330608×10−5 |
6 | 7.101890×10−4 | 3.316729×10−5 | 5.893619×10−6 | 4.357652×10−6 | 3.636630×10−6 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.46995×10−2 | 4.31885×10−3 | 1.18035×10−4 | 1.37984×10−4 | 2.41337×10−5 |
−6 | 5.42797×10−4 | 4.09318×10−5 | 2.07135×10−5 | 1.98629×10−5 | 1.96841×10−5 |
−4 | 6.67514×10−5 | 6.80546×10−5 | 6.95134×10−5 | 6.95535×10−5 | 6.95561×10−5 |
−2 | 2.85835×10−4 | 1.45990×10−4 | 1.76459×10−4 | 1.77204×10−4 | 1.77250×10−4 |
2 | 2.95911×10−4 | 1.38367×10−4 | 1.30731×10−4 | 1.30538×10−4 | 1.30526×10−4 |
4 | 1.37342×10−5 | 3.88842×10−5 | 4.33501×10−5 | 4.34606×10−5 | 4.34674×10−5 |
6 | 6.26818×10−4 | 3.35472×10−5 | 1.32485×10−5 | 1.21421×10−5 | 1.16359×10−5 |
8 | 3.16242×10−2 | 5.11491×10−3 | 2.48175×10−4 | 8.23161×10−5 | 4.64385×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞−norm | L∞-norm |
0.2 | 2.3012×10−6 | 1.0808×10−6 | 4.7416×10−5 | 1.9189×10−5 |
0.5 | 5.5695×10−6 | 3.5274×10−6 | 1.2312×10−4 | 5.6503×10−5 |
1 | 1.1864×10−5 | 1.0948×10−5 | 2.6261×10−4 | 1.9484×10−4 |
1.5 | 1.9400×10−5 | 3.6190×10−5 | 4.2096×10−4 | 5.2431×10−4 |
2 | 2.8162×10−5 | 2.7306×10−5 | 5.9999×10−4 | 5.3048×10−4 |
3 | 4.9735×10−5 | 6.9690×10−5 | 1.0324×10−3 | 9.9714×10−4 |
5 | 1.1395×10−4 | 1.1158×10−4 | 2.3020×10−3 | 1.9205×10−3 |
μ=0.25 | μ=0.05 | |||
τ | L2-norm | L∞-norm | L2-norm | L∞-norm |
0.2 | 7.99141×10−5 | 5.58252×10−4 | 6.09782×10−6 | 7.95415×10−5 |
0.4 | 7.80824×10−5 | 7.43978×10−4 | 5.60292×10−6 | 8.00388×10−5 |
0.6 | 7.89700×10−5 | 7.70075×10−4 | 1.19943×10−5 | 1.69741×10−4 |
0.8 | 8.83713×10−5 | 7.05602×10−4 | 1.71748×10−5 | 3.17626×10−4 |
1 | 9.50768×10−5 | 1.42997×10−3 | 1.92533×10−5 | 2.93372×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.56960×10−1 | 8.40339×10−2 | 1.84735×10−2 | 1.01369×10−3 | 4.19861×10−4 |
−6 | 2.71322×10−2 | 1.76316×10−2 | 6.73960×10−3 | 4.46380×10−4 | 3.73789×10−4 |
−4 | 2.11950×10−2 | 1.10186×10−3 | 1.70975×10−3 | 2.62175×10−4 | 2.59010×10−4 |
−2 | 2.46704×10−2 | 2.77489×10−3 | 2.77545×10−3 | 2.23101×10−4 | 2.23101×10−4 |
2 | 4.97110×10−3 | 4.97179×10−3 | 1.06710×10−3 | 1.06710×10−3 | 4.38632×10−4 |
4 | 5.34227×10−3 | 5.34505×10−3 | 1.77438×10−4 | 1.77435×10−4 | 1.77435×10−4 |
6 | 4.94204×10−3 | 4.96493×10−3 | 4.16657×10−4 | 4.16660×10−4 | 4.16660×10−4 |
8 | 3.40494×10−3 | 4.05882×10−3 | 6.38229×10−4 | 6.40161×10−4 | 1.80930×10−4 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−8 | 3.71171×10−1 | 1.18063×10−2 | 2.34982×10−3 | 1.69686×10−3 | 9.16820×10−5 |
−6 | 1.15766×10−2 | 1.73286×10−3 | 1.04057×10−3 | 4.97002×10−4 | 2.65137×10−5 |
−4 | 1.49387×10−3 | 5.60174×10−4 | 4.02871×10−5 | 4.88264×10−5 | 1.61169×10−5 |
−2 | 1.16173×10−3 | 7.66383×10−5 | 7.72591×10−5 | 7.72621×10−5 | 1.53091×10−5 |
2 | 8.10909×10−5 | 8.08199×10−5 | 8.08200×10−5 | 8.08200×10−5 | 1.84049×10−5 |
4 | 7.97483×10−5 | 8.20047×10−5 | 8.20047×10−5 | 4.23380×10−5 | 4.23380×10−5 |
6 | 3.48375×10−5 | 8.27213×10−5 | 8.27595×10−5 | 4.04893×10−5 | 2.11177×10−5 |
8 | 1.33879×10−3 | 5.81959×10−5 | 8.22777×10−5 | 3.90138×10−5 | 2.23853×10−5 |
[35] | Cubic Hermite splines | [35] | Cubic Hermite splines | |
τ | L2-norm | L2-norm | L∞-norm | L∞-norm |
0.2 | 1.3122×10−6 | 2.1142×10−6 | 1.2350×10−5 | 3.0019×10−5 |
0.5 | 6.0995×10−6 | 4.8016×10−6 | 5.1986×10−5 | 8.6760×10−5 |
1 | 2.1213×10−5 | 2.1057×10−5 | 6.3283×10−4 | 3.7204×10−4 |
1.5 | 3.2340×10−5 | 2.0231×10−5 | 8.5383×10−4 | 3.3954×10−4 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 1.16440×10−11 | 2.92000×10−9 | 9.47831×10−10 |
−15 | 1.37206×10−8 | 3.43787×10−6 | 3.07688×10−9 |
25 | 4.83920×10−10 | 3.38887×10−9 | 8.42926×10−10 |
30 | 1.40960×10−11 | 9.72520×10−11 | 4.76281×10−9 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 8.82590×10−11 | 5.43320×10−9 | 5.3765×10−10 |
−15 | 1.03995×10−7 | 6.39620×10−6 | 2.8219×10−9 |
25 | 1.98915×10−9 | 7.32706×10−9 | 9.7918×10−10 |
30 | 5.79430×10−11 | 2.10500×10−10 | 1.2574×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 2.86520×10−10 | 5.43180×10−9 | 6.26624×10−10 |
−15 | 3.37614×10−7 | 6.39128×10−6 | 2.31068×10−9 |
25 | 4.60456×10−9 | 7.32168×10−9 | 1.04420×10−9 |
30 | 1.34120×10−10 | 2.08880×10−10 | 2.55231×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 6.55230×10−10 | 5.43041×10−9 | 3.48943×10−9 |
−15 | 7.72060×10−7 | 6.38636×10−6 | 2.05489×10−9 |
25 | 8.43340×10−9 | 7.31631×10−9 | 9.50172×10−9 |
30 | 2.45660×10−10 | 2.07260×10−10 | 4.66049×10−8 |
ζ | [36] (ADM) | [36] (MQ) | Cubic Hermite splines |
−25 | 1.45683×10−6 | 5.42902×10−9 | 8.09276×10−9 |
−15 | 1.23638×10−9 | 6.38144×10−6 | 1.57301×10−9 |
25 | 1.35958×10−8 | 7.31090×10−10 | 3.57346×10−8 |
30 | 3.96040×10−10 | 2.05640×10−10 | 7.88424×10−8 |
ζ | n=20 | n=40 | n=100 | n=200 | n=400 |
−25 | 4.124052×10−2 | 3.071154×10−6 | 3.087868×10−8 | 9.088583×10−9 | 8.092759×10−9 |
−15 | 1.623276×10−3 | 8.594511×10−6 | 6.143898×10−6 | 2.501633×10−8 | 1.573009×10−9 |
15 | 6.569221×10−4 | 3.970126×10−4 | 5.117331×10−5 | 6.449049×10−8 | 6.246470×10−8 |
25 | 5.053591×10−4 | 4.519585×10−6 | 4.490637×10−7 | 5.625901×10−8 | 3.573463×10−8 |