import numpy as np import matplotlib.pyplot as plt alpha =0.1 x = np.array([[2104,3],[1600,3],[2400,3],[1416,2],[3000,4]]) t = np.array([400,300.369,232,540]) c0 = np.random.random() c1 = np.random.random() c2 = np.random.random() c = np.array([c1,c2]) epslion = 1e-4 e0 = 10 e1 = 10 e2 = 10 while e0 >= epslion or e1 >= epslion or e2 >= epslion: i = 0 e0 = (np.sum((x[i] * c[i])) + c0 * 1 - t[i]) e1 = e0 * x[1, i] e2 = e0 * x[2, i] c0 = c0 - alpha * e0 c1 = c1 - alpha * e1 c2 = c2 - alpha * e2 i += 1 print(c0,c1,c2)