In [5]:

    

import background


    

In [6]:

    

x = background.background(.7, .02, 1.0e-5, 0, 10.25)
y = background.background(.7, .02, 1.0e-5, 2, 10.25)
y.z0


    

    
Out[6]:





55.51603069472758

In [7]:

    

y.hubble(5)


    

    
Out[7]:





array(0.0012474669695379948)

In [8]:

    

z = linspace(0,10,10000)

plot(z, [x.hubble(i) for i in z])
plot(z, [y.hubble(i) for i in z])
figure()
plot(z, array([x.hubble(i) for i in z])/array([y.hubble(i) for i in z]))


    

    
Out[8]:





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






    













    

In [9]:

    

x.hubble(5)


    

    
Out[9]:





0.00073100846194857564

In [ ]: