Fit y^=Xw+b\hat{y} = Xw + by^=Xw+b by minimising mean squared error with batch gradient descent.
Implement fit(X, y, lr=0.01, epochs=1000) that starts from w = 0, b = 0 and returns (w, b) where w is a NumPy array of shape (n_features,).
fit(X, y, lr=0.01, epochs=1000)
w = 0
b = 0
(w, b)
w
(n_features,)