| Author |
Comment/Response |
djeong
|
 01/10/06 11:47pm
Hi,
I am using mathematica 5.2 student edition on the linux machine.
I entered
In[1]:= Om=0.27
In[2]:= Ov=1-Om
In[3]:= f[a_]:=Integrate[(Om/x+Ov*x^2)^(-3/2),{x,0,a}]
In[4]:= f[1]
When I did it in my mathematica, it gave me
Out[4]=0.943314+6.59085 10^-16 I
However, if I put the exactly same thing in my friend's mathematica (v.5.0 on window), I got
Out[4]=1.12594
Also, when I use "NIntegrate" instead of "Integrate", it gives 1.12594.
How can I cure this problem? Should I use only "NIntegrate" instead of "Integrate" for the definite integral in the future? Or is it just a bug for this particular integral?
Thank you in advance,
DJ
P.S.
Actually, v.5.0 can find the solution like this :
In[1]:=Integrate[(Om/x + Ov x^2)^(-3/2), {x, 0, a}, Assumptions -> {0 < a < 1}]
Out[1]=2.851112440442597*a^(5/2)*Hypergeometric2F1[5/6, 3/2, 11/6, -2.7037037037037033*a^3]
URL: , |
|
