Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: Could someone run this on Mathematica 7.0?

• To: mathgroup at smc.vnet.net
• Subject: [mg96745] Re: [mg96726] Could someone run this on Mathematica 7.0?
• From: Adriano Pascoletti <adriano.pascoletti at dimi.uniud.it>
• Date: Mon, 23 Feb 2009 05:04:13 -0500 (EST)
• References: <200902220813.DAA15166@smc.vnet.net>

Aaron,
with

In[8]:= $Version
Out[8]= "7.0 for Mac OS X x86 (64-bit) (November 11, 2008)"

I get

32.66443126862095


Adriano Pascoletti

2009/2/22 Aaron Fude <aaronfude at gmail.com>

> Hi,
>
> I'm not sure if such a request is appropiate on this newsgroup, but I
> know that Mathematica changed some aspects of Hypergeometric functions
> from version 6 to version 7. I have version 6 and I want to make sure
> that an error that I am looking for is due to Mathematica. The code is
> (thank you!)
>
> C0[x_][theta_] := Cos[Pi x]/Cos[theta] - 1
> C1[x_] := C0[x]'
> C2[x_] := C1[x]'
>
> psi[r_] := BesselJ[0, BesselJZero[0, 1]* r]/(
>  Sqrt[Pi] BesselJ[1, BesselJZero[0, 1]])
>
> L[x_][r_, theta_] :=
>  1/(1 + C0[x][
>     theta])^2 ((1 + C1[x][theta]^2/(1 + C0[x][theta])^2)*(psi''[r] +
>       1/r psi'[r])
>    - (C2[x][theta] (1 + C0[x][theta]) -
>      C1[x][theta]^2)/(1 + C0[x][theta])^2*1/r psi'[r])
>
> s = L[x][r, x*alpha] // Series[#, {x, 0, 2}] &;
>
> -Integrate[psi[r]*s*r, {alpha, -Pi, Pi}, {r, 0, 1}] +
>   Pi*psi[0]^2 + (1/3 Pi^3)*psi[0]^2 x^2 // Coefficient[#, x, 2] & //
>  N
>
>

• References:
  • Could someone run this on Mathematica 7.0?
    • From: Aaron Fude <aaronfude@gmail.com>