arxiv: v1 [math.na] 23 Apr 2018

Transcription

1 rxiv: v mth.na] 23 Apr 208 Solving generlized Abel s integrl equtions of the first nd second kinds vi Tylor-colloction method Eis Zrei, nd Smd Noeighdm b, Deprtment of Mthemtics, Hmedn Brnch, Islmic Azd University, Hmedn, Irn. b Deprtment of Mthemtics, Centrl Tehrn Brnch, Islmic Azd University, Tehrn, Irn. Abstrct The im of this pper is to present n efficient numericl procedure to pproximte the generlized Abel s integrl equtions of the first nd second kinds. For this reson, the Tylor polynomils nd the colloction method re pplied. Also, the error nlysis of presented method is illustrted. Severl exmples re pproximted nd the numericl results show the ccurcy nd efficiency of this method. keywords: Generlized Abel s integrl eqution, Colloction method, Tylor polynomils. Introduction The gret mthemticin Niels Abel, gve the inititive of integrl equtions in 823 in his study of mthemticl physics 3, 32, 33, 34]. Zeilon in 924, studied the generlized Abel s integrl eqution on finite segment 38]. The Abel s integrl equtions re the singulr form of Volterr integrl equtions. The singulr integrl equtions re the importnt nd pplicble kinds of integrl equtions which solved by mny uthors 9, 0,, 24, 28, 34]. These integrl equtions hve mny pplictions in vrious res such s simultneous dul reltions 29], stellr winds 8], wter wve 4], spectroscopic dt 3] nd the others 2, 3, 9]. In this pper, we investigte the following generlized Abel s integrl equtions of the first kind nd the second kind x (φ(x) φ(t)) αdt = g(x), 0 < α <, () x Φ(x) = g(x)+ (φ(x) φ(t)) αdt, 0 < α <, (2) Corresponding uthor, Tel.: , E-mil ddress: ei.zrei353@gmil.com. E-mil ddresses: s.noeighdm.sci@iuctb.c.ir, smdnoeighdm@gmil.com.

2 where is given rel vlue, g(x) nd φ(x) re known functions nd Φ(x) is n unknown function tht φ(t) is strictly monotoniclly incresing nd differentible function in some intervl < t < b, nd φ (t) 0 for every t in the intervl. Abel s integrl equtions() nd(2) hve pplictions in experimentl physics such s plsm dignostics, physicl electronics, nucler physics, optics nd strophysics 8, 9]. Also, Bker ], Wzwz 32, 33] nd Delves 5] studied the numericl tretment of singulr integrl equtions of the first nd second kinds. In recent yers, mny uthors solved the Abel s integrl equtions of the first nd second kinds by using different pplicble methods 6, 2, 5, 6, 24, 27, 37]. The properties of Abel s integrl equtions cn be found in 20, 25, 28, 33]. One of powerful nd efficient methods to estimte the integrl nd differentil equtions is the colloction method 7, 8, 2]. This method is one of the expnsion methods which is used with different bsis functions. In this pper, the colloction method nd the Tylor polynomils re pplied to estimte the generlized Abel s integrl equtions of the first nd second kinds. Knwl nd Liu in 7] presented the Tylor expnsion pproch for solving integrl equtions. It ws pplied for solving the Volterr- Fredholm integrl equtions 30, 3, 35], system of integrl equtions 22], Volterr integrl equtions 26], integro-differentil equtions 23, 36] nd the others 4, 7]. This pper is orgnized s follows: At first, in Section 2, the Tylor polynomils nd the colloction method re combined to estimte the generlized form of Abel s integrl equtions of the first nd second kinds. In Section 3, the error nlysis of presented method is illustrted nd in sequel, in Section 4, severl pplicble exmples re pproximted. Also, in this section the numericl results nd the bsolute errors for different vlues of x nd n re tbulted. Finlly, Section 5 is conclusion. 2 Min Ide In order to estimte the Abel s integrl equtions of the first nd second kinds, the Tylor polynomils of degree n t x = z is introduced s follows Φ n (x) = n j! Φ(j) (z)(x z) j, x,z b, (3) where the coefficients Φ (j) (z),j = 0,,...,n should be determined. To pproximte the Abel s integrl eqution of the first kind (), we rewrite it in the form n x Φ (j) (z)(t z) ] j j! α dt = g(x), 0 < α < nd x,z b. (4) φ(x) φ(t) In order to estimte the integrl eqution (4), the colloction points ( ) b x i = + i, i = 0,,2,...,n, (5) n 2

3 re substitute into Eq. (4) s follows n xi (t z) j j! φ(xi ) φ(t) Φ (j) (z) = g(x i ), 0 < α <. (6) Now, we cn write Eq. (6) in the form where A = x0 0! x 0! xn 0! (t z) 0 φ(x 0 ) φ(t) (t z) 0 φ(x ) φ(t). (t z) 0 φ(x n) φ(t)!!! x0 x xn Φ (0) (z) Φ () (z) X =. Φ (n) (z) AX = G, (7) (t z) φ(x 0 ) φ(t) (t z) φ(x ) φ(t) n!. (t z) φ(x n) φ(t) (n+) n! n! x0 x xn g(x 0 ) g(x ), G =. g(x n ) (t z) n φ(x 0 ) φ(t) (t z) n φ(x ) φ(t). (t z) n φ(x n) φ(t) (n+) (n+)(n+) By solving the system of equtions (7) nd substitute Φ (j) (z),j = 0,,...,n in Eq. (3) the pproximte solution of Abel s integrl eqution of the first kind () cn be obtined. In order to estimte the second kind form of Abel s integrl eqution (2), by substitute Eq. (3) into Eq. (2) we hve n x n j! Φ(j) (z)(x z) j j! Φ(j) (z)(t z) j (φ(x) φ(t)) α dt = g(x), 0 < α <, (8) nd by substitute colloction points (5) into Eq.(8) the following eqution is obtined n xi (x i z) j j! (t z) j ] (φ(x i ) φ(t)) αdt Φ (j) (c) = g(x i ), 0 < α <. (9)., 3

4 Finlly, we rewrite Eq. (9) in the mtrix form (B A)X = G where mtrices A,X,G were presented nd 0! (x 0 z) 0! (x 0 z) n! (x 0 z) n 0! (x z) 0! (x z) n! (x z) n B =.... 0! (x n z) 0! (x n z) n! (x n z) n (n+)(n+) By solving the obtined system nd determine the unknowns X nd substitute into Eq. (3), the pproximte solution of Eq. (2) cn be obtined. 3 Error nlysis In this section, the error nlysis theorem of this method is presented. Theorem : LetΦ(x) is n exct solution of Eqs. () nd(2), Φ n (x) isn pproximte solution which is obtined from n-th order Tylor-colloction method nd S n (x) = n Φ (i) (z)(x z) i i=0 i! is the Tylor polynomil of degree n t x = z then Φ(x) Φ n (x) M (i+)! Φ (i+) (ξ)+cmx e i (z), ξ,b], where M = mx (x z) i+,c = l, e i (z) = Φ (i) (z) Φ (i) (z). Proof: Assume tht R n (x) is the reminder term of n-th order Tylor polynomil S n t x = z which is given by Thus R n (x) = Φ(x) S n (x) = Φ (i+) (ξ) (i+)! (x z)i+,ξ,b]. R n (x) = Φ(x) S n (x) Φ (i+) (ξ) (i+)!.mx (x z)i+ = M (i+)! Φ (i+) (ξ), ξ,b]. (0) Now, one cn esily write tht Φ(x) Φ n (x) Φ(x) S n (x) + S n (x) Φ n (x) R n (x) + S n (x) Φ n (x). () On the other hnd we hve n ( Φ(i) S n (x) Φ n (x) = (z) Φ (z)) (i) i=0. (x z)i i! 4 E il E i. l C E i, (2)

5 where nd E i = (e 0 (z),e (z),...,e i (z),...,e n (z)), l = (l 0 (z),l (z),...,l i (z),...,l n (z)), e i (z) = Φ (i) (z) Φ (i) (z), l i (z) = (x z)i. i! By substitute Eqs. (0) nd (2) into Eq. (), the bound of bsolute error function cn be obtined s follows Φ(x) Φ n (x) M (i+)! Φ (i+) (ξ)+cmx e i (z), ξ,b]. 4 Numericl illustrtions In this section, severl exmples of generlized Abel s integrl equtions of the first nd second kinds re solved by using presented method 24, 32, 33]. Also, the bsolute error for different vlues of x nd the mximum error re presented in some tbles. The Mthemtic 8 ws pplied to numericl clcultions. Exmple : Consider the following Abel s integrl eqution of the first kind 32, 33] 2 3 πx3 = x 0 (3) x 2 t2dt, where the exct solution is Φ(x) = πx 3. For solving this integrl eqution by using the Tylor polynomils the following lgorithm is presented. Algorithm - Substitute the Tylor polynomil (3) in Eq. (3) for rbitrry vlue n = πx3 = 5 x Φ (j) (z)(t z) j dt. (4) j! 0 x 2 t 2 2- Substitute the colloction point x i from Eq. (5) in Eq. (4) s 2 3 πx3 i = 5 xi Φ (j) (z)(t z) j dt, 0 i n. (5) j! 0 x 2i t2 3- Construct the mtrix form AX = G where A = ,

6 nd G = Determine the unknowns X nd clculte the pproximte solution Φ n (x) Φ 5 (x) = x x x x x 5. The bsolute error for n = 5,7,9 nd different vlues of x is shown in Tble nd the mximum errors re demonstrte in Tble 2. Tble : Numericl results of Exmple for n = 5,7,9. x exct e 5 = Φ(x) Φ 5 (x) e 7 = Φ(x) Φ 7 (x) e 9 = Φ(x) Φ 9 (x) Tble 2: The mximum errors of Exmple. n = 5 n = 7 n = Exmple 2: Consider the generlized Abel s integrl eqution of the first kind 32]: 4 3 sin3 4 (x) = x 0 dt, (6) (sin(x) sin(t)) 4 with the exct solution Φ(x) = cosx. By using the presented method for n = 9 the pproximte solution is obtined s follows Φ 9 (x) = x x x x 4 +. Comprison between exct nd pproximte solutions for different vlues of x nd n = 5,7,9 re presented in Tble 3 nd the mximum errors re demonstrte in Tble 4. Exmple 3: In this exmple, the following first kind Abel s integrl eqution 32] is considered: 6 5 (ex ) 5 6 = x 0 6 (7) dt, (8) (e x e t ) 6

7 Tble 3: Numericl results of Exmple 2 for n = 5,7,9. x exct e 5 = Φ(x) Φ 5 (x) e 7 = Φ(x) Φ 7 (x) e 9 = Φ(x) Φ 9 (x) Tble 4: The mximum errors of Exmple 2. n = 5 n = 7 n = with exct solution Φ(x) = e x. By using presented method, we hve following pproximte solution Φ 9 (x) = x x 2 +. (9) In Tble 5 nd 6, the bsolute errors nd mximum vlues of bsolute error re presented for n = 5,7,9. Tble 5: Numericl results of Exmple 3 for n = 5,7,9. x exct e 5 = Φ(x) Φ 5 (x) e 7 = Φ(x) Φ 7 (x) e 9 = Φ(x) Φ 9 (x) Exmple 4: In this exmple, the second kind Abel s integrl eqution 32, 24] is considered s follows Φ(x) = x x 5 x5 2 dt, 0 (x t) 2 with the exct solution Φ(x) = x 2. Presented method for solving this exmple is exct for n = 5. Exmple 5: Let us consider the following integrl eqution 32]: Φ(x) = 2x 32 2 x x3 4 x 0 dt, (x t) 4 with the exct solution Φ(x) = 2x. In this exmple, the Tylor-colloction method for n = 3 is exct. 7

8 Tble 6: The mximum errors of Exmple 3. n = 5 n = 7 n = Acknowledgements This work hs been supported by Islmic Azd University, Hmedn brnch, reserch grnd. 5 Conclusions In this pper, Tylor polynomils to estimte the generlized form of Abel s integrl equtions of the first nd second kinds were used. Also, bound of bsolute error function ws illustrted. Severl pplicble exmples were solved by using Tylor-colloction method. Presented tbles nd obtined numericl results show the efficiency nd ccurcy of method. References ] CT.H. Bker, The numericl tretment of integrl equtions. Oxford: Clrendon Press; ] H. Brunner, : One hundred yers of Volterr integrl eqution of the first kind, Appl. Numer. Mth. 24 (997) ] C. J. Cremers, R.C. Birkebk, Appliction of the Abel integrl eqution to spectroscopic dt, Appl. Opt. 5 (966) ] S. De, B. N. Mndl, A. Chkrbrti, Use of Abel s integrl equtions in wter wve scttering by two surfce-piercing brriers, Wve Motion 47 (200) ] L.M. Delves, JL. Mohmed, Computtionl methods for integrl equtions. UK: Cmbridge University Press; ] M.A. Friborzi Arghi, H. Dei Ksmei, Numericl Solution of the Second Kind Singulr Volterr Integrl Equtions By Modified Nvot-Simpson s Qudrture, Int. J. Open Problems Compt. Mth. (2008). 7] M.A. Friborzi Arghi, Gh. Kzemi Gelin, Numericl solution of nonliner Hmmerstein integrl equtions vi Sinc-colloction method bsed on double exponentil trnsformtion, Mthemticl Sciences, 7 (203) -7. 8] M.A. Friborzi Arghi, S. Noeighdm, Fiboncci-regulriztion method for solving Cuchy integrl equtions of the first kind, Ain Shms Eng J. 8 (207)

9 9] M.A. Friborzi Arghi, S. Noeighdm, A novel technique bsed on the homotopy nlysis method to solve the first kind Cuchy integrl equtions rising in the theory of irfoils, Journl of Interpoltion nd Approximtion in Scientific Computing 206 () (206) -3. 0] M.A. Friborzi Arghi, S. Noeighdm, Homotopy nlysis trnsform method for solving generlized Abel s fuzzy integrl equtions of the first kind, IEEE (206). DOI: 0.09/CFIS ] M.A. Friborzi Arghi, S. Noeighdm, Homotopy regulriztion method to solve the singulr Volterr integrl equtions of the first kind, Jordn Journl of Mthemtics nd Sttistics (JJMS) () 208, -2. 2] M.A. Friborzi Arghi, S. Yzdni, A Method to Approximte Solution of the First kind Abel Integrl Eqution using Nvot s Qudrture nd Simpson s rule, Interntionl Journl of Industril Mthemtics, (2009) -. 3] R. Gorenflo, S. Vessell, Abel Integrl Equtions: Anlysis nd Applictions, Lecture Notes Mth., 46, Springer, Berlin, 99. 4] P. Hubsomboon, B. Novprteep, H. Kneko, On Tylor-series expnsion methods for the second kind integrl equtions, J. Comput. Appl. Mth. 234 (200) ] L. Hung, Y. Hung, X. F. Li, Approximte solution of Abel integrl eqution, Comput. Mth. Appl. 56 (2008) ] S. Jhnshhi, E. Bbolin, Delfim F.M. Torres, A. Vhidi,Solving Abel integrl equtions of first kind vi frctionl clculus, Journl of King Sud University-Science (27) ] R.P. Knwl, K.C. Liu, A Tylor expnsion pproch for solving integrl equtions, Int. J. Mth. Educ. Sci. Technol. 20 (989) ] O. Knill, R. Dgni, M. Vogel, A new pproch to Abel s integrl opertor nd its ppliction to stellr winds, Astronom. Astrophys. 274 (993) ] L. E. Kosrev, Applictions of integrl equtions of the first kind in experiment physics., Comput. Phys. Commun. 20 (980) ] P.K. Lmm, L. Elden, Numericl solution of first-kind Volterr equtions by sequentil Tikhonov regulriztion, SIAM J. Numer. Anl. 34 (4) (997) ] K. Mleknejd, K. Nedisl, Appliction of Sinc-colloction method for solving clss of nonliner Fredholm integrl equtions, Comput. Mth. Appl. 62 (20) ] K. Mleknejd, N. Aghzdeh, M. Rbbni, Numericl omputtionl solution of second kind Fredholm integrl equtions system by using Tylor-series expnsion method, Appl. Mth. Comput. 75 (2006)

10 23] K. Mleknejd, Y. Mhmoudi, Tylor polynomil solution of high-order nonliner Volterr- Fredholm integro-differentil equtions, Appl. Mth. Comput. 45 (2003) ] S. Noeighdm, E. Zrei, H. Brzegr Kelishmi, Homotopy nlysis trnsform method for solving Abel s integrl equtions of the first kind, Ain Shms Eng. J. 7 (206) ] A. D. Polynin, A.V. Mnzhirov, Hndbook of Integrl Equtions, CRC Press, Boc Rton, ] M. Sezer, Tylor polynomil solution of Volterr integrl equtions, Int. J. Mth. Educ. Sci. Technol. 25 (994) ] S. Sohrbi, Comprison Chebyshev wvelets method with BPFs method for solving Abel s integrl eqution, Ain Shms Engineering Journl (2) ] H.J. Teriele, Colloction method for wekly singulr second kind Volterr integrl equtions with non-smooth solution, IMA J. Numer. Anl. 2 (4) (982) ] J. Wlton, Systems of generlized Abel integrl equtions with pplictions to simultneous dul reltions, SIAM J. Mth. Anl. 0 (979) ] K. Wng, Q. Wng, Tylor polynomil method nd error estimtion for kind of mixed Volterr-Fredholm integrl equtions, Appl. Mth. Comput. 229 (204) ] K. Wng, Q. Wng, Tylor colloction method nd convergence nlysis for the Volterr- Fredholm integrl equtions, J. Comput. Appl. Mth. 260 (204) ] A.M. Wzwz, Liner nd nonliner integrl equtions: methods nd pplictions. Berlin: Higher Eduction, Beijing, nd Springer; ] A.M. Wzwz, A First Course in Integrl Equtions, World Scientific Publishing, Singpore, ] A.M. Wzwz, M.S. Mehnn, The combined Lplce-Adomin method for hndling singulr integrl eqution of het trnsfer, Int. J. Nonliner Sci. 0 (200) ] S. Ylcinbs, Tylor polynomil solutions of nonliner Volterr-Fredholm integrl equtions, Appl. Mth. Comput. 27 (2002) ] S. Ylcinbs, M. Sezer, The pproximte solution of high-order liner Volterr-Fredholm integro-differentil equtions in terms of Tylor polynomils, Appl. Mth. Comput. 2 (2000) ] C. Yng, An efficient numericl method for solving Abel integrl eqution, Appl. Mth. Comput. 227 (204) ] N. Zeilon, Sur quelques points de l theorie de l equtionintegrled Abel, Arkiv. Mt. Astr. Fysik. 8 (924) -9. 0