import numpy as np
import math
import random
sita0=np.random.random()
sita1=np.random.random()
sita2=np.random.random()
sita=np.array([sita1,sita2])
aerfa=0.0000001
x=np.array([[2104,3],[1600,3],[2400,3],[1416,2],[3000,4]])
t=np.array([400,330,369,232,540])
eps=1e-4
e0=9
e1=9
e2=9
for i in range(1):
i = 0
e0 = (np.sum((x[i, 0:1] * sita[i])) + sita0 * 1 - t[i])
e1 = e0 * x[i, 0]
e2 = e0 * x[i, 1]
sita0 = sita0 - aerfa * e0
sita1 = sita1 - aerfa * e1
sita2 = sita2 - aerfa * e2
i += 1
if(e0>=eps or e1>=eps or e0>=eps):
print (sita0,sita1,sita2)