This paper deals with the design fractional solution of Bessel equation. We obtain explicit solutions of the equation with the help of fractional calculus techniques. Using the N-fractional calculus operator Nν method, we derive the fractional solutions of the equation.
1. Introduction, Definitions, and Preliminaries
Fractional calculus has an important place in the field of math. Firstly, L’Hospital and Leibniz were interested in the topic in 1695, [1]. Fractional calculus is an area of applied mathematics that deals with derivatives and integrals of arbitrary orders and their applications in science, engineering, mathematics, economics, and other fields. The seeds of fractional derivatives were planted over 300 years ago. Since then many efficient mathematicians of their times, such as N. H. Abel, M. Caputo, L. Euler, J. Fourier, A. K. Grunwald, J. Hadamard, G. H. Hardy, O. Heaviside, H. J. Holmgren, P. S. Laplace, G. W. Leibniz, A. V. Letnikov, J. Liouville, and B. Riemann, have contributed to this field; all these references can be seen in [1–5]. The mathematics involved appeared very different applications of this field. Fractional calculus has been applied to almost every field of science. They are viscoelasticity, electrical engineering, electrochemistry, biology, biophysics and bioengineering, signal and image processing, mechanics, mechatronics, physics, and control theory. During the last decade, Samko et al. [4], Nishimoto [6–12], and Podlubny [3] have been helpful in introducing the field to engineering, science, economics and finance, and pure and applied field. Furthermore, there were many studies in this field [5, 13–16]. Various scientists have studied that concept. The progress in this field continues [2–4, 6–12, 17–20].
N-Fractional calculus is a very interesting method because this method is applied to singular equation. Note that fractional solutions can be obtained for kinds of singular equation via this method [6–12, 17]. In this paper, our aim is to apply the same way for singular Sturm-Liouville equation with Bessel potential and find fractional solutions of this equation. Furthermore, we give some applications and their graphs of fractional solutions of the equation.
Now, consider the following the Bessel equation:
(1)xd2zdx2+dzdx+[λx-p2x]z=0,0<x≤1,
where λ and p are real numbers. By means of the substitution y=xz (1) reduces to the form
(2)d2ydx2+(λ-p2-1/4x2)y=0.
Bessel equation for having the analogous singularity is given in [21].
The differintegration operators and their generalizations [6–11, 17, 18] have been used to solve some classes of differential equations and fractional differential equations.
Two of the most commonly encountered tools in the theory and applications of fractional calculus are provided by the Riemann-Liouville operator Rzυ(υ∈ℂ) and the Weyl operator Wzυ(υ∈ℂ), which are defined by [17, 19].
provided that the defining integrals in (3) and (4) exist, ℕ being the set of positive integers.
Definition 1 (cf. [6–10, 12, 20]).
Let
(5)D={D-,D+},C={C-,C+},
where C- is a curve along the cut joining two points z and -∞+iIm(z), C+ is a curve along the cut joining two points z and ∞+iIm(z), D- is a domain surrounded by C-, and D+ is a domain surrounded by C+. (Here D contains the points over the curve C.)
Moreover, let f=f(z) be a regular function in D(z∈D),
(6)fν(z)=(f(z))ν=Γ(ν+1)2πi∫Cf(t)dt(t-z)ν+1(ν∈ℝ∖ℤ-;ℤ-={…,-3,-2,-1}),f-n(z)=limν→-nfν(z)(n∈ℤ+),
where t≠z,
(7)-π≤arg(t-z)≤πforC-,0≤arg(t-z)≤2πforC+.
Then fν(z)(ν>0) is said to be the fractional derivative of f(z) of order ν and fν(z)(ν<0) is said to be the fractional integral of f(z) of order -ν, provided (in each case) that |fν(z)|<∞(ν∈ℝ).
Finally, let the fractional calculus operator (Nishimoto’s operator) Nν be defined by (cf. [6–10])
(8)Nν=(Γ(ν+1)2πi∫Cdt(t-z)ν+1)(ν∉ℤ-)
with
(9)N-n=limν→-nNν(n∈ℤ+).
We find it to be worthwhile to recall here the following useful lemmas and properties associated with the fractional differintegration which was defined earlier (cf. e.g., [6–10, 12]).
Lemma 2 (linearity property).
If the functions f(z) and g(z) are single-valued and analytic in some domain Ω⊆ℂ, then
(10)(h1f(z)+h2g(z))ν=h1fν(z)+h2gν(z)(ν∈ℝ;z∈Ω)
for any constants h1 and h2.
Lemma 3 (index law).
If the function f(z) is single-valued and analytic in some domain Ω⊆ℂ, then
(11)(fσ(z))ν=fσ+ν(z)=(fν(z))σ(fσ(z)≠0;fν(z)≠0;σ,ν∈ℝ;z∈Ω).
Lemma 4 (generalized Leibniz rule).
If the functions f(z) and g(z) are single-valued and analytic in some domain Ω⊆ℂ, then
(12)(f(z)·g(z))ν=∑n=0∞(νn)fν-n(z)·gn(z)(ν∈ℝ;z∈Ω),
where gn(z) is the ordinary derivative of g(z) of order n(n∈ℕ0:=ℕ∪{0}), being tacitly assumed (for simplicity) that g(z) is the polynomial part (if any) of the product f(z)g(z).
Property 1.
For a constant λ,
(13)(eλz)ν=λνeλz(λ≠0;ν∈ℝ;z∈ℂ).
Property 2.
For a constant λ,
(14)(e-λz)ν=e-iπνλνe-λz(λ≠0;ν∈ℝ;z∈ℂ).
Property 3.
For a constant λ,
(15)(zλ)ν=e-iπνΓ(ν-λ)Γ(-λ)zλ-ν(ν∈ℝ;z∈ℂ;|Γ(ν-λ)Γ(-λ)|<∞).
Now, let apply N-fractional method to nonhomogeneous Bessel equation.
2. The Nν-Method Applied to Bessel EquationTheorem 5.
Let y∈{y:0≠|yν|<∞;ν∈ℝ} and f∈{f:0≠|fν|<∞;ν∈ℝ}. We consider the nonhomogeneous Bessel equation:
(16)L[y,x,λ,p]=y2+y[λ-p2-1/4x2]=f(0<x≤1),
and it has particular solutions of the forms
(17)yı=xp+1/2e-iλx×{[(fx1/2-peiλx)-p-1/2e-2iλxxp-1/2]-1×e2iλxx-p-1/2∑}p-1/2,(18)yıı=xp+1/2eiλx×{[(fx1/2-pe-iλx)-p-1/2e2iλxxp-1/2]-1×e-2iλxx-p-1/2∑}p-1/2,(19)yııı=x-p+1/2e-iλx×{[(fx1/2+peiλx)p-1/2e-2iλxx-p-1/2]-1×e2iλxxp-1/2∑}-p-1/2,(20)yıv=x-p+1/2eiλx×{[(fx1/2+pe-iλx)p-1/2e2iλxx-p-1/2]-1×e-2iλxxp-1/2∑}-p-1/2,
where y2=d2y/dx2, y=y(z)(z∈ℂ), f=f(z) (an arbitrary given function), and p,λ are given constants.
Remark 6.
The cases p=0 of (19) and (20) coincide with those (17) and (18).
Proof.
Set
(21)y=xηψ,ψ=ψ(x).
Thus
(22)y1=ηxη-1ψ+xηψ1,y2=η(η-1)xη-2ψ+2ηxη-1ψ1+xηψ2.
Putting (21) and (22) in (16), we obtained
(23)ψ2xη+ψ12ηxη-1+ψη(η-1)xη-2+xηψ[λ-p2-1/4x2]=f.
With some rearrangement of the terms in (23), we have
(24)ψ2x+ψ12η+ψ[x-1(η2-η+14-p2)+λx]=fx1-η.
Here, we choose η such that
(25)η2-η+14-p2=0.
That is,
(26)η=12±p.
(I) Let η=p+(1/2). From (21) and (24), we have
(27)y=xp+(1/2)ψ,(28)ψ2x+ψ1(2p+1)+ψλx=fx(1/2)-p.
Set
(29)ψ=eμxϕ,ϕ=ϕ(x).
Rewrite (28) in the form
(30)(eμxϕ)2x+(eμxϕ)1(2p+1)+eμxϕλx=fx(1/2)-p.
At this point, differentiating eμxϕ two times,
(31)(eμxϕ)1=eμx(μϕ+ϕ1),(eμxϕ)2=eμx(μ2ϕ+2μϕ1+ϕ2),
and substituting from (29) and (31) in (30), we can express (30) as
(32)ϕ2x+ϕ1(2μx+2p+1)+ϕ(2pμ+μ+x(μ2+λ))=fx(1/2)-pe-μx.
Choose μ such that
(33)μ2+λ=0.
That is,
(34)μ=±λi.
(I) (i): For instance, taking μ=-λi, we have
(35)ψ=e-λixϕ,(36)ϕ2x+ϕ1[-λi2x+2p+1]+ϕ[-i(2p+1)λ]=fx(1/2)-pe-μx
from (29) and (32).
Applying the operator Nν to both members of (36), we find the following equality:
(37)[ϕ2x]ν+{ϕ1[-λi2x+2p+1]}ν+{ϕ[-i(2p+1)λ]}ν=[fx(1/2)-pe-μx]ν.
Using (3)–(12), we have
(38)[ϕ2x]ν=xϕ2+ν+νϕ1+ν,{ϕ1[-λi2x+2p+1]}ν=ϕ1+ν[-λi2x+2p+1]-λi2νϕν.
Making use of the relations (38), rewrite (37) in the following form:
(39)ϕ2+νx+ϕ1+ν[-λi2x+2p+1+ν]-ϕν[λi2ν+i(2p+1)λ]=[fx(1/2)-peiλx]ν.
Choose ν such that
(40)ν=-p-12.
We then have
(41)ϕ-p+3/2x+ϕ-p+1/2[p+12-2λix]=[fx(1/2)-peiλx]-p-1/2
from (39).
Next, writing
(42)ϕ-p+1/2=ω=ω(x),
we obtain the following equality from (41):
(43)ω1+ω[x-1(p+12)-2λi]=[fx(1/2)-peiλx]-p-1/2x-1.
This is an ordinary differential equation of the first order which has a particular solution,
(44)ω=[[fx(1/2)-peiλx]-p-1/2x-1e-2λixxp+1/2]-1×e2λixx-p-1/2.
Making use of the reverse process to obtain yı, we finally obtain the solution (17) from (44), (42), (35), and (27).
Inversely, (44) satisfies (43); then
(45)ϕ=ωp-1/2
satisfies (41). Therefore, (17) satisfies (16) because we have (27), (35), (44), and (45).
(I) (ii): In the case when μ=λi, we have
(46)ψ=eλixϕ,(47)ϕ2x+ϕ1[2λix+2p+1]+ϕ[(2p+1)λi]=fx(1/2)-pe-iλx
from (29) and (32).
Applying the operator Nν to both members of (47), we have
(48)ϕ2+νx+ϕ1+ν[2λix+2p+1+ν]+ϕν[λi(2p+1)i]=(fx(1/2)-pe-iλx)ν.
Choosing ν such that
(49)ν=-p-12
and replacing
(50)ϕ-p+1/2=ϑ=ϑ(x),
we then obtain
(51)ϑ1+ϑ[2λi+(p+12)x-1]=(fx(1/2)-pe-iλx)-p-1/2x-1
from (48). A particular solution of (51) is given by
(52)ϑ=[(fx1/2-pe-iλx)-p-1/2x-1e2λixxp+1/2dx]-1×e-2λixx-p-1/2.
Thus, we have (18) from (52), (50), (46), and (27).
(II) Let η=-p+(1/2).
With the help of the similar method in (I), replacing p by -p in (I) (i) and (I) (ii), we have other solutions (19) and (20) different from (17) and (18), respectively, if p≠0.
3. The Operator Nν-Method to a Homogeneous Bessel EquationTheorem 7.
If y∈℘,o just as in Theorem 5, then the homogeneous Bessel equation
(53)L[y,x,p]=y2+y[λ-p2-1/4x2]=0,(0<x≤1),
has solutions of the forms
(54)yı=kxp+1/2e-iλx{e2iλxx-p-1/2}p-1/2,(55)yıı=kxp+1/2eiλx{e-2iλxx-p-1/2}p-1/2,(56)yııı=kx-p+1/2e-iλx{e2iλxxp-1/2}-p-1/2,(57)yıv=kx-p+1/2eiλx{e-2iλxxp-1/2}-p-1/2,
for p≠0, where k is an arbitrary constant.
Remark 8.
In the case when p=0,(56) and (57) coincide with (54) and (55).
Proof.
When f=0 in Section 2, we have
(58)ω1+ω[(p+12)x-1-2λi]=0,(59)ϑ1+ϑ[(p+12)x-1+2λi]=0,
for μ=-iλ and μ=iλ, instead of (43) and (51).
Therefore, we get (54) for (58) and (55) for (59).
And, for η=-p+(1/2), replacing p by -p in (58) and (59), we have (56) and (57).
Theorem 9.
Let y∈℘o and f∈℘o just as in Theorem 5. Then the nonhomogeneous modified Sturm-Liouville equation (16) is satisfied by the fractional differintegrated functions
(60)y=yı+yı.
Proof.
It is clear by Theorems 5 and 7.
Application 1. If we substitute p=0,λ=1/4, and f=ix-1/2e-(i/2)x in (16), then we obtain the following equation:
(61)y2+y(14+14x2)=ix-1/2e-(i/2)x,
and its solution is
(62)y=x1/2e(-i/2)x{[(ix-1/2e-(i/2)xx1/2e(i/2)x)-1/2×e-ixx-1/2]-1eixx-1/2}-1/2.
By performing the necessary operations in (62), we get
(63)y=x1/2e(-i/2)x{[2ixπe-ixx-1/2]-1eixx-1/2}-1/2,
where Riemann Liouville operator is
(64)[i]-1/2=1Γ(1/2)∫0xix-tdt=2ixπ,y=x1/2e(-i/2)x(-2x-1/2π)-1/2,
and using the definitions of Riemann Liouville operator again, we obtain the following solution:
(65)(-2x-1/2π)-1/2=-1Γ(1/2)∫0x2t-1/2πx-tdt=-2,(66)y=-2x1/2e-(i/2)x.
Now, let us show that the last equality is the solution of (61):
(67)y2=e-(i/2)x[12x1/2+12x-3/2+ix-1/2].
Obviously, if (66) and (67) are put in (61), it is satisfied. The graph of the solution of (61) is given in Figures 1 and 2.
Rey.
Imy.
Application 2. If we substitute p=-1 and λ=0 in (53), then we obtain the following equation:(68)y2-34x2y=0(0<x≤1),
and its solution is
(69)y=kx-1/2{x1/2}-3/2.
We prove that yı is the solution of (67). With the help of Riemann Liouville operator,
(70)[x-1/2]-1/2=1Γ(3/2)∫0xt1/2x-tdt=πx24,(71)y=kx3/2π4.
Now, let us show that the last equality is the solution of (67),
(72)y2=3kx-1/2π16.
Obviously, if (71) and (72) are put into (68), it is satisfied. The graph of the solution of (68) is given in Figure 3.
4. Two Further Cases of Modified Bessel EquationTheorem 10.
In the similar way as in the previous sections, we can solve the following nonhomogeneous modified Bessel equation:
(73)y2+y[λ+(1/4)+p2x2]=f,y2+y[-λ+(1/4)+p2x2]=f,
which are obtained by replacing p by ip (-λ instead of λ) in (16); that is,
(74)y2+y[λ+(1/4)-(ip)2x2]=f,(75)y2+y[-λ+(1/4)-(ip)2x2]=f.
(i) Therefore, the solutions for (74) are given by replacing p by ip in (17), (18), (19), and (20) as follows:(76)y(ı)=xpi+1/2e-iλx×{[(fx1/2-pieiλx)-pi-1/2e-2iλxxpi-1/2]-1×e2iλxx-pi-1/2∑}pi-1/2,y(ıı)=xpi+1/2eiλx×{[(fx1/2-pie-iλx)-pi-1/2e2iλxxpi-1/2]-1×e-2iλxx-pi-1/2∑}pi-1/2,y(ııı)=x-pi+1/2e-iλx×{[(fx1/2+pieiλx)pi-1/2e-2iλxx-pi-1/2]-1×e2iλxxpi-1/2∑}-pi-1/2,y(ıv)=x-pi+1/2eiλx×{[(fx1/2+pie-iλx)pi-1/2e2iλxx-pi-1/2]-1×e-2iλxxpi-1/2∑}-pi-1/2.
(ii) In the same way, for the solutions for (75), substituting the relations (21), and (22) into (75), we have(77)ϕ2x+ϕ12ν+ϕ[(ν2-ν+14-(pi)2)x-1-λx]=fx1-ν.
Choose ν as follows:
(78)ν2-ν+14+p2=0.
That is
(79)ν=12±pi.
Let ν=pi+(1/2). From (21) and (77), we have
(80)y=xpi+(1/2)ϕ,(81)ϕ2x+ϕ1(2pi+1)-ϕλx=fx(1/2)-pi.
Next, set (29); then (81) is rewritten in the form
(82)(eμxψ)2x+(eμxψ)1(2ip+1)-eμxψλx=fx(1/2)-pi.
Substituting the relations (29) and (31) into (82), we have
(83)ψ2x+ψ1(2μx+2pi+1)+ψ[(μ2-λ)x+(2pi+1)μ]=fx(1/2)-ipe-μx.
Choose μ as follows:
(84)μ2-λ=0.
That is,
(85)μ=±λ.
(ii. 1) In the case when μ=-λ, we have
(86)ϕ=e-λxψ,(87)ψ2x+ψ1(-2λx+2pi+1)-ψ[λ(2pi+1)]=fx(1/2)-pieλx
from (29) and (83).
Applying the operator Nν to both members of (87), we then obtain
(88)(ψ2x)ν+[ψ1(-2λx+2pi+1)]ν+{ψ[-λ(2pi+1)]}ν=(fx(1/2)-pieλx)ν.
Using (4), (10), and (11), we have
(89)ψ2+νx+ψ1+ν(-2λx+2pi+1+ν)+ψν[-λ(2pi+1+2v)]=(fx(1/2)-pieλx)ν.
Choose ν such that
(90)ν=-pi-12.
We then have
(91)ψ2-pi-1/2x+ψ1-pi-1/2[-2λx+pi+12]=(fx(1/2)-pieλx)-pi-1/2
from (89).
Next, writing
(92)ψ1/2-pi=u=u(x),
we obtain
(93)u1+u[-2λ+(pi+12)x-1]=(fx(1/2)-pieλx)-pi-1/2x-1
from (91). This is an ordinary differential equation of the first order which has a particular solution
(94)u=[(fx(1/2)-pieλx)-pi-1/2e-2λxxpi-1/2]-1×e2λxx-pi-1/2.
We finally obtain the solution
(95)y(ı)=xpi+1/2e-λx×{[(fx1/2-pieλx)-pi-1/2e-2λxxpi-1/2]-1×e2λxx-pi-1/2∑}pi-1/2.
from (94), (92), (86) and (80).
(ii. 2) Similarly, in the case when μ=λ, we obtain
(96)y(ıı)=rpi+1/2eλx×{[(fx1/2-pie-λx)-pi-1/2e2λxxpi-1/2]-1×e-2λxx-pi-1/2∑}pi-1/2.
Let ν=-ip+(1/2). In the same way as in the procedure in (ii), replacing ip by –ip (ii. 1) and (ii. 2), we can obtain y(ııı) and y(ıv).
Theorem 11.
In the homogeneous case for (74) with f=0, using the solutions (54), (55), (56), and (57) and replacing p by pi, we obtain
(97)y(ı)=αxpi+1/2e-iλx(e2iλxx-pi-1/2)pi-1/2,y(ıı)=αxpi+1/2eiλx(e-2iλxx-pi-1/2)pi-1/2,y(ııı)=αx-pi+1/2e-iλx(e2iλxxpi-1/2)-pi-1/2,y(ıv)=αx-pi+1/2eiλx(e-2iλxxpi-1/2)-pi-1/2,
for p≠0, where α is an arbitrary constant.
5. Conclusion
The N-fractional calculus operator Nν-method is applied to the nonhomogeneous and homogeneous Bessel equation. Explicit fractional solutions of Bessel equations are obtained. Furthermore, similar solutions were obtained for the modified same equation by using the method.
Acknowledgment
The authors sincerely thank the reviewers for their valuable suggestions and useful comments.
AdamL.2004KilbasA. A.SrivastavaH. M.TrujilloJ. J.2006204Amsterdam, The NetherlandsElsevier (North-Holland) Science PublishersNorth-Holland Mathematical StudiesPodlubnyI.1999198New York, NY, USAAcademic PressSamkoS. G.KilbasA. A.MarichevO. I.1993Reading, UKGordon and Breach Science Publishers“Nauka i Tekhnika”, Minsk, (1987)OldhamK. B.SpanierJ.1974New York, NY, USAAkademic PressNishimotoK.1984IKoriyama, JapanDescartes PressNishimotoK.1987IIKoriyama, JapanDescartes PressNishimotoK.1989IIIKoriyama, JapanDescartes PressNishimotoK.1991IVKoriyama, JapanDescartes PressNishimotoK.1996VKoriyama, JapanDescartes PressNishimotoK.de RomeroS. S.N-fractional calculus operator Nν method to nonhomogeneous and homogeneous Whittaker Equations199691722NishimotoK.1991Koriyama, JapanDescartes PressCozM.RochusP.Translation kernels for velocity dependent interactions19761811223222402-s2.0-33744796368ErdelyA.MagnusW.OberhettingerF.TricomiF. G.1954IINew York, NY, USAMcGraw-HillErdelyA.Axially symmetric potentials and fractional integration196513216228WeinsteinA.Generalized axially symmetric potential theory1953592038RossB.1975New York, NY, USASpringerConference Proceedings held at the University of New Haven, June 1974Salinas De RomeroS.SrivastavaH. M.An application of the N-fractional calculus operator method to a modified Whittaker equation2000115111212-s2.0-0039592686SrivastavaH. M.BuschmanR. G.199279Dordrecht, The NetherlandsKluwer Academic PublishersKluwer Series on Mathematics and Its ApplicationsSrivastavaH. M.OwaS.NishimotoK.Some fractional differintegral equations198510623603662-s2.0-0001032572KoyunbakanH.PanakhovE. S.Transformation operator for singular Sturm-Liouville equations200314135143