Integrals of Bessel functions

*To*: "mathgroup" <mathgroup at yoda.ncsa.uiuc.edu>*Subject*: Integrals of Bessel functions*From*: "WILLIAM M. GOLDING" <golding at ccf4.nrl.navy.mil>*Date*: 21 Jan 91 16:39:00 EST

> Date: Wed, 9 Jan 91 9:26:57 PST > From: Bad Religion <tcrevier at jarthur.Claremont.edu> > To: mathgroup at yoda.ncsa.uiuc.EDU > Subject: Mathematica mailing list > Anyone know of a way to get mathematica to symbollically integrate > Bessel J functions? > -Tom Tom, It's certainly possible to get Mathematica to integrate symbolically any of the special functions, but you have to be able to add rules to the Integrate function. This is done with the sequence, Unprotect[Integrate] ... Add your definitions here ... Protect[Integrate] The next problem is what rules to add. That involves some work. I'd suggest writing down one particular integral involving Bessel functions that you're interested in and trying to get Mathematica to evaluate it correctly. Once you've done that, generalizations to other particular cases should be straight forward. You should also refer to pages 231-234 and 425-426 of the Mathematica book for examples of defining your own integrate function and of modifying the built-in Integrate function. Here is an example allowing you to do simple indefinite integrals of BesselJ functions. The definitions for StruveH are given in the package struve.m which is available in the archives of this site in the files /Symbolic/Mathematica/Packages/Programming_in_Mathematic1 and ...2 These files are the listings of the routines of the book by Roman Maeder. For testing, I just have all these commands placed in a file called bessint.m which can be input directly into Mathematica. <<"rmpack\\struve.m" (* This should refer to struve.m on your system *) Unprotect[Integrate] Integrate[BesselJ[0,z_],z_]:=z BesselJ[0,z] + Pi z/2 ( StruveH[0,z] BesselJ[1,z] - StruveH[1,z] BesselJ[0,z] ) Integrate[BesselJ[1,z_],z_]:= - BesselJ[0,z] Integrate[BesselJ[n_Integer,z_],z_]:= Integrate[BesselJ[n-2,z],z] - 2 BesselJ[n-1,z] /; n > 1 Protect[Integrate] Mathematica (MS-DOS 386/7) 1.2 (September 27, 1989) [With pre-loaded data] by S. Wolfram, D. Grayson, R. Maeder, H. Cejtin, S. Omohundro, D. Ballman and J. Keiper with I. Rivin, D. Withoff and T. Sherlock Copyright 1988,1989 Wolfram Research Inc. In[1]:= <<bessint.m (* the file of definitions described above *) Out[1]= {Integrate} (* Here is a sample of some of the lower order integer order bessel functions *) In[2]:= Integrate[BesselJ[0,x],x] Out[2]= x BesselJ[0, x] + Pi x (BesselJ[1, x] StruveH[0, x] - BesselJ[0, x] StruveH[1, x]) > ---------------------------------------------------------------- 2 In[3]:= Integrate[BesselJ[1,x],x] Out[3]= -BesselJ[0, x] In[4]:= Integrate[BesselJ[2,x],x] Out[4]= x BesselJ[0, x] - 2 BesselJ[1, x] + Pi x (BesselJ[1, x] StruveH[0, x] - BesselJ[0, x] StruveH[1, x]) > ---------------------------------------------------------------- 2 In[5]:= Integrate[BesselJ[3,x],x] Out[5]= -BesselJ[0, x] - 2 BesselJ[2, x] In[6]:= Integrate[BesselJ[4,x],x] Out[6]= x BesselJ[0, x] - 2 BesselJ[1, x] - 2 BesselJ[3, x] + Pi x (BesselJ[1, x] StruveH[0, x] - BesselJ[0, x] StruveH[1, x]) > ---------------------------------------------------------------- 2 In[7]:= Integrate[BesselJ[5,x],x] Out[7]= -BesselJ[0, x] - 2 BesselJ[2, x] - 2 BesselJ[4, x] (* As I was playing with these Bessel functions I found that mathematica already knew how to do integrals of the half integral order Bessel functions *) In[8]:= Integrate[BesselJ[5/2,x],x] -I -- Sqrt[2] Erf[Sqrt[-I] Sqrt[x]] 2 Sqrt[2] Cos[x] 2 Out[8]= ---------------- + -------------------------------- + Sqrt[Pi] Sqrt[x] Sqrt[-I] I - Sqrt[2] Erf[Sqrt[I] Sqrt[x]] 2 2 Sqrt[2] Sin[x] > ------------------------------ - ---------------- Sqrt[I] 3/2 Sqrt[Pi] x (* But you need to use for instance BesselJ[5/2,x] rather than BesselJ[2.5,x], otherwise you simply get the bessel function replaced this series under the integrate function *) In[9]:= Integrate[BesselJ[2.5,x],x] Out[9]= Integrate[(0.797884560802865 3. Cos[-4.71239 + x] 3. Sin[-4.71239 + x] > (0. + Cos[-4.71239 + x] - -------------------- - --------------------) 2 x x > ) / Sqrt[x], x] Anyways, No guarantees that the above works for all or almost all cases. But I think you see that it is possible to extend the integrate function to include the special functions to some degree. Mike Golding