Analytic formulas for the evaluation of the Pearcey integral
|
|
- Elwin Simon
- 6 years ago
- Views:
Transcription
1 Aalytic formulas for the evaluatio of the Pearcey itegral José L. López 1 ad Pedro J. Pagola 1 arxiv: v1 [math.na] 1 Ja Dpto. de Igeiería Matemática e Iformática, Uiversidad Pública de Navarra jl.lopez@uavarra.es, pedro.pagola@uavarra.es Abstract We ca fid i the literature several coverget ad/or asymptotic expasios of the Pearcey itegral P x, y i differet regios of the complex variables x ad y, but they do ot cover the whole complex x ad y plaes. The purpose of this paper is to complete this aalysis givig ew coverget ad/or asymptotic expasios that, together the kow oes, let the evaluatio of the Pearcey itegral i a large regio of the complex x ad y plaes. The accuracy of the approximatios derived i this paper is illustrated some umerical experimets. Moreover, the expasios derived here are simpler compared other kow expasios, as they are derived from a simple maipulatio of the itegral defiitio of P x, y. 010 AMS Mathematics Subject Classificatio: 33E0; 1A60. Keywords & Phrases: Pearcey itegral. Coverget ad asymptotic expasios. Watso lemma. 1 Itroductio The Pearcey itegral is the secod caoical diffractio itegrals [1, Sec. 36.1]. Apart from their mathematical importace i the uiform asymptotic approximatio of oscillatory itegrals [8], the caoical diffractio itegrals have physical applicatios i the descriptio of surface gravity waves [7], [13], bifurcatio sets, optics, quatum mechaics ad acoustics see [1, Sec. 36.1] ad refereces there i. The Cusp catastrophe or Pearcey itegral [1, p.777, eq ] is the itegral: P x, y := e it +xt +yt dt. 1 This itegral was first evaluated umerically by usig quadrature formulas i [1] i the cotext of the ivestigatio of the electromagetic field ear a cusp. I [1, Chap. 36] 1
2 we ca fid may properties such as symmetries, ilustrarive pictures, bifurcatio sets, scalig relatios, zeros, coverget series expasios, differetial equatios ad leadigorder asymptotic approximatios amog others. But we caot fid may details about aalytic approximatio formulas asymptotic expasios i particular or umerical evaluatio techiques. The itegral 1 exists oly for 0 arg x π ad real y. As it is idicated i [10], after a rotatio of the itegratio path through a agle of π/8 that removes the rapidly oscillatory term e it, the Pearcey itegral may be writte i the form P x, y = e iπ/8 P xe iπ/, ye iπ/8, P x, y := 0 e t xt cosytdt = 1 e t +iyt e xt dt. This itegral is absolutely coverget for all complex values of x ad y ad represets the aalytic cotiuatio of the Pearcey itegral P x, y to all complex values of x ad y [10]. Therefore, it is more coveiet to work the represetatio of the Pearcey itegral. We summarize below the more elemetary aalytic expasios coverget or asymptotic of P x, y that we ca fid i the literature ad give refereces to other more elaborated expasios. A coverget series expasio of P x, y may be foud i [1, p. 787, eqs ad 36.8.]: P x, y = Γ a x, y, 3 a 0 x, y = 1, a 1 x, y = y ad, for =, 3,,..., a x, y = 1 [y a 1x, y + x a x, y]. Aother coverget expasio of P x, y may be obtaied after a expasio of the cosie term i ad iterchage of sum ad itegral [11]: P x := P x, y = 1 +3/ Γ + 1 y! P x, x, y C, U + 1, 1; x 1 Γ + 1 M + 1, 1; x xγ + 3 M + 3, 3; x if Rx 0, if Rx < 0, 5 where Ma, b; z ad Ua, b; z are cofluet hypergeometric fuctios [9, Chap. 13] Apart from coverget expasios, we ca also fid i the literature several asymptotic expasios of P x, y. I [] we ca fid a asymptotic expasio of the Pearcey itegral whe x, y are ear the caustic 8x 3 7y = 0 that remais valid as x. The expasio is give i terms of Airy fuctios ad its derivatives ad the coefficiets are computed recursively. We refer the reader to [] for further details.
3 A exhaustive asymptotic aalysis of this itegral ca be foud i [10]. I particular, a complete asymptotic expasio is give i [10] by usig asymptotic techiques for itegrals applied to the itegral. The asymptotic aalysis of this itegral for large x is divided i two regios: arg x < π ad arg x > π. I the first regio we fid that [10] P x, y 1 x e y /x 1 Γ + 1/ M ; 1/; y. 6!x x The asymptotic expasio i the secod regio is a little bit more cumbersome, we refer to [10] for details. Aother coverget ad asymptotic expasio of P x, y for large x, derived from a differetial equatio satisfied by P x, y, may be foud i [5]. A complete asymptotic expasio for large y is derived i [6], where the asymptotic aalysis is divided i three differet regios: arg y π, π arg y < π ad π < arg y π. We give details here for the last regio ad refer the reader to [6] for further details. Whe π < arg y π we have that: 8 where P x, y π 1 e3 /3e iπ/3y/3 /3xy/3e iπ/3 +x /6+iπ/6 3 5/6 y 1/3 A x := m= +1/ b,m,k x := m b,m,k xc +m k x, xk m k/3 1 m k!m k! m! A x, 7 y/3 ad c + x = x 3 c +1x /3 3 c x, = 0, 1,,... /3 c 0 x = 1 ad c 1 x = x 3. 1/3 I [11] we ca fid the hyperasymptotic evaluatio of the Pearcey itegral for real values of x ad y usig Hadamard expasios. The Pearcey itegral is writte i terms of a ifiite series whose terms are Hadamard series. The terms of these Hadamard series are icomplete gamma fuctios. The aalytic expressio is sophisticated ad we refer the reader to [11] for details. I this paper we exploit the followig simple idea: the itegrad i is a expoetial. The, we maipulate the fuctio i the expoet i order to factorize the expoetial ad replace oe of the factors by its Taylor expasio at the origi. We cosider three differet possibilities. I the followig sectio we derive a complete coverget ad asymptotic expasio of P x, y for small x. I Sectio 3 we derive a complete asymptotic expasio of P x, y for large x. I Sectio we derive a complete asymptotic expasio of P x, y for large x ad y. Sectio 5 cotais some umerical experimets ad a few remarks. Through the paper we use the pricipal argumet arg z π, π] for ay complex umber z. 3
4 A coverget ad asymptotic expasio for small x Because P x, y = P x, y, out loss of geerality, we may restrict ourselves to Ry 0. A coverget expasio of P x, y that is also ad asymptotic expasio whe x 0 ca be obtaied by cosiderig the Taylor expasio of the expoetial at the origi: e xt = 1 x k t k k! + r xt, 8 where r t is the Taylor remaider. Replacig this expasio i the itegrad of the right had side of formula ad iterchagig sum ad itegral we fid: P k y := 1 P x, y = 1 x k P k y + R x, y, 9 k! t k e t +iyt dt, R x, y := 1 r xt e t +iyt dt. After some straightforward maipulatios, we obtai a explicit formula for the coefficiets P k y i 9 k = 0, 1,,...: P k y := k 1 + k Γ 1F 3 ; 1,, 3 y ; y 3 + k 3 + k 8 Γ 1F 3 ; 3, 5, 6 y ;. Itegratig by parts we derive the followig five terms recurrece for the coefficiets P k y: y k + 1k + 1 P k+ y = Pk+1 y P k y+k+1k+7/p k+ y, k = 0, 1,,... 8 P 0 y = 1 1 Γ 0F ; 1, 3 ; y y Γ 0F ; 5, 3 ; y, 56 P 1 y = 1 3 Γ 0F ; 1, 1 ; y y Γ 0F ; 3, 3 ; y, 56 P y = Γ 1F 3 ; 1, 1, 3 ; y y Γ 1F 3 ; 3, 5, 3 ; y 56 P 3 y = Γ 1F 3 ; 1, 1, 3 ; y y Γ 0F ; 3, 5, 3 ; y 56 Therefore, we have obtaied that, for ay complex x ad y, P x, y = 1,. P k y x k + R x, y. 10 k!
5 This expasio is coverget for ay value of x ad y ad is a asymptotic expasio of P x, y as x 0. Moreover, usig the Lagrage formula for the remaider i 8 we have r xt = r xt < e ξ Rx! xt 0 < ξ < t t! x if Rx 0, t e Rx t x if Rx < 0,! I the first case Rx 0, we immediately coclude that R x, y x! P i Iy. I secod case Rx < 0, we have R x, y = 1 r xt e t +iyt dt 1 x! t e Rx t e t Iyt dt. Writig the domiat term e t of itegrad i the form e t e t we ca rewrite the above itegral as 1 x! gt e t Iyt t dt gt := e t + Rx t. The fuctio gt attais it absolute maximum at t = Rx ad the gt g Rx = e R x. Therefore, we coclude that, for Rx < 0, R x, y x! R e x P i Iy. 3 A asymptotic expasio for large x Usig the Cauchy s residue theorem we fid that we ca replace the itegratio path, i by ay straight Γ := {u = se iσ : < s < }, σ < π/8: Defie the parameter α := u x 1/ e iθ/ we fid P x, y = 1 e iσ e iσ e u xu +iyu du. y x ad θ := arg x. After the chage of variable t = u x = P x, y = 1 e iσ+θ/ e t +iαt t /x e dt. 11 x e iσ+θ/ 5
6 From here, a asymptotic expasio of the itegral 11 for large x ca be obtaied cosiderig the Taylor expasio of the expoetial at the origi: e t /x = 1 1 k t k k!x k + r t x, 1 where r t is the Taylor remaider. From the Lagrage s formula for the remaider we have: t r = r x s e i σ = t!x e ξ ei σ 0 < ξ < s. The r t x < t. 13! x Replacig expasio 1 i the above itegral ad iterchagig sum ad itegral we fid y π e 1 x P x, y = 1 k x k! x Q kα + R k x, y, Q k α := e y x e iσ+θ/ π e iσ+θ/ t k e t +iαt dt, R x, y := e y x e iσ+θ/ t r e t +iαt dt. π e x iσ+θ/ 1 These itegrals exist whe arg x + σ < π/. For ay value of x such that arg x < 3π/, we ca chose σ σ < π/8 such that these itegrals are fiite. After some maipulatios ad usig the itegral represetatio of Hermite polyomials [eq , [3]] we fid H x = i e x π Q k α = k H k α e t t e ixt dt, k = 0, 1,,... Moreover, itroducig 13 i 1 we fid R x, y cos arg y arg x σ cosσ+arg x y x e! x cos +1/ σ + arg x H y siσ + arg y i. x cosσ + arg x It is clear that R x, y = Ox as x uiformly i y bouded y/ x ad arg x < 3π/. Therefore, we have obtaied that P x, y = π e y x x 1 1 k y k!x H k k + R x, y x, 15 6
7 is a asymptotic expasio of P x, y whe x uiformly i y bouded y/ x ad arg x < 3π/. This expasio is just the expasio 3. derived i [10]. The derivatio of 3. is much more cumbersome as it is obtaied from a cotour itegral represetatio of the Pearcey itegral i terms of a parabolic cylider fuctio. A error boud for the remaider is ot provided there. O the other had, a asymptotic expasio for arg x < π/ is give i [10]. A asymptotic expasio for large x ad y Defie the parameter γ := y. We ca write the Pearcey itegral i the form x P x, y = 1 e xft t dt, 16 phase fuctio ft := iγt t. The uique saddle poit of this phase fuctio is the poit t 0 = iγ. From [] we kow that the asymptotic aalysis of this itegral does ot require the applicatio of the stadard saddle poit method, but just a chage of variable. Usig the Cauchy s residue theorem we fid that we ca replace the itegratio path, i 16 by ay straight Γ := {iγ + t; t = re iγ : < r < }, σ < π/8: x P x, y = e γ e xt e iγ+t dt. 17 We write the above itegral i the form P x, y = e α x α e iσ x 6α t e e iσ htdt, 18 ht := e iγ3 t iγt 3 t. From here, a asymptotic expasio of the itegral 17 for large x ca be derived from Watso s lemma. From the Taylor expasio of the expoetial at the origi we have that ht = A γt, A γ := / k/3 j=0 iγ 3 k 3j iγ j 1 k. k!j! k 3j! Replacig this expasio i the itegral 18 ad iterchagig sum ad itegral we fid that /x y /16x P x, y e y A γp γ, x, P γ, x := t e x 6γ t dt, t = re iγ, < r <. 7
8 This itegral exists whe argx 6t = re iγ : < r < + σ < π/. For ay value of x ad y such that argx 3y /x < 3π/, we ca chose σ σ < π/8 such that this itegral is fiite: Γ + 1/ P γ, x = x 6γ. +1/ Therefore, we have obtaied that, whe x, /x y /16x P x, y e y x x 3 3y y x A Γ + 1/, 19 x x 3 3y uiformly i y bouded y/x ad argx 3y /x < 3π/. 5 Fial remarks ad umerical experimets Figure 1 summarizes the regios for x, y covered by the above metioed algorithms except [11] that is more sophisticated ad valid oly for real x, y ad by the three algorithms derived below i this paper. y x Small Moderate Large Small 3 Sec. [] 6 Sec. 3 Moderate Sec. 1 [] 6 Sec. 3 Large 7 7 Sec. Figure 1:I blue, the regios of x, y for which there are aalytic formulas for the computatio of P x, y available i the literature ad detailed i the itroductio sectio. The regios for which the ew formulas derived here are useful are dashed. The curve i the large x ad y box idicates the caustic. The expasio give i Sectio is ew ad valid for ay complex values of x ad y. The expasio give i Sectio 3 is just the oe give i eq. 3. of [10], although derived i a more straightforward way that permits a accurate error boud for the remaider. It is valid for arg x < 3π/ ad ay complex y. The expasio give i Sectio is also ew ad permits a wider regio for x, y tha the expasio give i []. The expasio 8
9 give i [] is restricted to the viciity of the caustic 8x 3 7y = 0, whereas the oe derived i Sectio above is valid for argx 3y /x < 3π/. Fially, we illustrate the accuracy of the above expasios. I the followig tables we show the relative error obtaied these algorithms for several values of x, y ad differet orders of the approximatio. As we do ot have at our disposal the exact value of the Pearcey itegral, we have take the umerical itegratio obtaied the program Mathematica workig precisio= 00 as the exact value of P x, y. x, y , i e iπ/, /5, i /10 i/8, i/0, i Table 1: Relative error for several values of x, y ad the umber of terms of the approximatio give by 10. x, y , i e iπ/, i, 1 i , i, i i, i Table : Relative error for several values of x, y ad the umber of terms of the approximatio give by 15. x, y , 3i e 5iπ/8, i, 0e iπ/ , 0i i, Table 3: Relative error for several values of x, y ad the umber of terms of the approximatio give by 19. 9
10 6 Ackowledgmets This research was supported by the Spaish Miistry of Ecoomía y Competitividad, project MTM The Uiversidad Pública de Navarra is ackowledged by its fiacial support. Refereces [1] M. V. Berry ad C. J. Howls, Itegrals coalescig saddles, i: NIST Hadbook of Mathematical Fuctios, Cambridge Uiversity Press, Cambridge, 010, pp Chapter 36. [] D. Kamiski, Asymptotic expasio of the Pearcey itegral ear the caustic, J. SIAM J. Math. Aal., , [3] T. H. Koorwider, R. Koekoek, R. F. Swarttouw, R. Wog, Orthogoal Polyomials, i: NIST Hadbook of Mathematical Fuctios, Cambridge Uiv. Press, 010, pp Chapter 18. [] J. L. López ad N. Temme, Uiform approximatios of Beroulli ad Euler polyomials i terms of hyperbolic fuctios, Stud. Appl. Math., , [5] J. L. López ad P. Pagola, Coverget ad asymptotic expasios of the Pearcey itegral. J. Math. Aal. Appl., 30, o. 1, [6] J. L. López ad P. Pagola, The Pearcey itegral i the highly oscillatory regio. Submitted. [7] Lord Kelvi, Deep water ship-waves, Phil. Mag., , [8] A. B. Olde Daalhuis, O the asymptotics for late coefficiets i uiform asymptotic expasios of itegrals coalescig saddles Meth. Appl. Aal., 7 o. 000, [9] A. B. Olde Daalhuis, Cofluet Hypergeometric Fuctios, i: NIST Hadbook of Mathematical Fuctios, Cambridge Uiversity Press, Cambridge, 010, pp Chapter 13. [10] R. B. Paris, The asymptotic behaviour of Pearcey s itegral for complex variables, Proc. Roy. Soc. Lodo Ser. A., 3 o , [11] R. B. Paris ad D. Kamiski, Hyperasymptotic evaluatio of the Pearcey itegral via Hadamard expasios, J. Comput. Appl. Math., , [1] T. Pearcey, The structure of a electromagetic field i the eighbourhood of a cusp of a caustic, Phi. Mag., , [13] F. Ursell, Itegrals a large parameter: several early coicidet saddle poits, Proc. Camb. Phil. Soc., 7 197, [1] R. Wog, Asymptotic approximatios of itegrals, Academic Press, New York,
The Pearcey integral in the highly oscillatory region
The Pearcey integral in the highly oscillatory region José L. López 1 and Pedro Pagola arxiv:1511.0533v1 [math.ca] 17 Nov 015 1 Dpto. de Ingeniería Matemática e Informática and INAMAT, Universidad Pública
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationThe Riemann Zeta Function
Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we
More informationC. Complex Numbers. x 6x + 2 = 0. This equation was known to have three real roots, given by simple combinations of the expressions
C. Complex Numbers. Complex arithmetic. Most people thik that complex umbers arose from attempts to solve quadratic equatios, but actually it was i coectio with cubic equatios they first appeared. Everyoe
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationThe Gamma function. Marco Bonvini. October 9, dt e t t z 1. (1) Γ(z + 1) = z Γ(z) : (2) = e t t z. + z dt e t t z 1. = z Γ(z).
The Gamma fuctio Marco Bovii October 9, 2 Gamma fuctio The Euler Gamma fuctio is defied as Γ() It is easy to show that Γ() satisfy the recursio relatio ideed, itegratig by parts, dt e t t. () Γ( + ) Γ()
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationTopics in Probability Theory and Stochastic Processes Steven R. Dunbar. Stirling s Formula Derived from the Gamma Function
Steve R. Dubar Departmet of Mathematics 23 Avery Hall Uiversity of Nebraska-Licol Licol, NE 68588-3 http://www.math.ul.edu Voice: 42-472-373 Fax: 42-472-8466 Topics i Probability Theory ad Stochastic Processes
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationFourier Series and their Applications
Fourier Series ad their Applicatios The fuctios, cos x, si x, cos x, si x, are orthogoal over (, ). m cos mx cos xdx = m = m = = cos mx si xdx = for all m, { m si mx si xdx = m = I fact the fuctios satisfy
More informationHarmonic Number Identities Via Euler s Transform
1 2 3 47 6 23 11 Joural of Iteger Sequeces, Vol. 12 2009), Article 09.6.1 Harmoic Number Idetities Via Euler s Trasform Khristo N. Boyadzhiev Departmet of Mathematics Ohio Norther Uiversity Ada, Ohio 45810
More informationj=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)
Problem. Compute the itegrals C r d for Z, where C r = ad r >. Recall that C r has the couter-clockwise orietatio. Solutio: We will use the idue Theorem to solve this oe. We could istead use other (perhaps
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationEvaluation of Some Non-trivial Integrals from Finite Products and Sums
Turkish Joural of Aalysis umber Theory 6 Vol. o. 6 7-76 Available olie at http://pubs.sciepub.com/tjat//6/5 Sciece Educatio Publishig DOI:.69/tjat--6-5 Evaluatio of Some o-trivial Itegrals from Fiite Products
More informationTR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT
TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the
More informationOn a class of convergent sequences defined by integrals 1
Geeral Mathematics Vol. 4, No. 2 (26, 43 54 O a class of coverget sequeces defied by itegrals Dori Adrica ad Mihai Piticari Abstract The mai result shows that if g : [, ] R is a cotiuous fuctio such that
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationCOMPUTING THE EULER S CONSTANT: A HISTORICAL OVERVIEW OF ALGORITHMS AND RESULTS
COMPUTING THE EULER S CONSTANT: A HISTORICAL OVERVIEW OF ALGORITHMS AND RESULTS GONÇALO MORAIS Abstract. We preted to give a broad overview of the algorithms used to compute the Euler s costat. Four type
More informationMonte Carlo Integration
Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce
More informationRotationally invariant integrals of arbitrary dimensions
September 1, 14 Rotatioally ivariat itegrals of arbitrary dimesios James D. Wells Physics Departmet, Uiversity of Michiga, A Arbor Abstract: I this ote itegrals over spherical volumes with rotatioally
More informationAn Asymptotic Expansion for the Number of Permutations with a Certain Number of Inversions
A Asymptotic Expasio for the Number of Permutatios with a Certai Number of Iversios Lae Clark Departmet of Mathematics Souther Illiois Uiversity Carbodale Carbodale, IL 691-448 USA lclark@math.siu.edu
More informationThe Gamma function Michael Taylor. Abstract. This material is excerpted from 18 and Appendix J of [T].
The Gamma fuctio Michael Taylor Abstract. This material is excerpted from 8 ad Appedix J of [T]. The Gamma fuctio has bee previewed i 5.7 5.8, arisig i the computatio of a atural Laplace trasform: 8. ft
More informationAnalytic Continuation
Aalytic Cotiuatio The stadard example of this is give by Example Let h (z) = 1 + z + z 2 + z 3 +... kow to coverge oly for z < 1. I fact h (z) = 1/ (1 z) for such z. Yet H (z) = 1/ (1 z) is defied for
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationSolutions to Homework 1
Solutios to Homework MATH 36. Describe geometrically the sets of poits z i the complex plae defied by the followig relatios /z = z () Re(az + b) >, where a, b (2) Im(z) = c, with c (3) () = = z z = z 2.
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationAMS Mathematics Subject Classification : 40A05, 40A99, 42A10. Key words and phrases : Harmonic series, Fourier series. 1.
J. Appl. Math. & Computig Vol. x 00y), No. z, pp. A RECURSION FOR ALERNAING HARMONIC SERIES ÁRPÁD BÉNYI Abstract. We preset a coveiet recursive formula for the sums of alteratig harmoic series of odd order.
More informationDirichlet s Theorem on Arithmetic Progressions
Dirichlet s Theorem o Arithmetic Progressios Athoy Várilly Harvard Uiversity, Cambridge, MA 0238 Itroductio Dirichlet s theorem o arithmetic progressios is a gem of umber theory. A great part of its beauty
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationREGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS
REGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS LIVIU I. NICOLAESCU ABSTRACT. We ivestigate the geeralized covergece ad sums of series of the form P at P (x, where P R[x], a R,, ad T : R[x] R[x]
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationk-generalized FIBONACCI NUMBERS CLOSE TO THE FORM 2 a + 3 b + 5 c 1. Introduction
Acta Math. Uiv. Comeiaae Vol. LXXXVI, 2 (2017), pp. 279 286 279 k-generalized FIBONACCI NUMBERS CLOSE TO THE FORM 2 a + 3 b + 5 c N. IRMAK ad M. ALP Abstract. The k-geeralized Fiboacci sequece { F (k)
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationSequences of Definite Integrals, Factorials and Double Factorials
47 6 Joural of Iteger Sequeces, Vol. 8 (5), Article 5.4.6 Sequeces of Defiite Itegrals, Factorials ad Double Factorials Thierry Daa-Picard Departmet of Applied Mathematics Jerusalem College of Techology
More informationA q 2 -Analogue Operator for q 2 -Analogue Fourier Analysis
JOURNAL OF MATEMATICAL ANALYSIS AND APPLICATIONS, 5758 997 ARTICLE NO AY975547 A -Aalogue Operator for -Aalogue Fourier Aalysis Richard L Rubi Departmet of Mathematics, Florida Iteratioal Uiersity, Miami,
More informationTHE N-POINT FUNCTIONS FOR INTERSECTION NUMBERS ON MODULI SPACES OF CURVES
THE N-POINT FUNTIONS FOR INTERSETION NUMBERS ON MODULI SPAES OF URVES KEFENG LIU AND HAO XU Abstract. We derive from Witte s KdV equatio a simple formula of the -poit fuctios for itersectio umbers o moduli
More informationBenaissa Bernoussi Université Abdelmalek Essaadi, ENSAT de Tanger, B.P. 416, Tanger, Morocco
EXTENDING THE BERNOULLI-EULER METHOD FOR FINDING ZEROS OF HOLOMORPHIC FUNCTIONS Beaissa Beroussi Uiversité Abdelmalek Essaadi, ENSAT de Tager, B.P. 416, Tager, Morocco e-mail: Beaissa@fstt.ac.ma Mustapha
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics ad Egieerig Lecture otes Sergei V. Shabaov Departmet of Mathematics, Uiversity of Florida, Gaiesville, FL 326 USA CHAPTER The theory of covergece. Numerical sequeces..
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More informationarxiv: v1 [math.co] 3 Feb 2013
Cotiued Fractios of Quadratic Numbers L ubomíra Balková Araka Hrušková arxiv:0.05v [math.co] Feb 0 February 5 0 Abstract I this paper we will first summarize kow results cocerig cotiued fractios. The we
More informationA PROOF OF THE TWIN PRIME CONJECTURE AND OTHER POSSIBLE APPLICATIONS
A PROOF OF THE TWI PRIME COJECTURE AD OTHER POSSIBLE APPLICATIOS by PAUL S. BRUCKMA 38 Frot Street, #3 aaimo, BC V9R B8 (Caada) e-mail : pbruckma@hotmail.com ABSTRACT : A elemetary proof of the Twi Prime
More informationDe la Vallée Poussin Summability, the Combinatorial Sum 2n 1
J o u r a l of Mathematics ad Applicatios JMA No 40, pp 5-20 (2017 De la Vallée Poussi Summability, the Combiatorial Sum 1 ( 2 ad the de la Vallée Poussi Meas Expasio Ziad S. Ali Abstract: I this paper
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationIntroduction of an Elementary Method to Express ζ(2n+1) in Terms of ζ(2k) with k 1
Itroductio of a Elemetary Method to Express ζ+ i Terms of ζk with k arxiv:85.3v [math-ph] 3 May 8 Kazuyuki FUJII ad Tatsuo SUZUKI Departmet of Mathematical Scieces Yokohama City Uiversity Yokohama, 36
More informationON THE LEHMER CONSTANT OF FINITE CYCLIC GROUPS
ON THE LEHMER CONSTANT OF FINITE CYCLIC GROUPS NORBERT KAIBLINGER Abstract. Results of Lid o Lehmer s problem iclude the value of the Lehmer costat of the fiite cyclic group Z/Z, for 5 ad all odd. By complemetary
More informationarxiv: v1 [math.nt] 5 Jan 2017 IBRAHIM M. ALABDULMOHSIN
FRACTIONAL PARTS AND THEIR RELATIONS TO THE VALUES OF THE RIEMANN ZETA FUNCTION arxiv:70.04883v [math.nt 5 Ja 07 IBRAHIM M. ALABDULMOHSIN Kig Abdullah Uiversity of Sciece ad Techology (KAUST, Computer,
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More informationThe Bilateral Laplace Transform of the Positive Even Functions and a Proof of Riemann Hypothesis
The Bilateral Laplace Trasform of the Positive Eve Fuctios ad a Proof of Riema Hypothesis Seog Wo Cha Ph.D. swcha@dgu.edu Abstract We show that some iterestig properties of the bilateral Laplace trasform
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationTERMWISE DERIVATIVES OF COMPLEX FUNCTIONS
TERMWISE DERIVATIVES OF COMPLEX FUNCTIONS This writeup proves a result that has as oe cosequece that ay complex power series ca be differetiated term-by-term withi its disk of covergece The result has
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationBernoulli numbers and the Euler-Maclaurin summation formula
Physics 6A Witer 006 Beroulli umbers ad the Euler-Maclauri summatio formula I this ote, I shall motivate the origi of the Euler-Maclauri summatio formula. I will also explai why the coefficiets o the right
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationIIT JAM Mathematical Statistics (MS) 2006 SECTION A
IIT JAM Mathematical Statistics (MS) 6 SECTION A. If a > for ad lim a / L >, the which of the followig series is ot coverget? (a) (b) (c) (d) (d) = = a = a = a a + / a lim a a / + = lim a / a / + = lim
More informationFall 2013 MTH431/531 Real analysis Section Notes
Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters
More informationNotes 19 Bessel Functions
ECE 638 Fall 17 David R. Jackso Notes 19 Bessel Fuctios Notes are from D. R. Wilto, Dept. of ECE 1 Cylidrical Wave Fuctios Helmholtz equatio: ψ + k ψ = I cylidrical coordiates: ψ 1 ψ 1 ψ ψ ρ ρ ρ ρ φ z
More informationCarleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.
Carleto College, Witer 207 Math 2, Practice Fial Prof. Joes Note: the exam will have a sectio of true-false questios, like the oe below.. True or False. Briefly explai your aswer. A icorrectly justified
More informationA Negative Result. We consider the resolvent problem for the scalar Oseen equation
O Osee Resolvet Estimates: A Negative Result Paul Deurig Werer Varhor 2 Uiversité Lille 2 Uiversität Kassel Laboratoire de Mathématiques BP 699, 62228 Calais cédex Frace paul.deurig@lmpa.uiv-littoral.fr
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationarxiv: v1 [math.ca] 29 Jun 2018
URAL MATHEMATICAL JOURNAL, Vol. 3, No., 207 arxiv:807.025v [math.ca] 29 Ju 208 EVALUATION OF SOME NON-ELEMENTARY INTEGRALS INVOLVING SINE, COSINE, EXPONENTIAL AND LOGARITHMIC INTEGRALS: PART II Victor
More informationShivley s Polynomials of Two Variables
It. Joural of Math. Aalysis, Vol. 6, 01, o. 36, 1757-176 Shivley s Polyomials of Two Variables R. K. Jaa, I. A. Salehbhai ad A. K. Shukla Departmet of Mathematics Sardar Vallabhbhai Natioal Istitute of
More informationThe r-generalized Fibonacci Numbers and Polynomial Coefficients
It. J. Cotemp. Math. Scieces, Vol. 3, 2008, o. 24, 1157-1163 The r-geeralized Fiboacci Numbers ad Polyomial Coefficiets Matthias Schork Camillo-Sitte-Weg 25 60488 Frakfurt, Germay mschork@member.ams.org,
More informationSome remarks for codes and lattices over imaginary quadratic
Some remarks for codes ad lattices over imagiary quadratic fields Toy Shaska Oaklad Uiversity, Rochester, MI, USA. Caleb Shor Wester New Eglad Uiversity, Sprigfield, MA, USA. shaska@oaklad.edu Abstract
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differetial Equatios CAAM 45 Sprig 005 Lecture 4 -step time-steppig methods: stability, accuracy Ruge-Kutta Methods, Istructor: Tim Warburto Today Recall AB stability regios
More informationMATH 1080: Calculus of One Variable II Fall 2017 Textbook: Single Variable Calculus: Early Transcendentals, 7e, by James Stewart.
MATH 1080: Calculus of Oe Variable II Fall 2017 Textbook: Sigle Variable Calculus: Early Trascedetals, 7e, by James Stewart Uit 3 Skill Set Importat: Studets should expect test questios that require a
More informationGENERALIZED HARMONIC NUMBER IDENTITIES AND A RELATED MATRIX REPRESENTATION
J Korea Math Soc 44 (2007), No 2, pp 487 498 GENERALIZED HARMONIC NUMBER IDENTITIES AND A RELATED MATRIX REPRESENTATION Gi-Sag Cheo ad Moawwad E A El-Miawy Reprited from the Joural of the Korea Mathematical
More informationOrthogonal polynomials derived from the tridiagonal representation approach
Orthogoal polyomials derived from the tridiagoal represetatio approach A. D. Alhaidari Saudi Ceter for Theoretical Physics, P.O. Box 374, Jeddah 438, Saudi Arabia Abstract: The tridiagoal represetatio
More informationA Note on Integrals & Hybrid Contours in the Complex Plane
A Note o Itegrals & Hybrid Cotours i the Complex Plae Joh Gill July 4 Abstract: Cotour itegrals ca be expressed graphically as simple vectors arisig from a secodary cotour. We start with a well-behaved
More informationON POINTWISE BINOMIAL APPROXIMATION
Iteratioal Joural of Pure ad Applied Mathematics Volume 71 No. 1 2011, 57-66 ON POINTWISE BINOMIAL APPROXIMATION BY w-functions K. Teerapabolar 1, P. Wogkasem 2 Departmet of Mathematics Faculty of Sciece
More informationMa 4121: Introduction to Lebesgue Integration Solutions to Homework Assignment 5
Ma 42: Itroductio to Lebesgue Itegratio Solutios to Homework Assigmet 5 Prof. Wickerhauser Due Thursday, April th, 23 Please retur your solutios to the istructor by the ed of class o the due date. You
More informationOn the Lebesgue constant for the Xu interpolation formula
O the Lebesgue costat for the Xu iterpolatio formula Le Bos Dept. of Mathematics ad Statistics, Uiversity of Calgary Caada Stefao De Marchi Dept. of Computer Sciece, Uiversity of Veroa Italy Marco Viaello
More informationThe log-behavior of n p(n) and n p(n)/n
Ramauja J. 44 017, 81-99 The log-behavior of p ad p/ William Y.C. Che 1 ad Ke Y. Zheg 1 Ceter for Applied Mathematics Tiaji Uiversity Tiaji 0007, P. R. Chia Ceter for Combiatorics, LPMC Nakai Uivercity
More informationThe natural exponential function
The atural expoetial fuctio Attila Máté Brookly College of the City Uiversity of New York December, 205 Cotets The atural expoetial fuctio for real x. Beroulli s iequality.....................................2
More information(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationSome p-adic congruences for p q -Catalan numbers
Some p-adic cogrueces for p q -Catala umbers Floria Luca Istituto de Matemáticas Uiversidad Nacioal Autóoma de México C.P. 58089, Morelia, Michoacá, México fluca@matmor.uam.mx Paul Thomas Youg Departmet
More informationInfinite Products Associated with Counting Blocks in Binary Strings
Ifiite Products Associated with Coutig Blocks i Biary Strigs J.-P. Allouche U. A. 226, Mathématiques et Iformatique J. O. Shallit Departmet of Computer Sciece Uiversité de Bordeaux Uiversity of Chicago
More informationSHARP INEQUALITIES INVOLVING THE CONSTANT e AND THE SEQUENCE (1 + 1/n) n
SHARP INEQUALITIES INVOLVING THE CONSTANT e AND THE SEQUENCE + / NECDET BATIR Abstract. Several ew ad sharp iequalities ivolvig the costat e ad the sequece + / are proved.. INTRODUCTION The costat e or
More informationPhysics 116A Solutions to Homework Set #9 Winter 2012
Physics 116A Solutios to Homework Set #9 Witer 1 1. Boas, problem 11.3 5. Simplify Γ( 1 )Γ(4)/Γ( 9 ). Usig xγ(x) Γ(x + 1) repeatedly, oe obtais Γ( 9) 7 Γ( 7) 7 5 Γ( 5 ), etc. util fially obtaiig Γ( 9)
More informationA Simplified Binet Formula for k-generalized Fibonacci Numbers
A Simplified Biet Formula for k-geeralized Fiboacci Numbers Gregory P. B. Dresde Departmet of Mathematics Washigto ad Lee Uiversity Lexigto, VA 440 dresdeg@wlu.edu Zhaohui Du Shaghai, Chia zhao.hui.du@gmail.com
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationA symbolic approach to multiple zeta values at the negative integers
A symbolic approach to multiple zeta values at the egative itegers Victor H. Moll a, Li Jiu a Christophe Vigat a,b a Departmet of Mathematics, Tulae Uiversity, New Orleas, USA Correspodig author b LSS/Supelec,
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More informationMechanical Quadrature Near a Singularity
MECHANICAL QUADRATURE NEAR A SINGULARITY 215 Mechaical Quadrature Near a Sigularity The purpose of this ote is to preset coefficiets to facilitate computatio of itegrals of the type I x~^fix)dx. If the
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationarxiv: v3 [math.nt] 24 Dec 2017
DOUGALL S 5 F SUM AND THE WZ-ALGORITHM Abstract. We show how to prove the examples of a paper by Chu ad Zhag usig the WZ-algorithm. arxiv:6.085v [math.nt] Dec 07 Keywords. Geeralized hypergeometric series;
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationNumerical integration of analytic functions
Numerical itegratio of aalytic fuctios Gradimir V. Milovaović, Dobrilo Ð Tošić, ad Miloljub Albijaić Citatio: AIP Cof. Proc. 1479, 146 212); doi: 1.163/1.4756325 View olie: http://dx.doi.org/1.163/1.4756325
More informationlim za n n = z lim a n n.
Lecture 6 Sequeces ad Series Defiitio 1 By a sequece i a set A, we mea a mappig f : N A. It is customary to deote a sequece f by {s } where, s := f(). A sequece {z } of (complex) umbers is said to be coverget
More informationContinued Fractions and Pell s Equation
Max Lah Joatha Spiegel May, 06 Abstract Cotiued fractios provide a useful, ad arguably more atural, way to uderstad ad represet real umbers as a alterative to decimal expasios I this paper, we eumerate
More information