Re: Could someone run this on Mathematica 7.0?
- To: mathgroup at smc.vnet.net
- Subject: [mg96748] Re: [mg96726] Could someone run this on Mathematica 7.0?
- From: Aaron Fude <aaronfude at gmail.com>
- Date: Mon, 23 Feb 2009 05:04:45 -0500 (EST)
- References: <200902220813.DAA15166@smc.vnet.net>
Thank you!
On Sun, Feb 22, 2009 at 4:18 AM, Adriano Pascoletti <
adriano.pascoletti at dimi.uniud.it> wrote:
> 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>
- Could someone run this on Mathematica 7.0?