[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
>
>
Prev by Date:
**Re: Dot or Inner ... but not quite**
Next by Date:
**Re: Dot or Inner ... but not quite**
Previous by thread:
**Re: Double Integration involving Struve and Neumann functions**
Next by thread:
**Re: Double Integration involving Struve and Neumann functions**
| |