Services & ResourcesWolfram ForumsMathGroup Archive

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

Could someone run this on Mathematica 7.0?

• To: mathgroup at smc.vnet.net
• Subject: [mg96726] Could someone run this on Mathematica 7.0?
• From: Aaron Fude <aaronfude at gmail.com>
• Date: Sun, 22 Feb 2009 03:13:00 -0500 (EST)

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

• Follow-Ups:
  • Re: Could someone run this on Mathematica 7.0?
    • From: Aaron Fude <aaronfude@gmail.com>
  • Re: Could someone run this on Mathematica 7.0?
    • From: Adriano Pascoletti <adriano.pascoletti@dimi.uniud.it>