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

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