Euler算法

import numpy as np
import matplotlib.pyplot as plt
x = 0
y = 1
a = np.ones(100)
b = np.ones(100)
t = 0
for i in range(1,100):
    t = y * 1.1 - 0.2 * (x/y)
    x = x + 0.1
    y = t
    
    a[i]=x
    b[i]=y
    

    
plt.plot(a,b)    
plt.show()