Traditional theories of term structure of interest rate consist of four major classical theories, including Pure Expectation Theory, Liquidity Preference Theory, Preferred Habitat Theory and Market Segmentation Theory. However, they cannot be well interpreted by the traditional static term structure of interest rate methods such as polynomial spline and exponential spline. To address problems on low precision and weak stability of traditional methods in constructing static interest rate term structure curve, in this paper, we introduce the tension interpolation spline based on a fourth-order differential equation with local tension parameters calculated by Generalized Reduced Gradient (GRG) algorithm. Our primary focus is to illustrate its better prediction effect and stability with an empirical study conducted using datum of treasury bonds. Then, we divided the datum into intra-sample datum for estimating tension parameters and out-of-sample datum for evaluating their robustness of predicting stochastics collected from Shanghai Stock Exchange on 2nd February, 2019. According to the principle of total least squares and total least absolute deviations, the result shows that the tension interpolation spline model has better precision and stronger stability in prediction of out-of-sample treasury bonds prices compared with the model established by polynomial spline and exponential spline. In addition, it can better explain the Liquidity Preference Theory, which confirms that it is suitable for constructing the static term structure of interest rates in the securities exchange market.
Citation: Xiangbin Qin, Yuanpeng Zhu. Static term structure of interest rate construction with tension interpolation splines[J]. AIMS Mathematics, 2024, 9(1): 240-256. doi: 10.3934/math.2024014
[1] | Azzh Saad Alshehry, Naila Amir, Naveed Iqbal, Rasool Shah, Kamsing Nonlaopon . On the solution of nonlinear fractional-order shock wave equation via analytical method. AIMS Mathematics, 2022, 7(10): 19325-19343. doi: 10.3934/math.20221061 |
[2] | Rasool Shah, Abd-Allah Hyder, Naveed Iqbal, Thongchai Botmart . Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis. AIMS Mathematics, 2022, 7(11): 19846-19864. doi: 10.3934/math.20221087 |
[3] | Yousef Jawarneh, Humaira Yasmin, M. Mossa Al-Sawalha, Rasool Shah, Asfandyar Khan . Fractional comparative analysis of Camassa-Holm and Degasperis-Procesi equations. AIMS Mathematics, 2023, 8(11): 25845-25862. doi: 10.3934/math.20231318 |
[4] | Ahmad El-Ajou, Rania Saadeh, Moawaih Akhu Dunia, Ahmad Qazza, Zeyad Al-Zhour . A new approach in handling one-dimensional time-fractional Schrödinger equations. AIMS Mathematics, 2024, 9(5): 10536-10560. doi: 10.3934/math.2024515 |
[5] | Manal Alqhtani, Khaled M. Saad, Rasool Shah, Thongchai Botmart, Waleed M. Hamanah . Evaluation of fractional-order equal width equations with the exponential-decay kernel. AIMS Mathematics, 2022, 7(9): 17236-17251. doi: 10.3934/math.2022949 |
[6] | Abdul Hamid Ganie, Fatemah Mofarreh, Adnan Khan . A novel analysis of the time-fractional nonlinear dispersive K(m, n, 1) equations using the homotopy perturbation transform method and Yang transform decomposition method. AIMS Mathematics, 2024, 9(1): 1877-1898. doi: 10.3934/math.2024092 |
[7] | Naveed Iqbal, Saleh Alshammari, Thongchai Botmart . Evaluation of regularized long-wave equation via Caputo and Caputo-Fabrizio fractional derivatives. AIMS Mathematics, 2022, 7(11): 20401-20419. doi: 10.3934/math.20221118 |
[8] | Yousef Jawarneh, Humaira Yasmin, M. Mossa Al-Sawalha, Rasool Shah, Asfandyar Khan . Numerical analysis of fractional heat transfer and porous media equations within Caputo-Fabrizio operator. AIMS Mathematics, 2023, 8(11): 26543-26560. doi: 10.3934/math.20231356 |
[9] | Ahmed M. Zidan, Adnan Khan, Rasool Shah, Mohammed Kbiri Alaoui, Wajaree Weera . Evaluation of time-fractional Fisher's equations with the help of analytical methods. AIMS Mathematics, 2022, 7(10): 18746-18766. doi: 10.3934/math.20221031 |
[10] | Muhammad Imran Liaqat, Sina Etemad, Shahram Rezapour, Choonkil Park . A novel analytical Aboodh residual power series method for solving linear and nonlinear time-fractional partial differential equations with variable coefficients. AIMS Mathematics, 2022, 7(9): 16917-16948. doi: 10.3934/math.2022929 |
Traditional theories of term structure of interest rate consist of four major classical theories, including Pure Expectation Theory, Liquidity Preference Theory, Preferred Habitat Theory and Market Segmentation Theory. However, they cannot be well interpreted by the traditional static term structure of interest rate methods such as polynomial spline and exponential spline. To address problems on low precision and weak stability of traditional methods in constructing static interest rate term structure curve, in this paper, we introduce the tension interpolation spline based on a fourth-order differential equation with local tension parameters calculated by Generalized Reduced Gradient (GRG) algorithm. Our primary focus is to illustrate its better prediction effect and stability with an empirical study conducted using datum of treasury bonds. Then, we divided the datum into intra-sample datum for estimating tension parameters and out-of-sample datum for evaluating their robustness of predicting stochastics collected from Shanghai Stock Exchange on 2nd February, 2019. According to the principle of total least squares and total least absolute deviations, the result shows that the tension interpolation spline model has better precision and stronger stability in prediction of out-of-sample treasury bonds prices compared with the model established by polynomial spline and exponential spline. In addition, it can better explain the Liquidity Preference Theory, which confirms that it is suitable for constructing the static term structure of interest rates in the securities exchange market.
The main objective of the paper is to investigate the qualitative behavior of a differential equation
(b(ϑ)((x(ϑ)+g(ϑ)x(τ(ϑ)))′′′)α)′+η∑i=1ϕi(ϑ)xβ(δi(ϑ))=0, | (1.1) |
where ϑ≥ϑ0>0 and η is a positive integer.
Let us define the corresponding function υ for the solution x as follows:
υ(ϑ)=x(ϑ)+g(ϑ)x(τ(ϑ)). | (1.2) |
Furthermore, we make the following supposition:
(M1) α and β are quotients of odd positive integers;
(M2) b∈C([ϑ0,∞),(0,∞)), b′(ϑ)≥0 and
∫ϑϑ01b1/α(ζ)dζ→∞ as ϑ→∞; | (1.3) |
(M3) τ,δi∈C([ϑ0,∞),R), δi(ϑ)≤ϑ, τ(ϑ)≤ϑ, τ′(ϑ)>0, and limϑ→∞τ(ϑ)=limϑ→∞δi(ϑ)=∞, i=1,2,...η;
(M4) ϕi,g∈C([ϑ0,∞),R), g(ϑ)>0, ϕi≥0 and ϕi is not identically zero for large ϑ, i=1,2,...η;
(M5) there exists a constant ε∈(0,1) such that
limϑ→∞(ϑτ(ϑ))3/ε1g(ϑ)=0. | (1.4) |
By a solution of (1.1), we mean a nontrivial function x∈C([ϑx,∞),R), ϑx≥ϑ0, which has the properties υ(ϑ)∈C3([ϑx,∞),R), b(ϑ)(υ′′′(ϑ))α∈C1([ϑx,∞),R) and x satisfies (1.1) on [ϑx,∞). We focus in our study on the solutions that satisfy sup{|x(ϑ)|:ϑ≥ϑa}>0, for every ϑa≥ϑx, and we assume that (1.1) possesses such solutions. Such a solution of (1.1) is called oscillatory if it is neither eventually positive nor eventually negative; otherwise, it is called nonoscillatory. An equation is called oscillatory if all of its solutions are oscillatory.
We frequently see repetitive movements used to illustrate different mechanical actions that occur in nature. In other words, oscillations abound in our universe. Examples of oscillatory motions include a ship moving up and down on the waves and a pendulum swing. Finding new necessary conditions for the oscillation or nonoscillation of the solutions of neutral functional differential equations, which is a component of so-called dynamical systems, has become increasingly important in recent years. As far as physicists and engineers are concerned, understanding and managing oscillations in various systems is their primary objective. Oscillation is a phenomenon that is observed in a variety of fields, including biology and economics, in addition to physics and mechanics.
One of the important differential equation branching problems is the oscillatory behavior of ordinary differential equations. The oscillatory problems to the wings of the plane can be modeled by the oscillatory problems of ordinary differential equations. Delay differential equations with deviating arguments introduce an additional layer of complexity by incorporating time delays into the modeling process. Unlike ordinary differential equations, these equations account for the influence of both current and past values of variables. In fact, differential equations with deviating arguments are widely used in physics, engineering, biology, economics, and more, making them an indispensable tool for understanding and predicting the behavior of complex phenomena, see [8,9,13,15,21,23,26,27]).
Bills and Schoenberg [12] examined certain oscillatory outcomes for a self-adjoint system of second-order equations. The oscillatory behavior of solutions to various classes of differential equations with a linear neutral term has been extensively studied in recent decades. However, there are few results about the oscillation of differential equations with nonlinear neutral terms; see, for example, Agarwal et al. and Grace and Graef [3,14,16].
Many researchers have studied the oscillatory properties of even-order differential equations on a larger scale than their odd-order counterparts. Different techniques and methods have been used to study the oscillation of different types of even-order differential equations. Illustrative and additional information can be found in references [4,11,17,18,19,20,25]. Specifically, we provide some detail.
Bazighifan and Ahmad [10] investigated the qualitative behavior of an even-order advanced differential equation
(b(ϑ)(x(n−1)(ϑ))α)′+η∑i=1ϕi(ϑ)f(x(δi(ϑ)))=0, | (1.5) |
where ϕi(ϑ)≥0, δi(ϑ)≥ϑ, and f(x)/xβ≥k>0 for x≠0. They established sufficient conditions for oscillation of (1.5) by utilizing the theory of comparison with first-order and second-order delay equations, as well as the Riccati substitution technique.
The oscillatory behavior of the differential equations
(b(ϑ)((x(ϑ)+g(ϑ)x(τ(ϑ)))′)α)′+η∑i=1ϕi(ϑ)xβ(δi(ϑ))=0 | (1.6) |
was taken into consideration by Abdelnaser et al. [1]. By introducing a new set of criteria, the researchers were able to prove that all solutions to Eq (1.6) oscillate.
Muhib et al. [22] considered a class of neutral differential equations
(b(ϑ)((x(ϑ)+ρ(ϑ)xℏ(σ1(ϑ))+g(ϑ)xγ(τ(ϑ)))(n−1))α)′+f(ϑ,x(δ(ϑ)))=0. | (1.7) |
They used Riccati transformations to present new conditions for the oscillation of (1.7), where ρ(ϑ)≥0, g(ϑ)≥0, γ and ℏ are ratios of odd positive integers with γ≥1, 0<ℏ<1, and f∈C([ϑ0,∞)×R,R) and there exists ϕ∈C([ϑ0,∞),(0,∞)) such that |f(ϑ,x)|≥ϕ(ϑ)|x|β. Below, we present one of the results in [22].
Theorem 1.1. Assume that
limϑ→∞g(ϑ)(ϑn−2∫ϑϑ01b1/α(ζ)dζ)γ−1=limϑ→∞ρ(ϑ)=0 | (1.8) |
holds. If
liminfϑ→∞1φ(ϑ)∫∞ϑαλNδn−2(ζ)δ′(ζ)b−1/α(ζ)φ(α+1)/α(ζ)dζ>α(α+1)(α+1)/α, |
for all λ∈(0,1) and N>0, then (1.7) is oscillatory, where
φ(ϑ)=ϵβ∫∞ϑϕ(u)Ω(u)du |
and
Ω(ϑ)={kβ−α1if β≥α,kβ−α2(ϑn−2)β−α(∫ϑϑ11b1/ℓ(ζ)dζ)β−αif β<α, |
for some ϵ∈(0,1) and k1,k2>0.
Agarwal et al. [2] studied the oscillation of a neutral differential equation
(x(ϑ)+g(ϑ)x(τ(ϑ)))(n)+ϕ(ϑ)x(δ(ϑ))=0, | (1.9) |
where n is even and g(ϑ)≥0. They established some sufficient conditions for oscillation of (1.9) using the Riccati transformation technique. Below, we present one of the results in [2].
Theorem 1.2. Let n≥4 be even (M3), (M4), and
δ(ϑ)≤τ(ϑ), 1−Υn−1(ϑ)g(τ−1(τ−1(ϑ)))≥0 |
hold. If the equation
((n−2)!λ0ϑn−2x′(ϑ))′+ϕ(ϑ)g∗(δ(ϑ))(τ−1(δ(ϑ)))n−1ϑn−1x(ϑ)=0 |
is oscillatory for some constant λ0∈(0,1), and the equation
x′′(ϑ)+∫∞ϑ(s−ϑ)n−3ϕ(s)g∗(δ(s))τ−1(δ(s))sds(n−3)!x(ϑ)=0 |
is oscillatory, then (1.9) is oscillatory, where
g∗(ϑ)=1g(τ−1(ϑ))(1−Υ(ϑ)g(τ−1(τ−1(ϑ)))), |
g∗(ϑ)=1g(τ−1(ϑ))(1−Υn−1(ϑ)g(τ−1(τ−1(ϑ)))) |
and
Υ(ϑ)=τ−1(τ−1(ϑ))τ−1(ϑ). |
The objective of this paper is to provide new results of oscillation (1.1) in canonical form, which would improve and extend some previous literature. In addition, an example is given that shows the applicability of the results we obtained.
The following notation will be used in the remaining sections of this work:
δ(ϑ)=min{δi(ϑ): i=1,2,...,η}. |
For the proof of our main results, we need to give the following lemmas:
Lemma 2.1. [5] Let f∈Cn([ϑ0,∞),(0,∞)), the derivative f(n)(ϑ) is of fixed sign and not identically zero on a subray of [ϑ0,∞), and there exists a ϑx≥ϑ0 such that f(n−1)(ϑ)f(n)(ϑ)≤0 for all ϑ≥ϑ1. If limϑ→∞f(ϑ)≠0, then for every λ∈(0,1) there exists a ϑλ≥ϑ1 such that
|f(ϑ)|≥λ(n−1)!ϑn−1|f(n−1)(ϑ)|, | (2.1) |
for all ϑ≥ϑλ.
Lemma 2.2. [7, Lemma 1] Let the function f satisfy f(n)(ϑ)>0, n=1,2,...,κ, and f(κ+1)(ϑ)≤0 eventually. Then, for every ε∈(0,1),
f(ϑ)f′(ϑ)≥εϑκ, | (2.2) |
eventually.
Lemma 2.3. [6] Let x be a positive solution of (1.1), and (1.3) holds. Then, (b(ϑ)(υ′′′(ϑ))α)′<0, we also find that there exist two potential cases eventually, which are as follows:
Case (1): υ(ϑ)>0,υ′(ϑ)>0,υ′′(ϑ)>0,υ′′′(ϑ)>0,υ(4)(ϑ)≤0,Case (2): υ(ϑ)>0,υ′(ϑ)>0,υ′′(ϑ)<0, υ′′′(ϑ)>0,υ(4)(ϑ)≤0. |
Lemma 2.4. [24] If y is a positive and strictly decreasing solution of the integral inequality
y(ϑ)≥∫∞ϑ(ν−ϑ)n−1(n−1)!f(ν,y(g1(ν)),y(g2(ν)),...,y(gm(ν)))dν, |
then there exists a positive solution x(ϑ) of the differential equation
(−1)nx(n)(ϑ)=f(ϑ,x(g1(ϑ)),x(g2(ϑ)),...,x(gm(ϑ))),,ϑ≥ϑ0 |
being such that x(ϑ)≤y(ϑ) for all large ϑ and satisfying limϑ→∞x(i)(ϑ)=0 monotonically (i=1,2,...,n−1), where f is a continuous function defined on [ϑ0,∞)×[0,∞)m and gj(ϑ) are continuous functions on the interval [ϑ0,∞) such that
limϑ→∞gj(ϑ)=∞(j=1,2,...,m). |
The function f=f(ϑ,u1,u2,...,um) is assumed to be increasing in each of u1,u2,...,um. Moreover, it is supposed that f is positive on [ϑ0,∞)×[0,∞)m and that the delays ϑ−gj(ϑ) are positive for ϑ≥ϑ0, i.e., gj(ϑ)<ϑ for every ϑ≥ϑ0 and (j=1,2,...,m).
We now present the main results of this paper.
Theorem 3.1. Assume that β≥1 and that there exists a positive function μ∈C1([ϑ0,∞),R) such that
μ(ϑ)<τ(ϑ), μ(ϑ)≤δi(ϑ), μ′(ϑ)>0 and limϑ→∞μ(ϑ)=∞. | (3.1) |
If
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1cβ−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0 | (3.2) |
and
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2cβ−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0 | (3.3) |
have no positive solutions, then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)), c>0 and ϵ1,ϵ2∈(0,1).
Proof. Assume that Eq (1.1) has a nonoscillatory solution x(ϑ), say x(ϑ)>0, x(δ(ϑ))>0, and x(τ(ϑ))>0 for ϑ≥ϑ1≥ϑ0. From (1.2), we find
x(ϑ)=υ(τ−1(ϑ))−x(τ−1(ϑ))g(τ−1(ϑ)) |
and so
x(ϑ)≥υ(τ−1(ϑ))g(τ−1(ϑ))−υ(τ−1(τ−1(ϑ)))g(τ−1(ϑ))g(τ−1(τ−1(ϑ))). | (3.4) |
We first consider what Case (1) holds. Since κ=3, in view of (2.2), for every ε∈(0,1), we get
υ(ϑ)υ′(ϑ)≥εϑκ≥εϑ3, | (3.5) |
now,
(υ(ϑ)ϑ3/ε)′=εϑ3/ευ′(ϑ)−3υ(ϑ)ϑ(3/ε)−1εϑ2(3/ε)=εϑυ′(ϑ)−3υ(ϑ)εϑ(3/ε)+1, | (3.6) |
using (3.5), we find
(υ(ϑ)ϑ3/ε)′=εϑυ′(ϑ)−3υ(ϑ)εϑ(3/ε)+1≤0. | (3.7) |
Since τ(ϑ)≤ϑ and τ′(ϑ)>0, (τ−1(ϑ))′>0 and furthermore ϑ≤τ−1(ϑ). Thus,
τ−1(ϑ)≤τ−1(τ−1(ϑ)). | (3.8) |
By using (3.7) and (3.8), it follows that
υ(τ−1(ϑ))(τ−1(ϑ))3/ε≥υ(τ−1(τ−1(ϑ)))(τ−1(τ−1(ϑ)))3/ε |
and so
(τ−1(τ−1(ϑ)))3/ευ(τ−1(ϑ))≥(τ−1(ϑ))3/ευ(τ−1(τ−1(ϑ))). | (3.9) |
From (3.4) and (3.9), we find
x(ϑ)≥υ(τ−1(ϑ))g(τ−1(ϑ))−(τ−1(τ−1(ϑ))τ−1(ϑ))3/ευ(τ−1(ϑ))g(τ−1(ϑ))g(τ−1(τ−1(ϑ)))≥υ(τ−1(ϑ))g(τ−1(ϑ))(1−(τ−1(τ−1(ϑ))τ−1(ϑ))3/ε1g(τ−1(τ−1(ϑ)))). | (3.10) |
From (M5), there exists a ϵ1∈(0,1) such that
(τ−1(τ−1(ϑ))τ−1(ϑ))3/ε1g(τ−1(τ−1(ϑ)))≤1−ϵ1. |
Using the above inequality in (3.10) gives
x(ϑ)≥υ(τ−1(ϑ))g(τ−1(ϑ))ϵ1. | (3.11) |
From (1.1) and (3.11), we have
(b(ϑ)(υ′′′(ϑ))α)′+η∑i=1ϕi(ϑ)υβ(τ−1(δi(ϑ)))gβ(τ−1(δi(ϑ)))ϵβ1≤0 |
and so
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1υβ(τ−1(δ(ϑ)))η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ)))≤0. | (3.12) |
In view of the fact that μ(ϑ)≤δ(ϑ) and υ′(ϑ)>0, inequality (3.12) becomes
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1υβ(τ−1(μ(ϑ)))η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ)))≤0. | (3.13) |
Since υ(ϑ)>0 and υ′(ϑ)>0, there exists a constant c>0 such that
υ(ϑ)≥c. | (3.14) |
From (3.13), (3.14), and β≥1, we find the following differential inequality:
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1cβ−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0, | (3.15) |
has a positive solution υ. This implies that (3.2) also has a positive solution, which contradicts our assumption.
Now, we consider what Case (2) holds. Since κ=1, in view of (2.2), for every ε∈(0,1), we get
υ(ϑ)υ′(ϑ)≥εϑ1, | (3.16) |
from which we see that
(υ(ϑ)ϑ1/ε)′=εϑ1/ευ′(ϑ)−υ(ϑ)ϑ(1/ε)−1εϑ2/ε=εϑυ′(ϑ)−υ(ϑ)εϑ1+(1/ε)≤0. | (3.17) |
By (3.8) and (3.17),
(τ−1(ϑ))1/ευ(τ−1(τ−1(ϑ)))≤(τ−1(τ−1(ϑ)))1/ευ(τ−1(ϑ)). | (3.18) |
Combining (3.4) and (3.18), we obtain
x(ϑ)≥υ(τ−1(ϑ))g(τ−1(ϑ))(1−(τ−1(τ−1(ϑ))τ−1(ϑ))1/ε1g(τ−1(τ−1(ϑ)))). | (3.19) |
From (M5), there exists a ϵ2∈(0,1) such that
(τ−1(τ−1(ϑ))τ−1(ϑ))1/ε1g(τ−1(τ−1(ϑ)))≤1−ϵ2, |
and using this in (3.19) implies
x(ϑ)≥ϵ2υ(τ−1(ϑ))g(τ−1(ϑ)). | (3.20) |
Using (3.20) in (1.1) yields
(b(ϑ)(υ′′′(ϑ))α)′+η∑i=1ϕi(ϑ)ϵβ2υβ(τ−1(δi(ϑ)))gβ(τ−1(δi(ϑ)))≤0 |
and so
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2υβ(τ−1(δ(ϑ)))η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ)))≤0. | (3.21) |
Since μ(ϑ)≤δ(ϑ) and υ′(ϑ)>0, then (3.21) becomes
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2υβ(τ−1(μ(ϑ)))η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ)))≤0. | (3.22) |
In view of (3.14) and β≥1, we find the following differential inequality:
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2cβ−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0, | (3.23) |
has a positive solution υ. This implies that (3.3) also has a positive solution, which contradicts our assumption. The proof is now complete.
Theorem 3.2. Assume that β<1 and there exists a positive function μ∈C1([ϑ0,∞),R) such that (3.1) holds. If
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1dβ−11(q3/ε(ϑ))β−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0 | (3.24) |
and
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2dβ−12(q1/ε(ϑ))β−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0, | (3.25) |
have no positive solutions, then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)), d1,d2>0 and ϵ1,ϵ2∈(0,1).
Proof. Assume that Eq (1.1) has a nonoscillatory solution x(ϑ), say x(ϑ)>0, x(δ(ϑ))>0, and x(τ(ϑ))>0 for ϑ≥ϑ1≥ϑ0.
We first consider what Case (1) holds. By performing the same steps as in the proof of Theorem 3.1, we arrive at (3.7) and (3.13). By (3.7), there exists a constant d1>0 such that
υ(ϑ)ϑ3/ε≤d1 |
and so
υ(ϑ)≤d1ϑ3/ε. | (3.26) |
Using (3.26) in (3.13) and applying the fact that β<1 yields
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1dβ−11(q3/ε(ϑ))β−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0. | (3.27) |
That is, (3.24) has a positive solution, a contradiction.
Now, we consider what Case (2) holds. By performing the same steps as in the proof of Theorem 3.1, we arrive at (3.17) and (3.22). By (3.17), there exists a constant d2>0 such that
υ(ϑ)ϑ1/ε≤d2 |
and so
υ(ϑ)≤d2ϑ1/ε. | (3.28) |
Using (3.28) in (3.22) and applying the fact that β<1 yields
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ2dβ−12(q1/ε(ϑ))β−1(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ(q(ϑ))≤0. | (3.29) |
That is, (3.25) has a positive solution, a contradiction. The proof is now complete.
Theorem 3.3. Assume that β≥1 and there exists a positive function μ∈C1([ϑ0,∞),R) such that (3.1) holds. If
y′(ϑ)+λ3!ϵβ1cβ−1q3(ϑ)b1/α(q(ϑ))(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))y1/α(q(ϑ))=0 | (3.30) |
and
ω′(ϑ)+ϵβ/α2c(β−1)/αε−1/αq−1/α(ϑ)(∫∞ϑ(∫∞u(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αb1/α(u)du)ω1/α(q(ϑ))=0 | (3.31) |
are oscillatory, for some constants λ,ε∈(0,1), then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)), c>0 and ϵ1,ϵ2∈(0,1).
Proof. Assume that Eq (1.1) has a nonoscillatory solution x(ϑ), say x(ϑ)>0, x(δ(ϑ))>0, and x(τ(ϑ))>0 for ϑ≥ϑ1≥ϑ0.
We first consider what Case (1) holds. By performing the same steps as in the proof of Theorem 3.1, we arrive at (3.15). Since υ(ϑ)>0 and υ′(ϑ)>0, we have limϑ→∞ υ(ϑ)≠0. Thus, by Lemma 2.1, we obtain
υ(ϑ)≥λ3!ϑ3υ′′′(ϑ), | (3.32) |
from which we see that
υ(q(ϑ))≥λ3!q3(ϑ)υ′′′(q(ϑ)). | (3.33) |
Using (3.33) in (3.15) yields
(b(ϑ)(υ′′′(ϑ))α)′+λ3!ϵβ1cβ−1q3(ϑ)(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ′′′(q(ϑ))≤0, |
With y(ϑ)=b(ϑ)(υ′′′(ϑ))α, we find y(ϑ) is a positive solution of the differential inequality
y′(ϑ)+λ3!ϵβ1cβ−1q3(ϑ)b1/α(q(ϑ))(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))y1/α(q(ϑ))≤0. | (3.34) |
It follows from Lemma 2.4 that the differential equation (3.30) also has a positive solution for all λ∈(0,1), but this contradicts our assumption on (3.30).
Now, we consider what Case (2) holds. By performing the same steps as in the proof of Theorem 3.1, we arrive at (3.16) and (3.23). Integrating (3.23) from ϑ to ∞ gives
(υ′′′(ϑ))α≥ϵβ2cβ−1υ(q(ϑ))b(ϑ)∫∞ϑ(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ |
and so
υ′′′(ϑ)≥ϵβ/α2c(β−1)/αυ1/α(q(ϑ))b1/α(ϑ)(∫∞ϑ(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/α. | (3.35) |
Integrating (3.35) from ϑ to ∞, we have
−υ′′(ϑ)≥ϵβ/α2c(β−1)/α(∫∞ϑ(∫∞u(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αb1/α(u)du)υ1/α(q(ϑ)) |
and so
υ′′(ϑ)+ϵβ/α2c(β−1)/α(∫∞ϑ(∫∞u(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αb1/α(u)du)υ1/α(q(ϑ))≤0. | (3.36) |
Using (3.16) in (3.36) yields
υ′′(ϑ)+ϵβ/α2c(β−1)/αε−1/αq−1/α(ϑ)(∫∞ϑ(∫∞u(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αb1/α(u)du)(υ′(q(ϑ)))1/α≤0. | (3.37) |
With ω(ϑ)=υ′(ϑ), we see that ω(ϑ) is a positive solution of the differential inequality
ω′(ϑ)+ϵβ/α2c(β−1)/αε−1/αq−1/α(ϑ)(∫∞ϑ(∫∞u(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αb1/α(u)du)ω1/α(q(ϑ))≤0, | (3.38) |
for every ε∈(0,1). We finalize the proof using the same method as outlined in Case (1). The proof is now complete.
Corollary 3.1. Let α=1 and β≥1 hold. Assume further that there exists a positive function μ∈C1([ϑ0,∞),R) such that (3.1) holds. If
limϑ→∞∫ϑq(ϑ)q3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ=∞ | (3.39) |
and
limϑ→∞∫ϑq(ϑ)q(ζ)(∫∞ζ(∫∞u(η∑i=1ϕi(ϱ)gβ(τ−1(δi(ϱ))))dϱ)b(u)du)dζ=∞, | (3.40) |
then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)).
Proof. We first consider what Case (1) holds. By performing the same steps as in the proof of Theorem 3.3, we arrive at (3.34). Integrating (3.34) from q(ϑ) to ϑ and then using α=1 and the fact that y is a decreasing function, we see that
∫ϑq(ϑ)y′(ζ)dζ≤−λ3!ϵβ1cβ−1y(q(ζ))∫ϑq(ϑ)q3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ |
and so
−y(q(ϑ))≤−λ3!ϵβ1cβ−1y(q(ϑ))∫ϑq(ϑ)q3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ, |
this can be expressed as follows:
3!λϵβ1cβ−1≥∫ϑq(ϑ)q3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ, |
so this contradicts (3.39).
Now, we consider what Case (2) holds. By performing the same steps as in the proof of Theorem 3.3, we arrive at (3.38). Integrating (3.38) from q(ϑ) to ϑ and then using α=1 and the fact that ω is a decreasing function, we see that
∫ϑq(ϑ)ω′(ζ)dζ≤−εϵβ2c(β−1)ω(q(ϑ))∫ϑq(ϑ)q(ζ)(∫∞ζ(∫∞u(η∑i=1ϕi(ϱ)gβ(τ−1(δi(ϱ))))dϱ)b(u)du)dζ |
and so
−ω(q(ϑ))≤−εϵβ2c(β−1)ω(q(ϑ))∫ϑq(ϑ)q(ζ)(∫∞ζ(∫∞u(η∑i=1ϕi(ϱ)gβ(τ−1(δi(ϱ))))dϱ)b(u)du)dζ, |
this can be expressed as follows:
1εϵβ2c(β−1)≥∫ϑq(ϑ)q(ζ)(∫∞ζ(∫∞u(η∑i=1ϕi(ϱ)gβ(τ−1(δi(ϱ))))dϱ)b(u)du)dζ, |
which contradicts (3.40). The proof is now complete.
Theorem 3.4. Assume that β<1 and there exists a positive function μ∈C1([ϑ0,∞),R) such that (3.1) holds. If
y′(ϑ)+ϵβ1dβ−11λ3!q(3(β−1)/ε)+3(ϑ)b1/α(q(ϑ))(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))y1/α(q(ϑ))=0 | (3.41) |
and
ω′(ϑ)+ϵβ/α2d(β−1)/α2ε−1/αq−1/α(ϑ)(∫∞ϑ(∫∞ζ(q1/ε(u))β−1(η∑i=1ϕi(u)gβ(τ−1(δi(u))))dub(ζ))1/αdζ)ω1/α(q(ϑ))=0 | (3.42) |
are oscillatory for some constants λ,ε∈(0,1), then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)), d1,d2>0 and ϵ1,ϵ2∈(0,1).
Proof. Assume that Eq (1.1) has a nonoscillatory solution x(ϑ), say x(ϑ)>0, x(δ(ϑ))>0, and x(τ(ϑ))>0 for ϑ≥ϑ1≥ϑ0.
We first consider what Case (1) holds. By performing the same steps as in the proof of Theorem 3.2, we arrive at (3.27). Since υ(ϑ)>0 and υ′(ϑ)>0, we have limϑ→∞υ(ϑ)≠0 and so by Lemma 2.1, we find (3.32) holds. Using (3.32) in (3.27) gives
(b(ϑ)(υ′′′(ϑ))α)′+ϵβ1dβ−11λ3!q(3(β−1)/ε)+3(ϑ)(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))υ′′′(q(ϑ))≤0. |
With y(ϑ)=b(ϑ)(υ′′′(ϑ))α, we see that y(ϑ) is a positive solution of the differential inequality
y′(ϑ)+ϵβ1dβ−11λ3!q(3(β−1)/ε)+3(ϑ)b1/α(q(ϑ))(η∑i=1ϕi(ϑ)gβ(τ−1(δi(ϑ))))y1/α(q(ϑ))≤0. | (3.43) |
It follows from Lemma 2.4 that the differential equation (3.41) also has a positive solution for all λ1∈(0,1), but this contradicts our assumption on (3.41).
Now, we consider what Case (2) holds. Then again (3.16) holds for every ε∈(0,1). By performing the same steps as in the proof of Theorem 3.2, we arrive at (3.29). Integrating (3.29) from ϑ to ∞, we obtain
−b(ϑ)(υ′′′(ϑ))α+ϵβ2dβ−12υ(q(ϑ))∫∞ϑ(q1/ε(ζ))β−1(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ≤0 |
and so
υ′′′(ϑ)≥ϵβ/α2d(β−1)/α2(1b(ϑ)∫∞ϑ(q1/ε(ζ))β−1(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ)1/αυ1/α(q(ϑ)). | (3.44) |
Integrating (3.44) from ϑ to ∞, we obtain
−υ′′(ϑ)≥ϵβ/α2d(β−1)/α2υ1/α(q(ϑ))∫∞ϑ(∫∞ζ(q1/ε(u))β−1(η∑i=1ϕi(u)gβ(τ−1(δi(u))))dub(ζ))1/αdζ |
and so
υ′′(ϑ)+ϵβ/α2d(β−1)/α2(∫∞ϑ(∫∞ζ(q1/ε(u))β−1(η∑i=1ϕi(u)gβ(τ−1(δi(u))))dub(ζ))1/αdζ)υ1/α(q(ϑ))≤0. | (3.45) |
With ω(ϑ)=υ′(ϑ) and using (3.16) in (3.45), we see that ω(ϑ) is a positive solution of
ω′(ϑ)+ϵβ/α2d(β−1)/α2ε−1/αq−1/α(ϑ)(∫∞ϑ(∫∞ζ(q1/ε(u))β−1(η∑i=1ϕi(u)gβ(τ−1(δi(u))))dub(ζ))1/αdζ)ω1/α(q(ϑ))≤0. |
We finalize the proof using the same method as outlined in Case (1). The proof is now complete.
Corollary 3.2. Let α=1 and β<1 hold. Assume further that there exists a positive function μ∈C1([ϑ0,∞),R) such that (3.1) holds. If
limϑ→∞∫ϑq(ϑ)q(3(β−1)/ε)+3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ=∞ | (3.46) |
and
limϑ→∞∫ϑq(ϑ)q(ϱ)(∫∞ϱ(1b(ζ)∫∞ζ(q1/ε(u))β−1(η∑i=1ϕi(u)gβ(τ−1(δi(u))))du)dζ)dϱ=∞, | (3.47) |
then every solution of (1.1) is oscillatory, where q(ϑ)=τ−1(μ(ϑ)).
Proof. The details of the proof are omitted as they are similar to those of Corollary 3.1.
We provide the following example to demonstrate our results:
Example 3.1. Consider the neutral differential equation
(x(ϑ)+ϑx(ϑA1))(4)+ϕ0ϑ2x(ϑA2)+ϕ0ϑ2x(ϑA3)=0, ϑ≥1, | (3.48) |
It is easy to verify that
∫∞ϑ01b1/α(s)ds=∞. |
Choosing μ(ϑ)=ϑ/A4 and A=max{A2,A3}, where A4>A1, A4≥A and A1,A2,A3,A4>1, then (3.1) holds. We also find that
τ−1(ϑ)=A1ϑ, q(ϑ)=A1A4ϑ, τ−1(δ1(ϑ))=A1A2ϑ andτ−1(δ2(ϑ))=A1A3ϑ. |
Let ε=1/4; we see that (1.4) holds.
Now, it is easy to check that the condition (3.39) is satisfied, where
limϑ→∞∫ϑq(ϑ)q3(ζ)b(q(ζ))(η∑i=1ϕi(ζ)gβ(τ−1(δi(ζ))))dζ =limϑ→∞∫ϑA1ϑ/A4(A1A4s)3(ϕ0s2(A1A2s)−1+ϕ0s2(A1A3s)−1)ds =limϑ→∞∫ϑA1ϑ/A4ϕ0(A1A4)3((A1A2)−1+(A1A3)−1)ds=∞. |
Moreover, we find that condition (3.40) is satisfied, where
limϑ→∞∫ϑq(ϑ)q(ζ)(∫∞ζ1b(u)(∫∞u(η∑i=1ϕi(ϱ)gβ(τ−1(δi(ϱ))))dϱ)du)dζ =limϑ→∞∫ϑA1ϑ/A4A1A4ζ(∫∞ζ(∫∞u(ϕ0ϱ3((A1A2)−1+(A1A3)−1))dϱ)du)dζ =limϑ→∞∫ϑA1ϑ/A4A1A4ζ(∫∞ζ(ϕ02u2((A1A2)−1+(A1A3)−1))du)dζ =limϑ→∞∫ϑA1ϑ/A4A1A4ζ(ϕ02ζ((A1A2)−1+(A1A3)−1))dζ=∞. |
Thus, using Corollary 3.1, every solution of (3.48) is oscillatory.
Remark 3.1. Consider the neutral differential equation
(x(ϑ)+ϑx(ϑ3))(4)+ϕ0ϑ2x(ϑ2)=0, ϑ≥1, | (3.49) |
as a special case of Eq (3.48), we see that Theorem 1.1 cannot be applied to (3.49) since
limϑ→∞g(ϑ)(ϑn−2∫ϑϑ01b1/α(ζ)dζ)γ−1≠0, |
accordingly, Theorem 1.1 fails to study the oscillation of (3.49).
Also, we see that Theorem 1.2 cannot be applied to (3.49) since δ(ϑ)=ϑ/2 is greater than τ(ϑ)=ϑ/3 for ϑ≥1. Accordingly, Theorem 1.2 fails to study the oscillation of (3.49).
Now, by using Corollary 3.1 and choosing μ(ϑ)=ϑ/4 and ε=1/4, we see that the conditions (3.39) and (3.40) are satisfied and therefore, Eq (3.49) is oscillatory.
By using comparison principles, we analyze the asymptotic behavior of solutions to a class of fourth-order neutral differential equations. We have obtained some new oscillation results for (1.1) in the case where ∫∞ϑ1/b1/α(ζ)dζ=∞. These results ensure that all solutions to the studied equation are oscillatory, and they also improve and extend some results from previous studies. It will be of interest to investigate the higher-order differential equations of the form
(b(ϑ)((x(ϑ)+g(ϑ)x(τ(ϑ)))(n−1))α)′+η∑i=1ϕi(ϑ)xβ(δi(ϑ))=0, |
where n≥4 is an even natural number.
The authors contributed equally to this work. Both of the authors have read and approved the final version of the manuscript for publication.
The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors acknowledge the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia, under project Grant No. 6018.
The authors declare no conflicts of interest.
[1] |
J. H. Mcculloch, Measuring the term structure of interest rates, J. Bus., 44 (1971), 19–31. https://doi.org/10.1086/295329 doi: 10.1086/295329
![]() |
[2] | J. H. Mcculloch, The tax-adjusted yield curve, J. Financ., 30 (1975), 811–830. https://doi.org/10.2307/2326860 |
[3] |
O. A. Vasicek, H. G. Fong, Term structure modeling using exponential splines, J. Financ., 37 (1982), 339–348. https://doi.org/10.1111/j.1540-6261.1982.tb03555.x doi: 10.1111/j.1540-6261.1982.tb03555.x
![]() |
[4] |
C. R. Nelson, A. F. Siegel, Parsimonious modeling of yield curves, J. Bus., 60 (1987), 473–489. https://doi.org/10.1086/296409 doi: 10.1086/296409
![]() |
[5] |
D. R. Chambers, W. T. Carleton, D. W. Waldman, A new approach to estimation of the term structure of interest rates, J. Financ. Quant. Anal., 19 (1984), 233–252. https://doi.org/10.2307/2331088 doi: 10.2307/2331088
![]() |
[6] | L. E. O. Svensson, Estimating and interpreting forward interest rates: Sweden 1992–1994, IMF Working Papers, International Monetary Fund, 1994 (1994), 1–76. https://doi.org/10.5089/9781451853759.001 |
[7] |
G. S. Shea, Interest rate term structure estimation with exponential splines: a note, J. Financ., 40 (1985), 319–325. https://doi.org/10.1111/j.1540-6261.1985.tb04952.x doi: 10.1111/j.1540-6261.1985.tb04952.x
![]() |
[8] |
J. M. Steeley, Estimating the gilt-edged term structure: basis splines and confidence intervals, J. Bus. Finan. Account., 18 (1991), 513–529. https://doi.org/10.1111/j.1468-5957.1991.tb01119.x doi: 10.1111/j.1468-5957.1991.tb01119.x
![]() |
[9] | B. H. Lin, D. A. Paxson, Valuing the "new-issue" quality option in bund futures, Rev. Futures Mark., 12 (1993), 349–388. |
[10] |
L. Barzanti, C. Corradi, A note on interest rate term structure estimation using tension splines, Insur. Math. Econ., 22 (1998), 139–143. https://doi.org/10.1016/S0167-6687(97)00035-8 doi: 10.1016/S0167-6687(97)00035-8
![]() |
[11] |
L. Andersen, Discount curve construction with tension splines, Rev. Deriv. Res., 10 (2007), 227–267. https://doi.org/10.1007/s11147-008-9021-2 doi: 10.1007/s11147-008-9021-2
![]() |
[12] |
J. Blomvall, Measurement of interest rates using a convex optimization model, Eur. J. Oper. Res., 256 (2017), 308–316. https://doi.org/10.1016/j.ejor.2016.05.053 doi: 10.1016/j.ejor.2016.05.053
![]() |
[13] |
Y. P. Zhu, Z. B. Chen, X. L. Han, C2 tension splines construction based on a class of sixth-order ordinary differential equations, Bull. Iran. Math. Soc., 48 (2022), 127–150. https://doi.org/10.1007/s41980-020-00505-3 doi: 10.1007/s41980-020-00505-3
![]() |
[14] | W. Sun, Y. X. Yuan, Feasible direction methods, In: Optimization theory and methods, Springer Optimization and Its Applications, Boston: Springer, 1 (2006), 493–531. https://doi.org/10.1007/0-387-24976-1_11 |
[15] |
B. G. Malkiel, E. F. Fama, Efficient capital markets: a review of theory and empirical work, J. Financ., 25 (1970), 383–417. https://doi.org/10.1111/j.1540-6261.1970.tb00518.x doi: 10.1111/j.1540-6261.1970.tb00518.x
![]() |
[16] |
D. G. Schweikert, An interpolation curve using a spline in tension, Stud. Appl. Math., 45 (1965), 312–317. https://doi.org/10.1002/sapm1966451312 doi: 10.1002/sapm1966451312
![]() |
[17] |
S. M. Schaefer, Measuring a tax-specific term structure of interest rates in the market for British government securities, Econ. J., 91 (1981), 415–438. https://doi.org/10.2307/2232594 doi: 10.2307/2232594
![]() |
[18] | D. F. Hendry, Predictive failure and econometric modelling in macroeconomics: the transactions demand for money, In: General-to-specific modelling, Edward Elgar, 2005,535–560. https://doi.org/10.1093/0198293542.003.0012 |
1. | M. Mossa Al-Sawalha, Rasool Shah, Kamsing Nonlaopon, Imran Khan, Osama Y. Ababneh, Fractional evaluation of Kaup-Kupershmidt equation with the exponential-decay kernel, 2023, 8, 2473-6988, 3730, 10.3934/math.2023186 | |
2. | Naveed Iqbal, Muhammad Tajammal Chughtai, Roman Ullah, Fractional Study of the Non-Linear Burgers’ Equations via a Semi-Analytical Technique, 2023, 7, 2504-3110, 103, 10.3390/fractalfract7020103 | |
3. | Aliyu Muhammed Awwal, Adamu Ishaku, Abubakar Sani Halilu, Predrag S. Stanimirović, Nuttapol Pakkaranang, Bancha Panyanak, Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery, 2022, 14, 2073-8994, 2375, 10.3390/sym14112375 | |
4. | M. Mossa Al-Sawalha, Osama Y. Ababneh, Rasool Shah, Nehad Ali Shah, Kamsing Nonlaopon, Combination of Laplace transform and residual power series techniques of special fractional-order non-linear partial differential equations, 2022, 8, 2473-6988, 5266, 10.3934/math.2023264 | |
5. | Saleh Alshammari, M. Mossa Al-Sawalha, Rasool Shah, Approximate Analytical Methods for a Fractional-Order Nonlinear System of Jaulent–Miodek Equation with Energy-Dependent Schrödinger Potential, 2023, 7, 2504-3110, 140, 10.3390/fractalfract7020140 | |
6. | M. Mossa Al-Sawalha, Osama Y. Ababneh, Rasool Shah, Amjad khan, Kamsing Nonlaopon, Numerical analysis of fractional-order Whitham-Broer-Kaup equations with non-singular kernel operators, 2022, 8, 2473-6988, 2308, 10.3934/math.2023120 | |
7. | Naveed Iqbal, Imran Khan, Rasool Shah, Kamsing Nonlaopon, The fuzzy fractional acoustic waves model in terms of the Caputo-Fabrizio operator, 2023, 8, 2473-6988, 1770, 10.3934/math.2023091 | |
8. | Naveed Iqbal, Muhammad Tajammal Chughtai, Nehad Ali Shah, Numerical simulation of fractional-order two-dimensional Helmholtz equations, 2023, 8, 2473-6988, 13205, 10.3934/math.2023667 | |
9. | Yong-Ju Yang, Guo-Li Han, Yu-Zhuo Yuan, Local fractional damped non-linear oscillation: Frequency estimation and energy consumption, 2024, 28, 0354-9836, 2143, 10.2298/TSCI2403143Y | |
10. | Nihar Ranjan Mallick, Snehashish Chakraverty, Rajarama Mohan Jena, Analysis of Large Membrane Vibrations Using Fractional Calculus, 2025, 9, 2504-3110, 219, 10.3390/fractalfract9040219 |