import matplotlib.pyplot as plt
import numpy as np
h=0.01
x = 0
xx=np.zeros(100)
yy=np.zeros(100)
y = 1
for i in range (100):
y=h*(y-2.0*x/y)+y
x=x+h
xx[i]=x
yy[i]=y
print(xx[i],yy[i])
plt.plot(xx,yy, ls='-', lw=2, label='OULA', color='blue')
plt.legend()
plt.show()