Arxiv:2003.00778V1 [Math.NA] 2 Mar 2020 Rpitsbitdt Elsevier to Submitted Preprint No.SR/WOS-A/PM-20/2018 File Their Vide .Introduction 1

Home , Polynomial sequence

arXiv:2003.00778v1 [math.NA] 2 Mar 2020 rpitsbitdt Elsevier to submitted Preprint No.SR/WOS-A/PM-20/2018 File their vide .Introduction 1. Keywords: f convenient and effective fruitful, the numerically of is accuracy method ensure analysis error and th convergence over, and differentiation of matrix operational s the of furnished utilization then and polynomials Lucas shifted polynomia propose Lucas for formula Rodrigues’ and function weight d order second value initial/boundary of solution for ar nique polynomials Lucas shifted non-familiar work, current In Abstract h hoyo aeesmr oefl ihalo hs prope these of all With powerful. more wavelets of theory the r eie nscin6 h aiiyo u ehdi checke is method section. our last of the validity The 6. section diffe in derived Pantograph are and Lane-Emden for procedure 5, section co In orthogonality their and to polynomials devoted Lucas 2, shifted section about follows: as is structure remaining The yet. wavelets Lucas introduced has Pantogr and tech Lane-Emden new the a of solutions construct for to transformation is work present two of the solutions behind the motivation in interested are we work, present solu the for In direction new a us provide theory wavelets kinds different methods, of solutions the get to ma friendly more are mials, interpolation better and smoothness quality The [1]. pat propagation, wave in applications their to due searcher URL: h aiu rprislk rhgnlt,cmatsuppor compact orthogonality, like properties various The e aee ehdbsdo hfe ua oyoil:Atau A polynomials: Lucas Shifted on based method wavelet New [email protected] hfe ua aee,oeainlmtie egtfunct weight , matrices operational wavelet, Lucas Shifted a colo ahmtc,Cmue nomto cecs Ce Sciences, Information & Computer Mathematics, of School Rks Kumar*), (Rakesh aehKumar* Rakesh K Srivastava) (K. upre yDT iityo cec n ehooy India Technology, and Science of Ministry DST, by Supported a en Koundal Reena , rsligdfeeta equations differential solving or feeta qain.Frti ueia cee ehv d have we scheme, numerical this For equations. ifferential rdc prtoa arxo h hfe ua wavelets. Lucas shifted the of matrix operational product e enrcgiin optrgahc n eia mg tech image medical and graphics computer recognition, tern fdfeeta qain.W a losyta te hnnu than other that say also can We equations. differential of nrdcd ehv osrce opttoa aee t wavelet computational a constructed have We introduced. e dto.I eto ,wvlt n hi rprisaepre are properties their and wavelets 4, section In ndition. itdLcsplnmaspoie ssitdLcswavelet Lucas shifted us provides polynomials Lucas hifted uatrn h aeestcnqebsdo rhgnlpol orthogonal on based technique wavelets the nufacturing p ifrnileutos otebs forkolde no knowledge, our of best the To equations. differential aph in ftedfeeta equations. differential the of tions noeve nLcsplnmas nscin3 ediscussed we 3, section In polynomials. Lucas on overview an ifrn id fscn re ifrnileutos Ou equations. differential order second of kinds different rpsdmto.Ilsrtv xmlsso htteprese the that show examples Illustrative method. proposed ynmrcleprmn nscin7 ocuini oein done is Conclusion 7. section in experiment numerical by d s ute,Lcsplnmasadterpoete r use are properties their and polynomials Lucas Further, ls. ,abtayrglrt n ihodrvnsigmmnsma moments vanishing order high and regularity arbitrary t, eta qaini icse.CnegneadErranaly Error and Convergence discussed. is equation rential iu ae notooa ua oyoil n wavelets and polynomials Lucas orthogonal on based nique ta nvriyo iahlPaeh hrmhl,India Dharamshala, Pradesh, Himachal of University ntral te,wvlt aeatatdtecnieaino vr r every of consideration the attracted have wavelets rties, o,Rdius oml,Tuapproach. Tau formula, Rodrigues’ ion, a .Srivastava K. , a approach hog WOS-A through ac ,2020 3, March eveloped merical sented. nology main r More- We . yno- ech- to d one sis ke nt e- 2. Overview of Lucas polynomials

A. F. Horadam [2] presented the polynomial sequence Ws(Θ) s N 0 deﬁned by the relation { } ∈ ∪{ }

Ws(Θ)= ϑ(Θ)Ws 1(Θ)+ υ(Θ)Ws 2(Θ), s 2,3,4,... , − − ∈{ }

K K K with initial conditions: W0(Θ)= χ0,W1(Θ)= χ1Θ , where ϑ(Θ)= χ2Θ , υ(Θ)= χ3Θ in which χ0, χ1, χ2, χ3 are constants s s and K = 0or1. If W0(Θ)= 2 and W1(Θ)= ϑ(Θ) then Ws(Θ) is expressed as Ws(Θ)= χ (Θ)+℘ (Θ), where χ(Θ)+℘(Θ)=

ϑ(Θ) and χ(Θ)℘(Θ)= υ(Θ). Here, Ws is recognized as Lucas polynomial sequences. The Lucas polynomials are − { } deﬁned by adopting ϑ(Θ)= Θ and υ(Θ)= 1 in Ws . It is given explicitly as follows: { } s s s 2 2 L∗(Θ)= 2− Θ Θ + 4 + Θ + Θ + 4 . (1) s − h p p i Alternatively, the Lucas polynomials Ls∗(Θ) of degree s can also be deﬁned by the following relation [3]:

2sinh(sθ), if s is odd, Ls∗(2sinhθ)=  (2) 2cosh(sθ), if s is even.

 2 3 Some of the Lucas polynomials are as follows: L0∗(Θ)= 2, L1∗(Θ)= Θ, L2∗(Θ)= Θ + 2, L3∗(Θ)= Θ + 3Θ, L4∗(Θ)= Θ4 + 4Θ2 + 2 and so on.

Proposition 2.1. The differential equation

d2y dy (Θ2 + 4) + Θ s2y = 0 (3) dΘ2 dΘ − is satisﬁed by the Lucas polynomial.

Proof. Proof is omitted.

The differential equation given by Eq. (3) is called the Lucas differential equation. From the above discussion, it is not hard to see the following proposition:

Proposition 2.2. If L0∗(Θ)= 2 and L1∗(Θ)= Θ, then three-term recurrence relation of Lucas polynomial Ls∗ is deﬁned as

Ls∗(Θ)= ΘLs∗ 1(Θ)+ Ls∗ 2(Θ), s 2,3,4,... , (4) − − ∈{ }

A sequence of polynomial fs(Θ) is said to be orthogonal on a rectiﬁable Jordan curve Γ if the following relation holds: { } 1 fs(Θ) fh(Θ)ω(Θ) dΘ = δsh, ℓ Γ | | Z where w(Θ) is a weight function deﬁned on Γ and ℓ is the perimeter of the boundary of Γ in the complex plane [4].

Lemma 2.1. The Lucas polynomials Ls∗(Θ) are orthogonal with respect to weight function w(Θ).

2 Proof.

0, if s = h, 6  Ls∗,Lh∗ = π/2, if s = h > 0, h i   π, if s = h = 0.   

Remark 2.1. The Lucas polynomials Ls∗(Θ) have precisely s zeroes in the form [5] (2 j + 1)π Θ j = 2ι cos ; j = 0,1,2,...,s 1. 2s −

Theorem 2.1. (Rodrigues’ Formula). The Lucas polynomials Ls∗(Θ) can be acquired through the Rodrigues’ formula given by:

s s! 2 1/2 d 2 s 1/2 L∗(Θ)= 2 Θ + 4 Θ + 4 − . (5) s (2s)! dΘs n o Proof. Proof is omitted.

Proposition 2.3. The generating function for Lucas polynomial is deﬁned as [6]

∞ 2 Θt L (Θ)th = − . ∑ s∗ 1 Θt t2 s=0 − −

3. Shifted Lucas polynomials

The shifted Lucas polynomials are equally important as Lucas polynomials. Few shifted Lucas polynomials are given below:

Q0∗(t)= 2,

Q∗(t)= 2t 2ι, 1 − 2 Q∗(t)= 4t 8ιt 2, 2 − − 3 2 Q∗(t)= 8t 24ιt t 18t + 2ι, 3 − − . .

Also it can be easily seen that the orthogonality condition for the shifted Lucas polynomials is as follows:

0, if s = h, 6 Qs∗,Qh∗ =  h i παs  2 , if s = h, where 

2, if s = 0, αs =  (6) 1, if s 1.  ≥  3 4. Lucas wavelets and shifted Lucas wavelet

Analogous to Legendre wavelets, Chebyshev wavelets , Lucas wavelets also have ϕh,s(Θ)= ϕ(k,h,s,Θ) four arguments. On the interval [0,2), the Lucas wavelets family of functions can be deﬁned as

k+1 hˆ 2 hˆ + 2 2 2 L (2k+1ιΘ hˆι), − Θ < , ∗s − 2k+1 ≤ 2k+1 ϕh,s(Θ)=  (7) 0, e o.w, where  1 , s = 0, √π L (Θ)=  (8) ∗s 2  L (Θ), s > 0,  π s∗ e r  s = 0,1,2,...,S 1; S is the maximum order of the Lucas polynomials; hˆ = 2(2h + 1) where h = 0,1,2,...,2k 1; k is any − − nonnegative integer and Ls∗(Θ) are Lucas polynomials of order s. Also, On the interval [0,2] shifted Lucas wavelet is deﬁned as

k+1 2 h h + 1 2 k 2 Qs∗(2 ιΘ 2hι), k 1 Θ < k 1 , παs − 2 ≤ 2 ϕh,s(Θ)=  r − − 0, o.w,

k  where h = 0,1,2,...,2 1 and αs is deﬁned in previous deﬁnitions. −

5. Shifted Lucas tau wavelet technique for Lane-Emden and Pantograph differential equations

A function ρ(Θ) L2 [0,2] can be deﬁned as a Lucas wavelet series in the following form ∈ wh ∞ ∞ E ρ(Θ)= ∑ ∑ h,sϕh,s(Θ), (9) h=0 s=0 where

2 Eh,s = ρ(Θ),ϕh,s(Θ) = ρ(Θ)ϕh,s(Θ)wh(Θ)dΘ. (10) wh(Θ) Z0

In equation (10), .,. denotes the inner product. h i The truncated Lucas wavelet series can be deﬁned as

2k 1 S 1 − − E T ρ(Θ)= ∑ ∑ h,sϕh,s(Θ)= E Ψ(Θ). (11) h=0 s=0

Here, E and Ψ(Θ) are 2kS 1 matrices. × Theorem 5.1. If Ψ(Θ) is Lucas wavelets vector then,

dΨ(Θ) = DΨ(Θ), (12) dΘ

4 be the the Ist derivative of the vector Ψ(Θ), where, D is 2kS 2kS matrix which is deﬁned by × F 0 D = (13) 0 F   Also, F is a square matrix of order S whose (ϒ,λ)th element is deﬁned by

k+1 αϒ 1 2 ι(ϒ 1) − , ϒ = 2,...,S, λ = 1,...,ϒ 1 and (ϒ + λ)odd, αλ 1 F = − − − (14) ϒ,λ  q 0, o.w.

Proof. The proof of this theorem can be obtain by using previous deﬁnitions.

Corollary 5.1. The differential operational matrix for hth-order can be attained by utilizing (12) as

d(h)Ψ(Θ) = D(h)Ψ(Θ), (15) dΘ(h) where D(h) is hth-order of the matrix D.

Further, we will introduce product operational matrix, which is utilized for solving the differential equations

T ˜ Ψ(Θ)Ψ (Θ)E ∼= EΨ(Θ). (16)

We consider the second order differential equation for the foundation of shifted Lucas wavelet technique as

ρ ′′ (Θ)= G(Θ,ρ(Θ),ρ(αΘ)), Θ (0,l), (17) ∈ with supplementary conditions

ρ(Θ) Θ 0 = A1, ρ ′ (Θ) Θ 0 = A2, (18) | = | =

ρ(Θ) Θ 0 = B1, ρ ′ (Θ) l = B2. (19) | = |Θ= In order to apply the shifted orthogonal Lucas wavelet technique, ﬁrstly we have to approximate the unknown function ρ(Θ) and ρ ′ (Θ), ρ ′′ (Θ) in equation (17), (18) and (19 and construct a residual function as:

ℜ(Θ)= ET D2Ψ(Θ) G Θ,ET Ψ(Θ),ET Ψ(αΘ) . (20) − k Considering the Tau method, equation (20) generates the (2 S 2) nonlinear equations with the unknown coefﬁcients En,m. − The supplementary conditions are approximated as:

T T ρ(Θ) Θ 0 = E Ψ(0)= A1, ρ ′ (Θ) Θ 0 = E DΨ(0)= A2, (21) | = | =

T T ρ(Θ) Θ 0 = E Ψ(0)= B1, ρ ′ (Θ) l = E DΨ(0)= B2, (22) | = |Θ= Thus, approximated initial and boundary conditions and (2kS 2) nonlinear equation generates set of 2kS nonlinear equations. − Further, these equations can be solved to determine the unknown component of vector E. Hence ρ(Θ) can be calculated.

5 6. Convergence and error analysis

Theorem 6.1. The Lucas wavelet series solution deﬁned in equation (9) converges towards ρ(Θ) by using Lucas wavelet.

Proof. proof is omitted.

Theorem 6.2. Suppose that ρ(Θ) L2 [0,2], as well as bounded second order derivative, say ρ (Θ) 6 N, then the function ∈ wh | ′′ | ρ(Θ) can be elaborated as an infinite sum of Lucas wavelets and the series is uniformly convergent to ρ(Θ). Additionally, the following inequalities are satisfied by the coefficients Eh,s defined in (10): 2N√π Eh,s , s > 1, h 0, (23) ≤ (h + 1)5/2(s2 1) ∀ ≥ − √π and Eh,1 max ρ′(Θ) , for s = 1, h 0. (24) ≤ (h + 1)3/2 0 Θ 2 ∀ ≥ ≤ ≤

Proof. Proof is omitted.

Lemma 6.1. Let Θ > s, then f (Θ) will be a continuous, positive and decreasing function. If f (k)= Ek along with ∑Eh, is ∞ ∞ convergent and ℜh = ∑ Ek, then ℜh f (Θ)dΘ [7]. k=h+1 ≤ h R Theorem 6.3. If ρ L2 deﬁned in (9) satisﬁes theorem (6.1) and ∈ wh 2k 1 S 1 − − E ρk,S(Θ)= ∑ ∑ h,sϕh,s(Θ) h=0 s=0 then

2 2 2 1 (S 2S)ln(S) S ln(S 2) + (2ln(S 2) 2)S + 2 ρ ρk,S w < N π − − − − − k − k h 25(2k 1) 25ln(2) 4S(S 2) s − − − is the error estimation, for S > 2.

Proof. Proof is omitted.

7. Illustrative examples

In this section, we consider the following test problems to show the validity and accuracy of the proposed scheme. Problem 1. Let us consider the nonlinear Lane-Emden problem [8] 6 ρ ′′ (Θ)+ ρ ′ (Θ)+ 14ρ(Θ)= 4ρ(Θ)ln(ρ(Θ)), 0 < Θ 1, (25) Θ − ≤ with

ρ(Θ) Θ 0 = 1, ρ ′ (Θ) Θ 0 = 0. (26) | = | =

Θ2 This problem has exact solution ρ(Θ)= e− . First of all to solve equation (25), we are used same transformation as we have applied in above example i.e. ρ(Θ) = ez(Θ). Then transformed form of equations (25) and (26) are 6 z′′ (Θ) + (z′ (Θ))2 + z′ (Θ)+ 14 = 4z(Θ), Θ 0 (27) Θ − ≥

6 with

z(Θ) Θ 0 = 0, z′ (Θ) Θ 0 = 0. | = | =

Also,

z(Θ)= Θ2. −

The values of component vector E for equation (27) are formed by utilizing the method described in section (5). In this case, we take k = 0 and S = 3. After some manipulations we get

z(Θ) 5.55112 10 17 Θ2 ρ(Θ)= e = e × − − , which is the approximate result for considered problem. Problem 2. Let us consider the linear Pantograph initial value problem [9]

3 Θ ρ ′′ (Θ) ρ(Θ) ρ + Θ2 2 = 0, 0 Θ 1, (28) − 4 − 2 − ≤ ≤ with

ρ(Θ) Θ 0 = ρ ′ (Θ) Θ 0 = 0. (29) | = | =

This problem has exact solution ρ(Θ)= Θ2 at γ = 2 given in [9]. In the proposed technique let k = 0 and S = 3. After performing some calculations we get

T 16 17 ρ(Θ)=E Ψ(Θ) = ( 4.996003610813204 10− + 8.088790061448873 10− ι) − × × 16 17 + (2.220446049250313 10− + 2.7193641226633042 10− ι)Θ × × 16 2 + (0.9999999999999994 4.8544160607900805 10− ι)Θ . − × which is the approximate solution for considered problem.

8. Conclusion

In the present article, we have developed a novel technique for numerical solutions of differential equations. The ba- sic foundations for this technique are shifted Lucas polynomials and Tau method. The shifted Lucas polynomials are well behaved orthonormal basis functions. To show the validity of NSLW technique, we have considered two different kind of differential equations. We hope that, the initiative of introducing the shifted Lucas polynomials will ﬁll the gap of liter- ature, related to Lucas polynomials. The computations associated with all the considered examples are done by software Mathematica 9.

References

[1] A. Constantmldes. Applied Numerical Methods with Personal Computers, McGraw-Hill, New York, NY, USA, 1987.

7 [2] A.F. Horadam. A Synthesis of Certain Polynomial Sequences in: Application of Fibonacci Number, vol. 6, Kluwer Academic Publishers, 1996, pp.215-229.

[3] A. Lachal. A trick around Fibonacci, Lucas and Chebyshev. arXiv preprint arXiv:1302.0357 (2013).

[4] G. Szeg¨o. Orthogonal polynomials, American Mathematical Society, Providence, R.I. 1975.

[5] K. Thomas. Fibonacci and Lucas numbers with applications. Vol. 1. John Wiley & Sons, 2017.

[6] A. Boussayoud, M. Kerada, N. Harrouche.On the k-Lucasnumbers and Lucas polynomials. Turkish Journal of Analysis and Number 5 (2017) 121-125.

[7] J. Stewart. Single Variable Essential Calculus: Early Transcendentals (Cengage Learning, 2012).

[8] H. Aminikhah, S. Moradian. Numerical Solution of Singular Lane-Emden Equation. Mathematical Physics Volume 2013,Article ID 507145, 9 pages.

[9] P. Rahimkhani, Y. Ordokhani, E. Babolian. Numerical solution of fractional pantograph differential equations by using generalized fractional-order Bernoulli wavelet. Journal of Computational and Applied Mathematics, 309 (2017), 493- 510.