import numpy as np
import matplotlib.pyplot as plt
x = np.array([[2104, 3], [1600, 3], [2400, 3], [1416, 2],[3000,4]])
t = np.array([400, 330, 365, 232, 540])
a = 0.1
b0 = np.random.random()
b1 = np.random.random()
b2 = np.random.random()
b = np.array([theta1, theta2])
eps = 1e-5
t0 = 999
t1 = 999
t2 = 999
while t0 >= eps or t1 >= eps or t2 >= eps:
i = 0
t0 = (np.sum((x[:, i] * b[i])) + b0 * 1 - t[i])
t1 = t0 * x[1, i]
t2 = t0 * x[2, i]
b0 = b0 - a * t0
b1 = b1 - a * t1
b2 = b2 - a * t2
i += 1
print(b0, b1, b2)
11153087938.716803 17844940701444.758 26767411052167.195