In [12]:

    

from get_background import get_hubble_rate
rhochi, da, h = get_hubble_rate(**{		'omega_m': 0.135 ,
					'H0':  70,
					'm_chi_over_H0' : 4,
					'chi0' : .5,
					'chidot0': 0,
					'r_s':147.5
					 })
junk, da2, h2 = get_hubble_rate(**{		'omega_m': 0.13924479,
					'H0': 72.2385516,
					'm_chi_over_H0' : 7,
					'chi0' : 0,
					'chidot0': 0,
					'r_s':147.5
					 })


    

In [13]:

    

z = linspace(0.1,10, 10000)
plot(z, da(z)/147.5)
scatter(array([0.35, 0.57, 2.4]), array([6.875, 9.191, 10.8]))
figure()
plot(z, h(z)*3e5)
y = array([12895, 14231, 3e5/9])/147.5
scatter(array([0.35, 0.57, 2.4]), y)
figure()
plot(1/z, rhochi(1/z))


    

    
Out[13]:





[<matplotlib.lines.Line2D at 0x107335c10>]






    













    













    

In [3]:

    

z = linspace(0,10,1000)
plot(z, h(z)/h2(z))
figure()
plot(z, da(z)/da2(z))


    

    




-c:4: RuntimeWarning: invalid value encountered in divide






    
Out[3]:





[<matplotlib.lines.Line2D at 0x106b3dcd0>]






    













    

In [ ]: