import pylab

r0 = int(10e15)
p0 = 1.7e-14
r = r0
rfinal = int(r*10000)

def p(r):
    return p0*pow(r/r0, -3)

x = list()
y = list()

while (r < rfinal):
    x.append(r)
    y.append(p(r))
    r = r+10e15
    print len(x)

pylab.plot(x,y)
pylab.xlim([10e15,10e17])
pylab.ylim([0,10e-19])
pylab.show()

t0 = int(10e7)
r = r0

def t(r):
    return t0*pow(r/r0, -2)

w = list()
z = list()

while (r < rfinal):
    w.append(r)
    z.append(t(r))
    r = r+10e15
    print len(w)

pylab.plot(w,z)
pylab.xlim([10e15,10e17])
pylab.ylim([0,10e7])
pylab.show()