2021-9-24

import numpy as np
import matplotlib.pyplot as plt

x=np.zeros(100)
y=np.zeros(100)

x[0]=0
y[0]=1

y_org=np.zeros(100)
y_org[0]=1

for i in range(1,100):
k_org=(y[i-1]-2x[i-1]/y[i-1])
y[i]=y[i-1]+0.1k_org
x[i] = 0.1 +x[i-1]
k_now = (y[i] - 2 * x[i] / y[i])
y[i]=y[i-1]+0.1*(0.5*(k_org+k_now))

for i in range(100):
y_org[i]=y[i-1]+0.1*(y[i-1]-2*x[i-1]/y[i-1])
x[i] = 0.1 + x[i - 1]

plt.subplot(1,2,1)
plt.scatter(x,y)
plt.subplot(1,2,2)
plt.scatter(x,y_org)
plt.show()