Matrix.
Get X /variables = c engine /file = *.
Print X.
Get Y/variables = mpg /file = *.
Print Y.
Compute xtx =T(x)*x.
Compute betahat = inv(xtx)*T(x)*y.
Print betahat.
Compute Yhat = X*betahat.
Print Yhat.
compute R = Y - Yhat.
Print R.
End matrix.