path: root/plot.py

import numpy as np
import matplotlib.pyplot as plt
from math import e

t = np.arange(0.001, 1, 0.001)
u = np.sqrt((4*e**2 + 4*e +1)*t**2+ (24*e**2-4*e+4)*t+ 36*e**2-24*e+4)/(2*(e-1)*t)-(3*e-1)/((e-1)*t)
print u[0]
print t[0]
plt.plot(t,u)
plt.show()