In [24]:

    

ind=100
r,M,tB,E=np.genfromtxt('model_fun.dat',unpack=True)
t,t1=np.genfromtxt('t.dat',unpack=True)

R_r_a = np.genfromtxt('R_r_a.dat')
R_r = np.genfromtxt('R_r.dat')

clf()
plot(t[:],R_r_a[:,ind],'b')
plot(t[:],R_r[:,ind],'r+')

a=R_r_a[:,ind]
b=R_r[:,ind]


    

In [25]:

    

ind=100
r,M,tB,E=np.genfromtxt('model_fun.dat',unpack=True)
t,t1=np.genfromtxt('t.dat',unpack=True)

R_t_a = np.genfromtxt('R_t_a.dat')
R_t = np.genfromtxt('R_t.dat')

clf()
plot(t[:],R_t_a[:,ind],'b-')
plot(t[:],R_t[:,ind],'r+')

a=R_r_a[:,ind]
b=R_r[:,ind]
ratio = a/b
print "ratio_max={0}".format(ratio.max())
print "ratio_min={0}".format(ratio.min())
#plot(ratio)


    

    




ratio_max=1.0
ratio_min=1.0

In [26]:

    

ind=100
r,M,tB,E=np.genfromtxt('model_fun.dat',unpack=True)
t,t1=np.genfromtxt('t.dat',unpack=True)

R_t_a = np.genfromtxt('R_t_a.dat')
R_rt_a = np.genfromtxt('R_rt_a.dat')
R_rt = np.genfromtxt('R_rt.dat')

clf()
#plot(t[:],R_t_a[:,ind],'+')
plot(t[:],R_rt_a[:,ind],'b-')
plot(t[:],R_rt[:,ind],'r+')

a=R_rt_a[:,ind]
b=R_rt[:,ind]
ratio=a/b
#plot(ratio)
#duza ruznica dla t bliskiego 0
print "ratio_max={0}".format(ratio.max())
print "ratio_min={0}".format(ratio.min())

#a=R_t_a[:,ind]
#b=R_rt_a[:,ind]
#ratio=a/b
#plot(ratio)


    

    




ratio_max=1.07399758793
ratio_min=1.0

In [3]:

    

expk0 = np.genfromtxt('exp_k0.dat')
exp = np.genfromtxt('exp_a.dat')
t1,t = np.genfromtxt('t.dat',unpack=True)
ind=100
clf()
plot(t[:],expk0[:,ind],'b-')
plot(t[:],exp[:,ind],'r-')


    

    
Out[3]:





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

In [4]:

    

rok0 = np.genfromtxt('rho_k0.dat')
ro = np.genfromtxt('rho_a.dat')
t1,t = np.genfromtxt('t.dat',unpack=True)
ind=100
clf()
plot(t[:],rok0[:,ind],'b-')
plot(t[:],ro[:,ind],'r-')


    

    
Out[4]:





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

In [ ]: