import numpy as np import matplotlib.pyplot as plt k = np.array([1]) x=[0] for i in range(1,101,1): a = [i/10] x = np.append(x,a) #x = np.array([0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2]) y = np.array([1]) for i in range(0, 100): k1 = k[i] y1 = 1.1 * y[i] - 0.2 * x[i] / y[i] k2 = y1 - 2 * x[i+1] / y1 k3 = np.array([(k1 + k2) / 2]) k = np.append(k, k3) y2 = np.array([y[i] + k[i+1] * 0.1]) y = np.append(y, y2) plt.plot(x,y,'g-s') plt.show()