import numpy as np
x = np.array([[5, 6], [1, 4], [7, 8], [2, 3]])
t = np.array([5, 1, 2, 6])
alpha = 0.0001
m0 = np.random.random()
m1 = np.random.random()
m2 = np.random.random()
m = np.array([m1, m2])
eps = 1e-4
e0 = 99
e1 = 99
e2 = 99
while e0 >= eps or e1 >= eps or e2 >= eps:
i = 0
e0 = (np.sum((x[:, i] * m[i])) + m0 - t[i])
e1 = e0 * x[1, i]
e2 = e0 * x[2, i]
m0 = m0 - alpha * e0
m1 = m1 - alpha * e1
m2 = m2 - alpha * e2
i += 1
print(m0, m1, m2)
0.7937814248441828 0.14726026635657302 0.9167388517279859