[Date Index] [Thread Index] [Author Index]

Re: Double Integration involving Struve and Neumann functions

• To: mathgroup at smc.vnet.net
• Subject: [mg79633] Re: Double Integration involving Struve and Neumann functions
• From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
• Date: Wed, 1 Aug 2007 04:50:15 -0400 (EDT)
• Organization: Uni Leipzig
• References: <f8n2v8$lsl$1@smc.vnet.net>
• Reply-to: kuska at informatik.uni-leipzig.de

Hi,

and that phi(ze) is not the proper Mathematica syntax and you mean

psi[ze_]:=Sqrt[2/Le] Cos[(Pi ze)/Le];
psi[zh_]:=Sqrt[2/Lh] Cos[(Pi zh)/Lh];

does not matter ??

Please post your code in correct syntax ..
Than try a single integral and if Mathematica does not find it
read carefull
http://www.amazon.de/Table-Integrals-Products-Academic-Press/dp/0123736374/ref=sr_1_1/028-4565458-4869346?ie=UTF8&s=books-intl-de&qid=1185884008&sr=1-1

if this does not help try some approximation method using series 
expansions ..

Regards
   Jens
Sooraj R wrote:
> Hi all
> I need to evaluate an integral of the form given below
> 
> Eb=hbar^2/(2 mu R^2) - 4 charge^2/(epsilon R^2)
> Integrate[(psi(ze))^2 (psi(zh))^2 (Abs[ze - zh])
> ((Pi/2) (H1[2 Abs[(ze - zh)]/R] - 
>            Nu1[2 Abs[(ze - zh)]/R]) - 1), {ze, -Le/2,
> Le/2}, {zh, -Lh/2, 
>    Lh/2}]
> 
> I define
> Eb=Table[0,{16}];
> charge = 1.602 10^-19;
> hbar = 1.054571 10^-34;
> m0 = 9.1019 10^-31;
> mu=0.05 m0;
> epsilon=8.85418 10^-12
> Le =
> {100,105,110,115,120,125,130,135,140,145,150,155,160,165,170};
> Lh =
> {82,87,92,97,102,107,112,117,122,127,132,137,142,147,152,157};
> 
> psi(ze)=Sqrt[2/Le] Cos[(Pi ze)/Le]; and 
> psi(zh)=Sqrt[2/Lh] Cos[(Pi zh)/Lh];
> are the wave functions of the system.
> R is a parameter and here it can be treated as a
> constant.(anyway the result will be a function of R)
> 
> Le and Lh consists of lists..each having 16 values.
> I need to evaluate the integral for these 16 values of
> Le and Lh.
> H1 and Nu1 are the first order Struve functions and
> Neumann functions.
> For ze, the limit of integration is from -Le/2 to Le/2
> and for zh, the limit is from -Lh/2 to Lh/2.
> 
> I coded the integral as follows.
> 
> 
> For[i = 1, 1 < 17, i++, 
>  Eb[[i]] = 
>   hbar^2/(2 mu R^2) - 
>    4 charge^2/(epsilon R^2) Integrate[(Sqrt[
>           2/Le[[i]]] Cos[(Pi ze)/Le[[i]]])^2 (Sqrt[
>           2/Lh[[i]]] Cos[(Pi zh)/Lh[[i]]])^2 (Abs[
>         ze - zh]) ((Pi/2) (StruveH[1, 2 Abs[(ze -
> zh)]/R] - 
>            BesselY[1, 2 Abs[(ze - zh)]/R]) - 1), {ze,
> -Le[[i]]/2, 
>       Le[[i]]/2}, {zh, -Lh[[i]]/2, Lh[[i]]/2} // N]]
> 
> 
> Mathematica was evaluating it for almost two
> days...and I had to quit the kernel to stop it.(Abort
> evaluation was not working out)
> 
> I was thinking the problem is because of the For loop,
> so decided to check whether it is able to evaluate for
> a particular value of Le and Lh.
> 
> To my sorrow, I found that after running for almost 24
> hrs, mathematica returned the integral unevaluated.
> 
> Can someone help me how to solve this integral?. The
> Cosine wave functions are the simplest case, and I
> need to calculate Eb for systems having Airy functions
> as the wavefunctions, and I'm scared to imagine such a
> situation.
> 
> Any help is deeply appreciated.
> Many thanks in advance
> Best regards
> 
> Sooraj
> 
> 
> 
> 
> 
> 
>       ____________________________________________________________________________________
> Fussy? Opinionated? Impossible to please? Perfect.  Join Yahoo!'s user panel and lay it on us. http://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7 
> 
>