double loop
- To: mathgroup at smc.vnet.net
- Subject: [mg113006] double loop
- From: maria giovanna dainotti <mariagiovannadainotti at yahoo.it>
- Date: Sun, 10 Oct 2010 06:43:47 -0400 (EDT)
Dear All,
I cleared up the notebook but when I plot the picture I don't recover the
contour plot desidered. The colours should be the different correlation
coefficient and it is not the case.
I will be grateful if you could help
Maria
Â
h=0.71
ckms=299792.5
HubE[z_,Om_,Ol_]:=SQRT(Om*(1+z)^3+Ol)dL[z_,Om_,Ol_,h_]:=3.0856*1018*106*ckms/(100*h)*(1+z)*NIntegrate[1/HubE[z,Om,Ol],{z,0,z},MaxRecursion
->100]
TabOmegalambda=Table[i,{i,0.1,1.1,0.05}]
TabOmegaM=Table[i,{i,0.0,1.0,0.05}]
LxOmegalambda={};
Do[
Ol=TabOmegalambda[[k]];
LxOmegaM={};
Do[
Om=TabOmegaM[[j]];
DataGood={};
Do[
idGRB=Dataspectrum[[i,1]];
z=Dataspectrum[[i,2]];
Espectrum=Dataspectrum[[i,3]];
Fx=Dataspectrum[[i,4]];
beta=Dataspectrum[[i,5]];
EisonewOm=4*p*dL[z,Om,Ol,h]2*(1+z)-2*Espectrum;
Lx=4*p*dL[z,Om,Ol,h]2*(1+z)-(1+beta)*Fx;
AppendTo[DataGood,{EisonewOm,Lx,z,idGRB,Om,Oï?¬}],{i,1,Length[Dataspectrum]}];
CLxEiso=N[Correlation[DataGood[[All,1]],DataGood[[All,2]]]];
AppendTo[LxOmegaM,{Om,Ol,CLxEiso}],{j,1,Length[TabOmegaM]}];
If[Ol+Om<2.2,AppendTo[LxOmegalambda,{LxOmegaM}]],{k,1,Length[TabOmegalambda]}]
Flatten[LxOmegalambda,2]
ListContourPlot[Flatten[LxOmegalambda,2],PlotRange
->All,FrameLabel
->{"Omega_M ","Omega_lambda"}]