Create linfit.py

Fit a straight line `f(x) = a + bx` to points `(x, y)`.  Returns  coefficients `a` and `b` that minimize the squared error.