Next Article in Journal
An Interactive Estimation of the Distribution Algorithm Integrated with Surrogate-Assisted Fitness
Previous Article in Journal
Two Residual Attention Convolution Models to Recover Underexposed and Overexposed Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Addressing Volterra Partial Integro-Differential Equations through an Innovative Extended Cubic B-Spline Collocation Technique

1
Department of Mathematics and Statistics, College of Science, King Faisal University, P.O. Box 400, Al-Ahsa 31982, Saudi Arabia
2
Department of Mathematics, University of Sargodha, Sargodha 40100, Pakistan
*
Authors to whom correspondence should be addressed.
Symmetry 2023, 15(10), 1851; https://doi.org/10.3390/sym15101851
Submission received: 30 August 2023 / Revised: 22 September 2023 / Accepted: 23 September 2023 / Published: 2 October 2023
(This article belongs to the Section Mathematics)

Abstract

:
This paper introduces a novel collocation scheme based on an extended cubic B-spline for approximating the solution of a second-order partial integro-differential equation. The proposed scheme employs new extended cubic B-splines to discretize the second-order derivatives in the spatial domain, while discretization of spatial derivatives of lower orders is achieved using extended cubic B-spline functions. Temporal derivatives are discretized using the forward difference formula. The stability of the algorithm is assessed using the von Neumann stability method to ensure that error magnification is avoided. Furthermore, convergence analysis of the scheme is provided. Numerical experiments are conducted to validate the effectiveness and efficiency of the proposed scheme. The free parameter is optimized using L 2 and L norms. The computed results are compared with those obtained from various standard numerical schemes found in the literature. Mathematical 12 is used to obtain numerical results.

1. Introduction

Integro-differential equations (IDEs) play a crucial role in modeling dynamical systems across various fields in the applied sciences. They combine both derivatives and integrals of an unknown function, offering practical utility across a wide spectrum of domains, including mathematical modeling, biomedical studies, chemical bonding analysis, aerospace exploration, and industrial mathematics [1,2,3,4,5,6,7,8]. While analytical solutions exist for some IDE-based systems, numerical techniques are often necessary to tackle complex problems effectively. Although IDEs have not received as much attention as more well-known equations like differential equations and integral equations, researchers have been actively exploring numerical approaches for solving IDEs in recent years. Notably, several studies have contributed to presenting numerical solutions for IDEs using different numerical schemes. Desch, Grimmer, and Schappacher [9] discussed considerations for linear integrodifferential equations, emphasizing the fundamental aspects of these equations. Day [10] introduced a note on the numerical solution of integro-differential equations, highlighting early attempts to tackle these equations computationally. Ebadi et al. [11] employed the tau method to present a numerical approach for solving nonlinear Volterra integral-differential equations. Siraj-ul-Islam, Aziz, and Fayyaz et al. [12] proposed a novel approach using Haar wavelets for numerical solutions of integro-differential equations, showcasing the potential of wavelet-based methods. Von et al. [13] introduced wavelet discretizations for parabolic integral-differential equations. Zalik and Vinogradov [14] presented the quasi-elastic method, offering a solution strategy for a specific class of integro-differential equations. Rostami and Maleknejad [15] introduced the projection method as a powerful tool for numerically solving partial integro-differential equations (PIDEs), emphasizing its utility in addressing complex mathematical problems. Long et al. [16] proposed a quasi-wavelet method (QWM) tailored for a specific class of PIDEs, highlighting the potential of wavelet techniques in achieving accurate solutions. Akram et al. [17] developed and analyzed a new extended cubic B-spline (NECuBS) approximation method for solving nonlinear time fractional Klein–Gordon equations, showcasing its robustness and versatility. In this paper, our focus is on approximating a second-order PIDE using a new extended cubic-B-spline-based collocation method. The target PIDE is described by the equation:
s ( z , t ) t = 0 t ( t y ) ψ 2 s ( z , y ) z 2 d y + p ( z , t ) , a z b , t 0 , 0 < ψ < 1 ,
subject to the initial condition,
s ( z , 0 ) = γ ( z )
and the boundary conditions,
s ( a , t ) = 0 , s ( b , t ) = 0 , t 0 .
Here, ψ represents a parameter between 0 and 1, and a, b, γ ( z ) is the initial condition. p ( z , t ) is the right-hand side of PIDE. This term represents the immediate effects and external forces acting on the system. The study of symmetries in PIDEs (1)–(3) is a fundamental topic in the field of mathematical physics and has significant applications in various areas, including mechanics and plasma physics. For example, in plasma physics, understanding the symmetries of IDEs can lead to insights into the behavior of plasma waves and instabilities. In this study, we have employed a novel approximation technique based on the extended cubic B-spline basis. B-spline-based methods have gained widespread popularity in modern times due to their remarkable performance [17,18,19,20]. A spline function of order n is a piecewise polynomial function of degree ( n 1 ) . The term “B-spline” is derived from “basis spline” and was first introduced by I. J. Schoenberg. Additionally, credit is due to Pierre B’ezier at Renault and Birkhoff, Garabedian, and de Boor at General Motors for their work in the early 1960s or late 1950s, which laid the foundation for B-splines and led to numerous research publications (see, for example, [21,22,23]).
The extended cubic B-spline, a fourth-order B-spline with a degree of three, incorporates an additional independent parameter that allows for localized control. Our novel proposed approximation, founded on the basis of the extended cubic B-spline, introduces a unique foundation for handling second-order spatial derivatives [17]. This technique greatly supports accurate and efficient solutions for the specific problem under consideration. Numerous alternative B-spline-centered collocation strategies have been applied to address similar issues within the available literature. Ali et al. [18] introduced a quartic B-spline (QBS) collocation technique, demonstrating its effectiveness in solving PIDEs with weak kernel functions. Gholamian and Nadjafi [19] discussed the application of cubic B-spline (CuBS) collocation methods for solving a class of PIDEs, highlighting their potential in various engineering and scientific contexts. Ali et al. [20] proposed a computational modeling approach based on trigonometric cubic B-spline functions, offering approximate solutions for second-order PIDEs. Akram et al. [24] investigated the application of extended cubic B-splines in solving time fractional telegraph equations, highlighting their role in solving wave-like phenomena. Wasim, Abbas, and Iqbal [25] introduced a novel extended B-spline approximation technique tailored for second-order singular BVPs arising in physiology. The study emphasizes the significance of this method in modeling physiological processes and addressing complex boundary value problems. Stoer and Bulirsch’s book [26] provides a comprehensive overview of numerical techniques, numerical stability, and algorithm design, making it a foundational reference in the field. De Boor [21] explored the convergence properties of odd-degree spline interpolation, contributing to the understanding of spline-based approximation methods. Hall [22] focused on error bounds for spline interpolation, providing insights into the accuracy of spline-based approximations. Prenter’s book [23] offers a comprehensive exploration of spline functions and their applications in approximation theory, numerical analysis, and mathematical modeling, making it a valuable resource for researchers and practitioners in the field.
The structure of this paper is organized as follows: Section 2 covers the discretization of space and time. Section 3 and Section 4 delve into the stability and convergence analysis of the proposed method, respectively. Section 5 presents a comparison of numerical outcomes with various other methods documented in the existing literature. Finally, Section 6 provides a concise summary of the study’s findings.

2. Description of the Method

Adhering to customary notations, consider Δ k = t T as the time increment and h = b a N as the spatial step size, where T and N stand as positive integers. Let us define t f = f Δ k for 0 f T and z n = n h for 0 n N , which establish the time and spatial domain partitions, respectively. The spatial domain a z b is uniformly discretized using knots z j , and the range a z b is segmented into N subintervals [ z j , z j + 1 ] of uniform length h, with j = 0 , 1 , 2 , , N 1 , where a = z 0 < z 1 < < z n 1 < z N = b . In order to derive an approximate solution S ( z , t ) for the exact solution s ( z , t ) of Equation (1), we require
S ( z , t ) = j = 1 N + 1 D j ( t ) B j 4 ( z , σ ) .
In this context, the quantities D j ( t ) represent time-dependent variables that require computation, while B j 4 ( z , σ ) denote the extended cubic B-spline basis functions as detailed in [24].
B j 4 ( z , σ ) = 1 24 h 4 4 h ( 1 σ ) ( z z j ) 3 + 3 σ ( z z j ) 4 , z [ z j , z j + 1 ) ( 4 σ ) h 4 + 12 h 3 ( z z j + 1 ) + 6 h 2 ( 2 + σ ) ( z z j + 1 ) 2 12 h ( z z j + 1 ) 3 3 σ ( z z j + 1 ) 4 , z [ z j + 1 , z j + 2 ) ( 4 σ ) h 4 + 12 h 3 ( z j + 3 z ) + 6 h 2 ( 2 + σ ) ( z j + 3 z ) 2 12 h ( z j + 3 z ) 3 3 σ ( z j + 1 z ) 4 , z [ z j + 2 , z j + 3 ) 4 h ( 1 σ ) ( z j + 4 z ) 3 + 3 σ ( z j + 4 z ) 4 , z [ z j + 3 , z j + 4 ) 0 , otherwise .
Due to the inherent property of localized support within the extended cubic B-spline (ECuBS), exclusively, the basis functions B j 1 4 ( z , σ ) , B j 4 ( z , σ ) , and B j + 1 4 ( z , σ ) remain intact at the grid point z j . As a result, the approximation S f at f t h time level is given as
S ( z , t f ) = S f = j = 1 N + 1 D j f ( t ) B j 4 ( z , σ ) .
The variables denoted as D j f ( t ) are determined through the application of collocation conditions to B j 4 ( z , σ ) , in conjunction with the provided initial and boundary conditions. This process leads to the derivations of the approximations for both S f and its crucial derivatives, yielding the following results:
S f = δ 1 D j 1 f + δ 2 D j f + δ 1 D j + 1 f , ( S z ) f = δ 3 D j 1 f + δ 4 D j f + δ 3 D j + 1 f ,
where δ 1 = 4 σ 24 ,   δ 2 = 8 + σ 12 ,   δ 3 = 1 2 h ,   δ 4 = 0 .
The novel estimation for the second-order spatial derivative is provided in reference [25] as
( S 0 f ) z z = 1 24 h 2 ( 2 ( 14 σ ) D 1 f + 3 ( 3 σ 22 ) D 0 f + 8 ( 7 2 σ ) D 1 f + 14 ( σ 2 ) D 2 f + 6 ( 2 σ ) D 3 f + ( σ 2 ) D 4 f ) , ( S j f ) z z = 1 24 h 2 ( ( 2 σ ) D j 2 f + 4 ( 4 + σ ) D j 1 f + 4 ( 4 + σ ) D j + 1 f + ( 2 σ ) D j + 2 f ) 6 ( 6 + σ ) D j f , j = 1 , 2 , , N 1 . ( S N f ) z z = 1 24 h 2 ( ( σ 2 ) D N 4 f + 6 ( 2 σ ) D N 3 f + 14 ( σ 2 ) D N 2 f + 8 ( 7 2 σ ) D N 1 f + 3 ( 3 σ 22 ) D N f + 2 ( 14 σ ) D N + 1 f ) .
Applying the forward difference scheme for discretizing the time derivative in Equation (1), we acquire
s ( z , t ) t S f + 1 ( z ) S f ( z ) Δ k .
The right-hand side of Equation (1) can be expressed as
0 t ( t y ) ψ s z z ( z , y ) d y + p ( z , t ) = 0 t f + 1 ( t f + 1 y ) ψ s z z ( z , y ) d y + p ( z , t f + 1 ) .
The first component on the right-hand side of Equation (10) is discretized with respect to time as follows:
0 t f + 1 ( t f + 1 y ) ψ s z z ( z , y ) d y = 0 t f + 1 ( y ) ψ s z z ( y , t f + 1 y ) d y , = r = 0 f t r t r + 1 ( y ) ψ s z z ( z , t f + 1 y ) d y , = r = 0 f s z z ( y , t f r + 1 ) t r t r + 1 ( y ) ψ d y , = Δ k 1 ψ 1 ψ r = 0 f s z z ( z , t f r + 1 ) [ ( r + 1 ) 1 ψ ( r ) 1 ψ ] , = Δ k 1 ψ 1 ψ r = 0 f c r s z z ( z , t f r + 1 ) ,
where c r = ( r + 1 ) 1 ψ ( r ) 1 ψ . Note that c 0 = 1 . Thus, Equation (1) becomes
s f + 1 ( z ) s f ( z ) Δ k = Δ k 1 ψ 1 ψ r = 0 f c r s z z f r + 1 ( z ) + p ( z , t f + 1 ) .
Let A = Δ k 2 ψ 1 ψ , so that the last equation becomes
s f + 1 ( z ) A c 0 s z z f + 1 ( z ) = s f ( z ) + A r = 1 f c r s z z f r + 1 ( z ) + Δ k p ( z , t f + 1 ) .
The discretization of the spatial derivative is executed using Equation (9), resulting in the reduction of (11) to
S f + 1 ( z , σ ) A S z z f + 1 ( z , σ ) = S f ( z , σ ) + A r = 1 f c r S z z f r + 1 ( z , σ ) + Δ k p f + 1 ( z ) .
For z = z j , where j = 0 , 1 , 2 , , N , we have
S f + 1 ( z j , σ ) A S z z f + 1 ( z j , σ ) = S f ( z j , σ ) + A r = 1 f c r S z z f r + 1 ( z j , σ ) + Δ k p f + 1 ( z j ) .
Setting
P j = S f ( z j , σ ) + A r = 1 f c r S z z f r + 1 ( z j , σ ) + Δ k p f + 1 ( z j ) ,
we can express this for values of j ranging from 0 to N,
S f + 1 ( z j , σ ) A S z z f + 1 ( z j , σ ) = P j
To derive the new scheme using (8), we need to split (14) for the knots j = 0 , j = 1 , 2 , , N 1 and j = N as
S f + 1 ( z 0 ) A S z z f + 1 ( z 0 ) = P 0 ,
S f + 1 ( z j ) A S z z f + 1 ( z j ) = P j , j = 1 , , N 1 ,
S f + 1 ( z N ) A S z z f + 1 ( z N ) = P N
Now, substituting values of S ( z ) and S z z ( z ) from Equations (7) and (8), respectively, in above three equations, we obtain a system of N + 1 equations in N + 3 unknowns. For the unique solution to this system, the boundary conditions (3) are used to eliminate unknowns D 1 and D N + 1 so that a consistent system is obtained, which can be represented in matrix form as
Z n e w D n e w = P n e w ,
where
Z n e w = α 1 ˜ α 2 ˜ α 4 α 5 α 6 0 0 0 β 1 ˜ β 2 ˜ β 2 β 1 0 0 0 0 0 β 1 β 2 β 3 β 4 β 5 0 0 0 0 β 1 β 2 β 3 β 4 β 5 0 0 0 0 0 β 5 β 4 β 3 β 2 β 1 0 0 0 0 0 β 1 β 2 β 2 ˜ β 1 ˜ 0 0 0 α 6 α 5 α 4 α 2 ˜ α 1 ˜ ,
D n e w = D 0 f + 1 D 1 f + 1 D N 1 f + 1 D N f + 1 T ,
P n e w = P 0 P 1 P 2 P N 1 P N T ,
where α 1 = δ 1 A ( 14 σ 24 h 2 ) , α 2 = δ 2 A ( 3 σ 22 8 h 2 ) , α 3 = δ 1 A ( 7 2 σ 3 h 2 ) , α 4 = A ( 7 σ 2 12 h 2 ) , α 5 = A ( 2 σ 4 h 2 ) , α 6 = A ( σ 2 24 h 2 ) , β 1 = A ( 2 σ 24 h 2 ) , β 2 = δ 1 A ( 4 + σ 6 h 2 ) , β 3 = δ 2 + A ( 6 + σ 4 h 2 ) , α 1 ˜ = ( δ 2 δ 1 ) α 1 + α 2 , α 2 ˜ = α 3 α 1 , β 1 ˜ = ( δ 2 δ 1 ) β 1 + β 2 , β 2 ˜ = β 3 β 1 . To commence the iterations, the initial vector is evaluated by following the same procedure of the previous scheme.

Initial Vector

The initial vector, D 0 = D 0 0 D 1 0 D N 1 0 D N 0 T , can be found by initial condition and boundary values of derivative of initial condition as s ( z j , 0 ) = γ ( z j ) , j = 0 , 1 , 2 , 3 , , N . This system produces an ( N + 1 ) × ( N + 1 ) matrix system of the form
A D 0 = B ,
where
A = δ 3 δ 4 δ 3 0 0 δ 1 δ 2 δ 1 0 0 0 δ 1 δ 2 δ 1 0 0 0 δ 1 δ 2 δ 1 0 0 δ 3 δ 4 δ 3
D 0 = D 0 0 D 1 0 D N 1 0 D N 0 T and B = γ ( z 0 ) γ ( z 1 ) γ ( z N ) γ ( z N + 1 ) T . Upon acquiring the initial vector D 0 , the subsequent recurrence relation governs the time progression of vectors D n , enabling the computation of the approximate solution.

3. Stability Analysis

The von Neumann stability method is employed to demonstrate the stability of the proposed scheme, as discussed in [17], ensuring that errors do not amplify. Given that the presented scheme employs distinct approximation values at various spatial points, including cases for j = 0 , j = N , and j = 1 , , N 1 , it becomes necessary to establish stability for all these specific scenarios. Beginning with the case of j = 0 , we have
S f + 1 ( s 0 ) A S z z f + 1 ( z 0 ) = S f ( z ) + A r = 1 f c r S z z f r + 1 ( z 0 ) .
Using Equations (7) and (8), the above equation reduces to
( δ 1 D 1 f + 1 + δ 2 D 0 f + 1 + δ 1 D 1 f + 1 ) A ( 1 24 h 2 ( 2 ( 14 σ ) D 1 f + 1 + 3 ( 3 σ 22 ) D 0 f + 1 + 8 ( 7 2 σ ) D 1 f + 1 + 14 ( σ 2 ) D 2 f + 1 + 6 ( 2 σ ) D 3 f + 1 + ( σ 2 ) D 4 f + 1 ) ) = ( δ 1 D 1 f + δ 2 D 0 f + δ 1 D 1 f ) + A r = 1 f c r ( 1 24 h 2 ( 2 ( 14 σ ) D 1 f r + 1 + 3 ( 3 σ 22 ) D 0 f r + 1 + 8 ( 7 2 σ ) D 1 f r + 1 + 14 ( σ 2 ) D 2 f r + 1 + 6 ( 2 σ ) D 3 f r + 1 + ( σ 2 ) D 4 f r + 1 ) ) ,
α 1 D 1 f + 1 + α 2 D 0 f + 1 + α 3 D 1 f + 1 + α 4 D 2 f + 1 + α 5 D 3 f + 1 + α 6 D 4 f + 1 = ( δ 1 D 1 f + δ 2 D 0 f + δ 1 D 1 f ) + A r = 1 f c r ( 1 h 2 μ 1 D 1 f r + 1 + μ 2 D 0 f r + 1 + μ 3 D 1 f r + 1 + μ 4 D 2 f r + 1 + μ 5 D 3 f r + 1 + μ 6 D 4 f r + 1 ) ) ,
where μ 1 = ( 14 σ 24 ) , μ 2 = ( 3 σ 22 8 ) , μ 3 = ( 7 2 σ 3 ) , μ 4 = ( 7 σ 2 12 ) , μ 5 = ( 2 σ 4 ) , and μ 6 = ( σ 2 24 ) . By introducing the Fourier mode D j f = ξ f exp ( ι j ϕ h ) into Equation (22), where ϕ represents the mode number, h is the step size, ξ denotes the growth factor, and ι signifies the imaginary unit, we arrive at
E ξ f + 1 e x p ι ϕ h = F ξ f e x p ι ϕ h + A h 2 r = 1 f c r G ξ f r + 1 e x p ι ϕ h ,
where
E = α 1 ξ f + 1 e x p ι ( 1 ) ϕ h + α 2 ξ f + 1 e x p ι ( 0 ) ϕ h + α 3 ξ f + 1 e x p ι ( 1 ) ϕ h + α 4 ξ f + 1 e x p ι ( 2 ) ϕ h + α 5 ξ f + 1 e x p ι ( 3 ) ϕ h + α 6 ξ f + 1 e x p ι ( 4 ) ϕ h , F = δ 1 ξ f e x p ι ( 1 ) ϕ h + δ 2 ξ f e x p ι ( 0 ) ϕ h + δ 1 ξ f e x p ι ( 1 ) ϕ h , G = μ 1 ξ f r + 1 e x p ι ( 1 ) ϕ h + μ 2 ξ f r + 1 e x p ι ( 0 ) ϕ h + μ 3 ξ f r + 1 e x p ι ( 1 ) ϕ h + μ 4 ξ f r + 1 e x p ι ( 2 ) ϕ h + μ 5 ξ f r + 1 e x p ι ( 3 ) ϕ h + μ 6 ξ f r + 1 e x p ι ( 4 ) ϕ h .
Now, rearrange the terms of (23) to obtain
ξ f ( F E + A G E h 2 ) ξ f 1 A G E h 2 r = 2 f c r ξ f r = 0 .
Set b 1 = F E A G E h 2 , b r = A G E h 2 r = 2 f c r in (25) to obtain
ξ f + b 1 ξ f 1 + b 2 ξ f 2 + + b f 1 ξ + b f = 0 .
From Equation (24), it is clear that E and F are positive and G 0 so that the coefficients b 1 , b 2 , b f are positive. For further procedure, we need to follow next theorem
Theorem 1
([26], page 320). For any ξ j , which represents a root of a polynomial p ( ξ ) = a 0 ξ n + + a n with the condition that a 0 0 , it is always true that | ξ j | max 1 , j = 1 n | a j a 0 | .
Demonstrating the stability requires establishing that all the roots ξ j of (26) meet the condition | ξ j | 1 . Since, from Theorem (1), a 0 = 1 and a j > 0 , j = 1 , , f , we have
j = 1 f | a j a 0 | = | ( F + A G h 2 r = 1 f c r ) E | ,
where
r = 1 f c r = r = 2 f [ ( r + 1 ) 1 ψ ( r ) 1 ψ ] = ( f + 1 ) 1 ψ 1 .
Assume R ψ = ( f + 1 ) 1 ψ 1 so that (27) becomes
j = 1 f a j = | F + A G R ψ h 2 E | .
Starting from the definition of G presented in (24), upon setting G = 0 within the same equation, the outcome is that h equates to 0. Consequently, the transformation of (25) occurs, yielding
ξ f F E ξ f 1 = 0 This further leads to the deduction that ξ = 0 or ξ = 1 . Thus, the prerequisite condition for stability ( | ξ j | 1 ) is satisfied. Moving forward, if we consider the scenario where G 0 ( G < 0 ) , the progression from (28) is as follows:
F A G R ψ h 2 < E .
In this manner, the condition of stability ( | ξ j | 1 ) is met. Upon simplification of the preceding inequality, we acquire
cos σ h > ( F + A G + δ 2 h 2 2 δ 1 h 2 )
which implies the unconditional stability of the presented technique for j = 0 . In a similar pattern, the presented method can be proved stable for j = N and j = 1 , , N 1 .

4. Convergence Analysis

This section contains the convergence analysis of the proposed algorithm. For this purpose, the convergence analysis is performed separately for spatial and temporal directions. For spatial convergence, we present the following theorem.
Theorem 2.
When s ^ ( z ) represents the precise solution of Equations (1)–(3), and b ^ ( z ) stands for the B-spline collocation approximation to s ^ ( z ) , the method exhibits secondary order convergence. In such a case, the following holds:
s ^ ( z ) b ^ ( z ) η h 2 ,
where η takes the form of a finite constant, given by η = κ 0 Ł h 2 + R .
Proof. 
Suppose s ^ ( z ) precisely satisfies Equations (1)–(3). In this context, s ^ ( z ) is estimated by the approximation b ^ ( z ) , thereby leading to
b ^ ( z ) = j = 1 N + 1 D ^ j ( t ) B j ( z , σ ) ,
where D ^ = ( D ^ 1 , D ^ 0 , , D ^ N + 1 ) . Additionally, let us assume that b ˜ ( z ) represents the computed cubic B-spline approximation for b ^ ( z ) as follows:
b ˜ ( x ) = j = 1 N + 1 D j ˜ ( t ) B j ( z , σ ) ,
where D ˜ = ( D ˜ 1 , D ˜ 0 , , D ˜ N + 1 ) . In order to estimate the error | s ^ ( z ) b ^ ( z ) | , it becomes essential to ascertain the individual errors | s ^ ( z ) b ˜ ( z ) | and | b ˜ ( z ) b ^ ( z ) | . The computation of b ˜ ( z ) and b ^ ( z ) necessitates the evaluation of values for the vectors D ^ and D ˜ through a pair of linear equations,
Z n e w D ^ = P ^ ,
and
Z n e w D ˜ = P ˜ .
Now, by subtracting (34) from (33), we obtain
Z n e w ( D ˜ D ^ ) = P ˜ P ^ .
By defining matrix Z n e w according to the formulation in Equation (18), the arrangement ensures that Z n e w possesses a strict diagonal dominance, making it nonsingular. Consequently,
( D ˜ D ^ ) = Z n e w 1 ( P ˜ P ^ ) .
Taking the infinity norm of the above equation, we obtain
( D ˜ D ^ ) Z n e w 1 ( P ˜ P ^ ) .
Let us think of τ i ( 0 i N ) as the summation of the i-th row of the matrix Z n e w = [ z i j ] ( N + 1 ) × ( N + 1 ) . Consequently, we arrive at the following relation:
τ 0 = j = 0 N a 0 j = α 1 ˜ + α 2 ˜ + α 4 + α 5 + α 6 = α ,
τ 1 = j = 0 N a 1 j = β 1 ˜ + β 2 ˜ + β 1 + β 2 = β ^ ,
τ i = j = 0 N a i j = β 3 + 2 β 1 + 2 β 2 = β , i = 2 , , N 2 .
τ N 1 = j = 0 N a ( N 1 ) j = β 1 ˜ + β 2 ˜ + β 1 + β 2 = β ^ ,
τ N = j = 0 N a ( N ) j = α 1 ˜ + α 2 ˜ + α 4 + α 5 + α 6 = α .
It is a widely recognized fact in matrix theory that
j = 0 N z i j 1 = 1 , i = 0 , 1 , , N .
Here, z i j 1 are the entries of Z n e w 1 . Thus,
Z n e w 1 = j = 0 N | z i j 1 | 1 τ ,
where τ = min 0 i N τ i = min { α , β , β ^ } . Substituting (38) into (36)
( D ˜ D ^ ) 1 τ ( P ˜ P ^ ) .
Using (13), the upper bound of ( P ˜ P ^ ) is computed for all values of 0 i N and for j = 0 . as
| P i ˜ P ^ i | | s i ˜ s ^ i | + Δ k | g ˜ i f + 1 g ^ i f + 1 | + A h 2 r = 1 f | c r | ( μ 1 | D ˜ 1 f r + 1 D ^ 1 f r + 1 | + μ 2 | D ˜ 0 f r + 1 D ^ 0 f r + 1 | + μ 3 | D ˜ 1 f r + 1 D ^ 1 f r + 1 | + μ 4 | D ˜ 2 f r + 1 D ^ 2 f r + 1 | + μ 5 | D ˜ 3 f r + 1 D ^ 3 f r + 1 | + μ 6 | D ˜ 4 f r + 1 D ^ 4 f r + 1 | ) .
Using the following theorem,
Theorem 3
([21,22]).  If P ( z ) c 4 [ a , b ] , | P 4 ( z ) | Ł , z [ a , b ] and the interval [ a , b ] is partitioned by Δ = { a = z 0 < z 1 < < z N = b } into segments of length h, if b ( z ) stands as the unique spline function that performs interpolation for P ( z ) at the points z 0 , z 1 , , z N , there exists a certain constant κ j with the property that | P ( l ) b ( l ) | κ l · Λ h 4 l for l = 0 , 1 , 2 , 3 .
Then, (40) becomes
| s i ˜ s ^ i | = | b ˜ i ( z ) b ^ i ( z ) | κ 0 Ł h 4 .
Furthermore, c r r = 1 f forms a decreasing sequence of positive values, and c r remains less equal 1 for 1 r n . Consequently, by utilizing (41) and considering g i f + 1 ˜ = g i f ^ + 1 , the Equation (40) transforms to
| P i ˜ P ^ i | κ 0 Ł h 4 + A h 2 r = 1 f d r .
In the above equation, we have taken
( μ 1 | D ˜ 1 f r + 1 D ^ 1 f r + 1 | + μ 2 | D ˜ 0 f r + 1 D ^ 0 f r + 1 | + μ 3 | D ˜ 1 f r + 1 D ^ 1 f r + 1 | + μ 4 | D ˜ 2 f r + 1 D ^ 2 f r + 1 | + μ 5 | D ˜ 3 f r + 1 D ^ 3 f r + 1 | + μ 6 | D ˜ 4 f r + 1 D ^ 4 f r + 1 | ) d r .
Let r = 1 f d r = D f and κ 0 Ł h 4 + A h 2 D f = R f , and then (42) becomes
| P i ˜ P ^ i |   R f .
Using (43) in (39), we obtain
( D ˜ D ^ ) 1 τ R f = R h 2 .
where R h 2 = 1 τ R f = max { α , β , β ^ } . To advance, we require the following theorem.
Theorem 4
([23]).  The set of B-splines B 1 , B 0 , B 1 , , B N 1 , B N , B N + 1 adheres to the following inequality
| i = 1 N + 1 B i ( z ) | 1 , 0 z 1 .
Now, by subtracting (32) from (31), we have
b ˜ ( z ) b ^ ( z ) = j = 1 N + 1 ( D j ˜ D ^ j ) B j 4 ( z ) .
Applying the infinity norm to both sides leads to
b ˜ ( z ) b ^ ( z ) = j = 1 N + 1 ( D j ˜ D ^ j ) B j 4 ( z ) , | j = 1 N + 1 B j 4 ( z ) | ( D j ˜ D ^ j ) ,
that is,
b ˜ ( z ) b ^ ( z ) R h 2 .
From Theorem (3) and Equation (41), we have
s ^ ( z ) b ˜ ( z ) κ 0 Ł h 4 ,
Thus, from (47) and (48)
s ^ ( z ) b ^ ( z ) s ^ ( z ) b ˜ ( z ) + b ˜ ( z ) b ^ ( z ) , κ 0 Ł h 4 + R h 2 , = σ h 2 ,
where σ = κ 0 Ł h 2 + R . Now, to achieve temporal convergence, we can utilize a Taylor expansion on Equation (14), resulting in
( s f ( z ) + Δ k s t f ( z ) + Δ k 2 2 s t t f ( z ) + ) T ( s z z f ( z ) + Δ k s z z t f ( z ) + Δ k s z z t f ( z ) + Δ k 2 2 ! s z z t t f ( z ) + ) = s f ( z ) + A r = 1 f c r s z z f r + 1 ( z ) + Δ ( g f ( z ) + Δ p t f ( z ) ) .
Rearranging terms in the above equation, we obtain
Δ k ( s t f g f ) A ( c 0 s z z f + r = 1 f c r s z z f r + 1 ) + Δ s z z t f + Δ k 2 2 ! s z z t t f + Δ k 2 2 ! ( s z z f + p t f ) = O ( Δ k ) .
Supposing that s ( z , t ) represents the exact solution, while s f ( z , t ) corresponds to the approximated solution of Equations (1)–(3), it follows from (49) that
s ( z , t ) s f ( z , t ) ϱ ( k + h 2 ) .
Here, ϱ denotes a finte constant value. □

5. Numerical Results

The effectiveness and accuracy of the proposed approach are verified in this section through the application to different test problems, employing the L 2 and L error norms as defined by
L 2 = S S f 2 = h j = 0 N | ( S ( z j , t f ) S j f ) 2 |
and
L = S S f = max j | S ( z j , t f ) S j f | .
Example 1.
Consider (1) with
p ( z , t ) = 2 t 1 2 π ( π 5 2 sin ( π z ) 4 t 3 2 π sin ( 2 π z ) ) 2 π 5 2 t 2 sin ( 2 π z ) ,
subject to the boundary conditions,
s ( 0 , t ) = s ( 1 , t ) = 0 , 0 t < 1 ,
and the initial condition
s ( z , 0 ) = sin ( π z ) , z [ 0 , 1 ] .
The exact solution of this problem is s ( z , t ) = sin ( π z ) 4 t 5 2 π sin ( 2 π z ) with ψ = 1 2 . The introduced scheme is applied to the above problem and the numerical results are acquired using Mathematica 12. In Table 1, the computed errors of the present scheme are tabulated and compared to those presented in [16,17,18,19] for various time stages. Figure 1 depicts a close comparison of exact and approximate solutions. In Figure 2, 3D error profiles of the approximate and exact solutions are displayed, which shows the exactness of our methodology. Figure 3 displays the error functions for both 2D and 3D cases. The numerical solution for the scenario where h = 0.05 , t = 1 , and Δ k = 0.01 is provided as
S ( z , 1 ) = 1.73472 × 10 18 1.53732 z + 0.0026847 z 2 + 4.07744 z 3 + 318.921 z 4 , z [ 0 , 1 20 ) 0.0197548 2.79779 z + 28.2192 z 2 244.176 z 3 + 920.358 z 4 , z [ 1 20 , 1 10 ) 0.311013 12.0329 z + 130.518 z 2 684.653 z 3 + 1417.75 z 4 , z [ 1 10 , 3 20 ) 1227.72 + 5589.46 z 9523.68 z 2 + 7211 z 3 2049.49 z 4 , z [ 17 20 , 9 10 ) 990.477 + 4256.43 z 6837.66 z 2 + 4877.41 z 3 1305.73 z 4 , z [ 9 10 , 19 20 ) 446.381 + 1803.22 z 2707.92 z 2 + 1799.53 z 3 448.446 z 4 , z [ 19 20 , 1 ) .
Example 2.
Consider the Equation (1) with exact solution s ( z , t ) = ( t + 1 ) 2 sin π z subject to initial condition s ( z , 0 ) = sin π z and boundary conditions s ( 0 , t ) = s ( 1 , t ) = 0 . Moreover, p ( z , t ) is to be chosen accordingly when ψ = 0.5 .
The proposed technique is implemented for the above problem and the approximate solutions are obtained when h = 0.01 , Δ k = 10 5 for optimized value of σ . Table 2 illustrates the comparison between the calculated errors using the presented method and those from the references [18,19] across various time levels. Figure 4 showcases a substantial disparity between the approximated and actual solutions. Furthermore, Figure 5 provides a remarkably close comparison of the 3D graphs between the approximate and exact solutions. The presentation of 2D and 3D errors can be observed in Figure 6. Additionally, the numerical solution for the case where h = 0.05 , t = 1 , and Δ k = 0.01 is given as
S ( z , 1 ) = 1.11022 × 10 16 + 12.5337 z 9.69194 × 10 7 z 2 61.5212 z 3 + 614.183 z 4 , z [ 0 , 1 20 ) 0.0381009 + 10.096 z + 54.8197 z 2 548.301 z 3 + 1827.33 z 4 , z [ 1 20 , 1 10 ) 0.637203 9.04529 z + 269.596 z 2 1497.86 z 3 + 2995.54 z 4 , z [ 1 10 , 3 20 ) 1758.86 8018.7 z + 13749.2 z 2 10484.3 z 3 + 2995.54 z 4 , z [ 17 20 , 9 10 ) 1343.99 5784.17 z + 9373.93 z 2 6761.04 z 3 + 1827.33 z 4 , z [ 9 10 , 19 20 ) 565.196 2284.7 z + 3500.54 z 2 2395.21 z 3 + 614.183 z 4 , z [ 19 20 , 1 ) .
Example 3.
Consider (1) with exact solution
s ( z , t ) = ( t + 1 ) cos π z , z [ 0 , 1 ] .
The assessment of initial and boundary conditions is based on Equation (52). The selection of the function p ( z , t ) should align accordingly.
The introduced methodology is employed for the above problem to obtain the approximate solution. The computed errors of the presented scheme are compared with those presented in [18,19] for h = 0.05 , Δ t = 10 4 , ψ = 1 3 , and for optimized value of σ for Example 3 in Table 3. In Figure 7, a close contrast of exact and approximate solutions is depicted, while 3D comparison of exact and approximate solutions is plotted in Figure 8. The graphs are in tremendous affirmation. Figure 9 shows the 2D and 3D error graphs. The following numerical solution corresponds to the values h = 0.05 , t = 1 , and Δ k = 0.01 .
S ( z , 1 ) = 2 0.0101225 z 0.48932 z 2 373.465 z 3 + 3735.06 z 4 , z [ 0 , 1 20 ) 2.09164 5.49694 z + 108.786 z 2 1092.76 z 3 + 3642.92 z 4 , z [ 1 20 , 1 10 ) 2.78381 26.0803 z + 310.986 z 2 1730.41 z 3 + 3461.22 z 4 , z [ 1 10 , 3 20 ) 2018.49 + 9249.53 z 15887.1 z 2 + 12114.5 z 3 3461.22 z 4 , z [ 17 20 , 9 10 ) 2655.54 + 11505.5 z 18688 z 2 + 13478.9 z 3 3642.92 z 4 , z [ 9 10 , 19 20 ) 3363.09 + 13818.8 z 21289.5 z 2 + 14566.8 z 3 3735.06 z 4 , z [ 19 20 , 1 ) .
Example 4.
Consider (1) with initial condition
s ( z , 0 ) = 2 sin 2 π z ,
and the boundary conditions
s ( 0 , t ) = 0 , s ( 1 , t ) = 0 .
This equation has the exact solution
s ( z , t ) = 2 ( t 2 + t + 1 ) sin 2 π z ,
p ( z , t ) is to be chosen with ψ = 1 4 . The offered scheme is implemented to obtain the numerical results of Example 4. Errors of the proposed algorithm for h = 0.05 , Δ k = 0.0001 , and Q = 50 , 100 , , 300 are evaluated. Table 4 compares calculated errors of the present scheme with those computed in [18,19]. Figure 10 plots an excellent comparison of approximate and exact solutions for different times. Figure 11 displays the three-dimensional comparison between the approximate and exact solutions. The visualization of both the 2D and 3D error functions is presented in Figure 12. Additionally, the numerical solution corresponding to h = 0.05 , t = 1 , and Δ k = 0.01 is given by
S ( z , 1 ) = 3.46945 × 10 18 0.0136697 z + 0.299777 z 2 + 2326.44 z 3 23265.4 z 4 , z [ 0 , 1 20 ) 0.538757 + 32.0261 z 629.067 z 2 + 6293.66 z 3 20979.8 z 4 , z [ 1 20 , 1 10 ) 3.86934 + 127.617 z 1498.46 z 2 + 8326.4 z 3 16653.7 z 4 , z [ 1 10 , 3 20 ) 9702 + 44504.9 z 76441.4 z 2 + 58288.4 z 3 16653.7 z 4 , z [ 17 20 , 9 10 ) 15283.7 + 66264.4 z 107627 z 2 + 77625.6 z 3 20979.8 z 4 , z [ 9 10 , 19 20 ) 20938.6 + 86081.5 z 132613 z 2 + 90735 z 3 23265.4 z 4 , z [ 19 20 , 1 ) .

6. Conclusions

This study presents an innovative extended cubic B-spline collocation technique as a powerful numerical tool for addressing second-order PIDE. Our proposed approach harnesses the flexibility and accuracy of B-spline functions, providing an effective means of approximating solutions to Volterra PIDEs. A new extended B-spline-based approximation is used to approximate second-order spatial derivatives. The temporal derivatives are discretized by the finite difference formula. The suggested scheme’s stability and convergence analysis are provided to substantiate accuracy and validity. As the understanding of dynamic processes with memory effects continues to grow in importance across various domains, the need for efficient and reliable numerical techniques for solving Volterra PIDEs becomes increasingly crucial. We believe that our proposed method can contribute significantly to advancing research in these areas. The obtained results show that the proposed method offers high accuracy and robustness, making it a valuable addition to the toolkit of researchers and practitioners working in fields where Volterra PIDEs play a significant role. The obtained results are compared with those acquired using TCuBS, CuBS, QBS, and QWM already available in the literature.

Author Contributions

Conceptualization, I.A. and M.Y.; Data curation, S.K.; Formal analysis, S.K.; Funding acquisition, I.A.; Resources, I.A.; Software, M.Y. and S.K.; Writing—original draft, S.K.; Writing—review and editing, M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT No. 4017].

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT No. 4017].

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NECuBSNew extended cubic B-Spline
TCuBSTrigonometric cubic B-Spline
CuBSCubic B-Spline
QBSQuartic B-Spline
QWMQuasi-Wavelet Method

References

  1. Lakshmikantham, V.; Rama, M. Theory of Integro-Differential Equations; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  2. Thiem, H.R. A model for spatio spread of an epidemic. Math. Biol. 1977, 4, 337–351. [Google Scholar] [CrossRef] [PubMed]
  3. Sloan, I.; Thome, V. Time Discretization of an Integro-Differential Equation of Parabolic Type. SIAM J. Numer. Anal. 1986, 23, 1052–1061. [Google Scholar] [CrossRef]
  4. Shakeri, F.; Dehghan, M. A high order finite volume element method for solving elliptic partial integro-differential equations. Appl. Numer. Math. 2013, 65, 105–118. [Google Scholar] [CrossRef]
  5. Cushing, J.M. Integrodifferential Equations and Delay Models in Population Dynamics; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  6. Lin, Y.; Xu, C. Finite difference/spectral approximations for the time-fractional diffusion equation. J. Comput. Phys. 2007, 225, 1533–1552. [Google Scholar] [CrossRef]
  7. Dehghan, M.; Salehi, R. The numerical solution of the non-linear integro-differential equations based on the meshless method. J. Comput. Appl. Math. 2012, 236, 2367–2377. [Google Scholar] [CrossRef]
  8. McLean, W.; Thome, V. Numerical solution of an evolution equation with a positive-type memory term. J. Aust. Math. Soc. Ser. B 1993, 35, 23–70. [Google Scholar] [CrossRef]
  9. Desch, W.; Grimmer, R.; Schappacher, W. Some considerations for Linear integrodifferential equations. J. Math. Anal. 1984, 104, 219–234. [Google Scholar] [CrossRef]
  10. Day, J.T. Note on the numerical solution of Integro-Differential equations. J. Comput. 1967, 9, 394–395. [Google Scholar] [CrossRef]
  11. Ebadi, G.; Rahimiardabili, M.; Shahmorad, S. Numerical solution of the nonlinear Volterra integro-differential equations by the tau method. Appl. Math. Comput. 2007, 188, 1580–1586. [Google Scholar]
  12. Siraj-ul-Islam, A.I.; Fayyaz, M. A new approach for numerical solution of integro-differential equations via Haar wavelets. Int. J. Comput. Math. 2013, 90, 1971–1989. [Google Scholar] [CrossRef]
  13. Von Petersdorff, T.; Schwab, C. Wavelet discretizations of parabolic integrodifferential equations. SIAM J. Numer. Anal. 2003, 41, 159–180. [Google Scholar] [CrossRef]
  14. Zalik, R.; Vinogradov, A. The Quasi-Elastic method of solution for a class of integrodifferential equations. J. Integral Equ. Appl. 1993, 5, 277–295. [Google Scholar] [CrossRef]
  15. Rostami, Y.; Maleknejad, K. Numerical solution of partial Integro-Differential equations by using projection method. Mediterr. J. Math. 2017, 14, 113. [Google Scholar] [CrossRef]
  16. Long, W.T.; Xu, D.; Zeng, X.Y. Quasi wavelet based numerical method for a class of partial integro-differential equation. Appl. Math. Comput. 2012, 218, 11842–11850. [Google Scholar] [CrossRef]
  17. Akram, T.; Abbas, M.; Riaz, M.B.; Ismail, A.I.; Ali, N.H.M. Development and analysis of new approximation of extended cubic b-spline to the nonlinear time fractional klein–gordon equation. Fractals 2020, 28, 2040039. [Google Scholar] [CrossRef]
  18. Ali, A.; Ahmad, S.; Shah, S.I.; Haq, F. A Quartic B-spline collocation technique for the solution of partial integro-differential equations with a weakly kernel. Sci. Int. 2015, 27, 2953–2958. [Google Scholar]
  19. Gholamian, M.; Saberi-Nadjafi, J. Cubic B-splines collocation method for a class of partial integro-differential equation. Alex. Eng. J. 2017, 57, 2157–2165. [Google Scholar] [CrossRef]
  20. Ali, A.; Khan, K.; Haq, F.; Shah, S.I. A Computational Modeling Based on Trigonometric Cubic B-Spline Functions for the Approximate Solution of a Second Order Partial Integro-Differential Equation. In World Conference on Information Systems and Technologies; Rocha, Á., Adeli, H., Reis, L., Costanzo, S., Eds.; Springer International Publishing: Cham, Switzerland, 2019; Volume 930, pp. 844–854. [Google Scholar]
  21. de Boor, C. On the convergence of odd degree spline interpolation. J. Approx. Theory 1968, 1, 452–463. [Google Scholar] [CrossRef]
  22. Hall, C.A. On error bounds for spline interpolation. J. Approx. Theory 1968, 1, 209–218. [Google Scholar] [CrossRef]
  23. Prenter, P.M. Spline and Variational Methods; Wiley: New York, NY, USA, 1975. [Google Scholar]
  24. Akram, T.; Abbas, M.; Ismail, A.I.; Ali, N.H.M.; Baleanu, D. Extended cubic B-splines in the numerical solution of time fractional telegraph equation. Adv. Differ. Equ. 2019, 2019, 365. [Google Scholar] [CrossRef]
  25. Wasim, I.; Abbas, M.; Iqbal, K.I. A new extended B-spline approximation technique for second order singular boundary value problems arising in physiology. J. Math. Comput. Sci. 2019, 19, 258–267. [Google Scholar] [CrossRef]
  26. Stoer, J.; Bulirsch, R. Introduction to Numerical Analysis, 3rd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
Figure 1. Exact and estimated solutions (represented by triangles, stars, and circles) for Example 1 at different time instances, with a fixed value of h = 0.01 .
Figure 1. Exact and estimated solutions (represented by triangles, stars, and circles) for Example 1 at different time instances, with a fixed value of h = 0.01 .
Symmetry 15 01851 g001
Figure 2. Approximated solution (on the left) and precise solution (on the right) for Example 1 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 2. Approximated solution (on the left) and precise solution (on the right) for Example 1 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g002
Figure 3. Error patterns in both two-dimensional and three-dimensional contexts for Example 1 under the parameter settings h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 3. Error patterns in both two-dimensional and three-dimensional contexts for Example 1 under the parameter settings h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g003
Figure 4. Exact and estimated solutions (represented by triangles, stars, and circles) for Example 2 at different time instances with a fixed value of h = 0.01 .
Figure 4. Exact and estimated solutions (represented by triangles, stars, and circles) for Example 2 at different time instances with a fixed value of h = 0.01 .
Symmetry 15 01851 g004
Figure 5. Approximated solution (on the left) and exact solution (on the right) for Example 2 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 5. Approximated solution (on the left) and exact solution (on the right) for Example 2 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g005
Figure 6. Error profiles in both two-dimensional and three-dimensional contexts for Example 2 under the parameter settings h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 6. Error profiles in both two-dimensional and three-dimensional contexts for Example 2 under the parameter settings h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g006
Figure 7. Exact and estimated (depicted by triangles, stars, and circles) solutions for Example 3 across different time points with a step size of h = 0.01 .
Figure 7. Exact and estimated (depicted by triangles, stars, and circles) solutions for Example 3 across different time points with a step size of h = 0.01 .
Symmetry 15 01851 g007
Figure 8. Approximated solution (on the left) and accurate solution (on the right) for Example 3 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 8. Approximated solution (on the left) and accurate solution (on the right) for Example 3 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g008
Figure 9. The estimated solution (on the left) and the exact solution (on the right) for Example 3 with h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 9. The estimated solution (on the left) and the exact solution (on the right) for Example 3 with h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g009
Figure 10. Exact and estimated (indicated by triangles, stars, and circles) solutions for Example 4 at different time instances with a step size of h = 0.01 .
Figure 10. Exact and estimated (indicated by triangles, stars, and circles) solutions for Example 4 at different time instances with a step size of h = 0.01 .
Symmetry 15 01851 g010
Figure 11. The approximated solution (on the left) and the exact solution (on the right) for Example 4 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 11. The approximated solution (on the left) and the exact solution (on the right) for Example 4 under the conditions h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g011
Figure 12. Two-dimensional and three-dimensional representations of error profiles for Example 4 with the parameter values h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Figure 12. Two-dimensional and three-dimensional representations of error profiles for Example 4 with the parameter values h = 1 60 , t = 0.1 , and Δ k = 0.01 .
Symmetry 15 01851 g012
Table 1. Comparative analysis of errors for Example 1 with Δ k = 10 5 and N = 10 .
Table 1. Comparative analysis of errors for Example 1 with Δ k = 10 5 and N = 10 .
QNECuBSTCuBS [20]CuBS [19]QBS [18]QWM [16]
  L 2 L L 2 L L L L
508.40 ×   10 8 1.13 ×   10 7 5.96 ×   10 6 8.93 ×   10 6 1.24 ×   10 6 1.18 ×   10 4 1.58 ×   10 3
1501.48 ×   10 7 2.02 ×   10 7 3.09 ×   10 5 4.42 ×   10 5 6.34 ×   10 6 6.75 ×   10 4 7.89 ×   10 3
2501.93 ×   10 7 2.62 ×   10 7 6.65 ×   10 5 9.45 ×   10 5 1.36 ×   10 5 1.40 ×   10 3 1.61 ×   10 2
3502.29 ×   10 7 3.11 ×   10 7 1.10 ×   10 4 1.56 ×   10 4 2.25 ×   10 5 2.51 ×   10 3 2.53 ×   10 2
4502.59 ×   10 7 3.53 ×   10 7 1.60 ×   10 4 2.28 ×   10 4 3.28 ×   10 5 3.70 ×   10 3 3.46 ×   10 2
Table 2. Comparative analysis of errors for Example 2 with h = 0.01 and Δ k = 10 5 .
Table 2. Comparative analysis of errors for Example 2 with h = 0.01 and Δ k = 10 5 .
QNECuBSTCuBS [20]CuBS [19]
L 2 L L 2 L L 2 L
1505.43 ×   10 9 7.68 × 10 9 2.82 ×   10 8 3.99 ×   10 8 3.93 ×   10 8 5.56 ×   10 8
2006.13 ×   10 9 8.67 ×   10 9 4.56 ×   10 8 6.45 ×   10 8 6.26 ×   10 8 8.86 ×   10 8
2506.43 ×   10 9 9.06 ×   10 9 6.59 ×   10 8 9.31 ×   10 8 8.97 ×   10 8 1.27 ×   10 8
3006.38 ×   10 9 9.03 ×   10 9 8.86 ×   10 8 1.25 ×   10 7 1.20 ×   10 7 1.70 ×   10 7
3506.01 ×   10 9 8.50 ×   10 9 1.14 ×   10 7 1.61 ×   10 7 1.53 ×   10 7 2.17 ×   10 7
4005.34 ×   10 9 7.56 ×   10 9 1.41 ×   10 7 1.99 ×   10 7 1.89 ×   10 7 2.67 ×   10 7
Table 3. Comparative analysis of errors for Example 3 with h = 0.05 and Δ k = 10 4 .
Table 3. Comparative analysis of errors for Example 3 with h = 0.05 and Δ k = 10 4 .
QNECuBSTCuBS [20]CuBS [19]
L 2 L L 2 L L 2 L
504.42 ×   10 8 1.73 ×   10 8 2.73 ×   10 7 4.19 ×   10 7 3.48 ×   10 7 5.34 ×   10 7
1001.39 ×   10 7 1.73 ×   10 7 8.48 ×   10 7 1.26 ×   10 6 1.08 ×   10 6 1.61 ×   10 6
1502.70 ×   10 7 1.73 ×   10 7 1.64 ×   10 6 2.42 ×   10 6 2.09 ×   10 6 3.08 ×   10 6
2004.31 ×   10 7 6.27 ×   10 7 2.61 ×   10 6 3.82 ×   10 6 3.33 ×   10 6 4.87 ×   10 6
2506.16 ×   10 7 9.04 ×   10 7 3.73 ×   10 6 5.51 ×   10 6 4.76 ×   10 6 7.02 ×   10 6
3008.24 ×   10 7 1.20 ×   10 6 4.99 ×   10 6 7.33 ×   10 6 6.36 ×   10 6 9.35 ×   10 6
Table 4. Comparison of errors and assessment of error norms for Example 4 with h = 0.05 and Δ k = 0.0001 .
Table 4. Comparison of errors and assessment of error norms for Example 4 with h = 0.05 and Δ k = 0.0001 .
QNECuBSTCuBS [20]CuBS [19]
L 2 L L 2 L L 2 L
505.58 ×   10 7 8.51 ×   10 7 2.29 ×   10 6 4.26 ×   10 6 2.46 ×   10 6 4.53 ×   10 6
1001.06 ×   10 6 1.50 ×   10 6 7.89 ×   10 6 1.24 ×   10 5 8.44 ×   10 6 1.32 ×   10 5
1501.52 ×   10 6 1.97 ×   10 6 1.61 ×   10 5 2.50 ×   10 5 1.72 ×   10 5 2.67 ×   10 5
2001.97 ×   10 6 2.29 ×   10 6 2.66 ×   10 5 3.94 ×   10 5 2.84 ×   10 5 4.19 ×   10 5
2502.43 ×   10 6 2.53 ×   10 6 3.90 ×   10 5 5.74 ×   10 5 4.17 ×   10 5 6.12 ×   10 5
3002.89 ×   10 6 3.43 ×   10 6 5.33 ×   10 5 7.76 ×   10 5 5.69 ×   10 5 8.29 ×   10 5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ali, I.; Yaseen, M.; Khan, S. Addressing Volterra Partial Integro-Differential Equations through an Innovative Extended Cubic B-Spline Collocation Technique. Symmetry 2023, 15, 1851. https://doi.org/10.3390/sym15101851

AMA Style

Ali I, Yaseen M, Khan S. Addressing Volterra Partial Integro-Differential Equations through an Innovative Extended Cubic B-Spline Collocation Technique. Symmetry. 2023; 15(10):1851. https://doi.org/10.3390/sym15101851

Chicago/Turabian Style

Ali, Ishtiaq, Muhammad Yaseen, and Sana Khan. 2023. "Addressing Volterra Partial Integro-Differential Equations through an Innovative Extended Cubic B-Spline Collocation Technique" Symmetry 15, no. 10: 1851. https://doi.org/10.3390/sym15101851

APA Style

Ali, I., Yaseen, M., & Khan, S. (2023). Addressing Volterra Partial Integro-Differential Equations through an Innovative Extended Cubic B-Spline Collocation Technique. Symmetry, 15(10), 1851. https://doi.org/10.3390/sym15101851

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop