MathGroup Archive 2010

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

Search the Archive

Re: solving an integral

  • To: mathgroup at smc.vnet.net
  • Subject: [mg113487] Re: solving an integral
  • From: "Sjoerd C. de Vries" <sjoerd.c.devries at gmail.com>
  • Date: Sun, 31 Oct 2010 02:06:33 -0500 (EST)

Maria,

One of the problems in your code is the use of SQRT where it should be Sqrt.
Also, in one instance you're using round parentheses instead of square
brackets for the function call. The variable Ok isn't defined either.

Cheers Sjoerd

From: maria giovanna dainotti [mailto:mariagiovannadainotti at yahoo.it] 
Sent: Wednesday 27 October 2010 18:36
To: sjoerd.c.devries at gmail.com
Cc: mathematica group
Subject: [mg113487] 

 

Dear  Sjoerd,

I have a problem in solving the following integral

h=0.71
ckms=299792.5

SQRT[Abs[1-Om-Ol]]>0

1-Om-Ol<0

A=3.0856*10^18*10^6*ckms/(100*h)

Om=Table[i,{i,0.0,1.0,0.01}]

Ol=Table[i,{i,0.0,1.0,0.01}]

z=Table[i,{i,0.2,8.26,0.08}]

 HubE[z_,Om_,Ol_,Ok_]:=SQRT(Om*(1+z)^3+Ol+(1-Om-Ol)(1+z)^2)

dL[z_,Om_,Ol_,h_]:=3.0856*10^18*10^6*ckms/(100*h)*(1+z)*1/(SQRT[Abs[Ok]])NIn
tegrate[1/HubE[z,Om,Ol,Ok],{z,0,z},MaxRecursion

->100]
The integral diverge even though I have putted the right condition not to
make it diverging.

I am very grateful for your help.

Cheers

Maria


 

 


  • Prev by Date: Re: Importing "Plaintext" from PDF
  • Next by Date: Re: Importing "Plaintext" from PDF
  • Previous by thread: solving an integral
  • Next by thread: GraphicsData["PICT" ... in cell does not display