Deriving the closed form solution for multiple linear regression, with a scikit-learn example in code.