from sympy import *
init_printing()
X = Matrix([[1, 10],
[1, 8],
[1, 13],
[1, 9],
[1, 11],
[1, 14],
[1, 6],
[1, 4],
[1, 12],
[1, 7],
[1, 5]])
y = Matrix([[8.04],
[6.95],
[7.58],
[8.81],
[8.33],
[9.96],
[7.24],
[4.26],
[10.84],
[4.82],
[5.68]])
X
y
v = Inverse(X.T * X) * X.T * y
v
#v = (b, m), or (intercept, slope)
#So, the best fit line is given by (approximately) y = 0.5x + 3