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