       Re: General--Different Results

• To: mathgroup at smc.vnet.net
• Subject: [mg69445] Re: General--Different Results
• From: dimmechan at yahoo.com
• Date: Wed, 13 Sep 2006 04:00:56 -0400 (EDT)
• References: <ee64m2\$7jl\$1@smc.vnet.net>

```Hello.

***My version of Mathematica does not recognise the character Ïµ,
below. So I replace it with m.

-((4*Pi*Subscript[Ïµ, 0]*n^2*Subscript[k, 0]*A^2)*Abs[Î½]*
Integrate[r^2*BesselJ[Î½, Îº*r]*BesselJ[Î½ - 1, Îº*r],
{r, 0, a}])/((Îº*c)*Î½)

\$Version
5.2 for Microsoft Windows (June 20, 2005)

***Evaluating the integral I get a conditional result (as it shoulb be
since your integral depends on parameters)

-((4*Pi*Subscript[m, 0]*n^2*Subscript[k, 0]*A^2)*Abs[Î½]*
Integrate[r^2*BesselJ[Î½, Îº*r]*BesselJ[Î½ - 1, Îº*r],
{r, 0, a}])/((Îº*c)*Î½) //Timing

{30.141000000000002*Second, (-4*a*A^2*n^2*Pi*Abs[Î½]*If[Im[1/(a*Îº)] >=
1 || 1 + Im[1/(a*Îº)] <= 0 || Im[1/(a*Îº)] == 0 ||
Re[1/(a*Îº)] != 0, (a*(a*Îº)^(2*Î½)*HypergeometricPFQ[{1/2 + Î½},
{2*Î½, 2 + Î½}, -(a^2*Îº^2)])/
(4^Î½*Îº*Gamma[Î½]*Gamma[2 + Î½]), Integrate[a^2*r^2*BesselJ[-1 +
Î½, a*r*Îº]*BesselJ[Î½, a*r*Îº], {r, 0, 1},
Assumptions ->  !(Im[1/(a*Îº)] >= 1 || 1 + Im[1/(a*Îº)] <= 0 ||
Im[1/(a*Îº)] == 0 || Re[1/(a*Îº)] != 0)]]*
Subscript[k, 0]*Subscript[m, 0])/(c*Îº*Î½)}

parameters.

Cheers
Dimitris

aerolapo at hotmail.com wrote:
> I'va got problem in solving this equation,
>
> -((4*Pi*Subscript[&#1013;, 0]*n^2*Subscript[k, 0]*A^2)*Abs[&#957;]*
>
>   Integrate[r^2*BesselJ[&#957;, &#954;*r]*BesselJ[&#957; - 1, &#954;*r],
>    {r, 0, a}])/((&#954;*c)*&#957;)
>
> because sometimes (the first two or three times I calculate it after opening Mathematica) the result is:
>
> -((4^(1 - &#957;)*a^3*A^2*n^2*Pi*(a*&#954;)^(2*&#957; - 1)*abs[&#957;]*
>
>    HypergeometricPFQ[{&#957; + 1/2}, {2*&#957;, &#957; + 2}, (-a^2)*&#954;^2]*
>    Subscript[k, 0]*Subscript[&#1013;, 0])/(c*&#954;*&#957;*Gamma[&#957;]*
>    Gamma[&#957; + 2]))
>
> and some other times: (after two or three identical calculation the result change and it is always this one)
>
> -((2*a^2*A^2*n^2*Pi*Abs[&#957;]*(&#957;*BesselJ[&#957;, a*&#954;]^2 -
>
>     (&#957; - 1)*BesselJ[&#957; - 1, a*&#954;]*BesselJ[&#957; + 1, a*&#954;])*
>    Subscript[k, 0]*Subscript[&#1013;, 0])/(c*&#954;^2*&#957;))
>
> Note that I run exaclty the same procedure in both cases (only the expression on the notebook, I press shift+enter, but if I cancel the result and retype shift+enter it sometimes changes).
>
> I don't know if the two of them are equivalent, but, is it possible there is a mistake and only one of them is correct? Is there a method to test the equivalence? If they are equivalent, which is the rule followed by Mathematica to choose the result to return?
>
> Thank you very much.
>
> Jacopo
>
> Link to the forum page for this post:
> http://www.mathematica-users.org/webMathematica/wiki/wiki.jsp?pageName=Special:Forum_ViewTopic&pid=13453#p13453
> Posted through http://www.mathematica-users.org [[postId=13453]]

```

• Prev by Date: Re: Locating common subexpressions
• Next by Date: Re: Evaluating a Meijer G-function
• Previous by thread: Re: General--Different Results
• Next by thread: Re: Re: Color names and the 10 elementary colors?