Browse code

6.5 works

Michi authored on03/06/2012 17:33:04
Showing1 changed files
... ...
@@ -242,7 +242,7 @@ for (i in (1:values_A)){
242 242
 		integration_MA[i]=I
243 243
 	}
244 244
 
245
-integration_MB=array(0, c(values_B, values_B))
245
+integration_MB=array(0, c(values_B+1, values_B+1))
246 246
 for (i in (1:values_B)){
247 247
 	for (ii in (1:values_B)){
248 248
 		I=Myintegrate(give_function_MB(i,ii),15.5,16.5)
... ...
@@ -250,5 +250,5 @@ for (i in (1:values_B)){
250 250
 	}
251 251
 }
252 252
 
253
-contour(Omegam_vec_B,w_vec,integration_MB)
253
+contour(Omegam_vec_B,w_vec,integration_MB,xlab='Omegam',ylab='w',col = "blue")
254 254