  We are presenting a method of linear regression based on Gram-Schmidt
orthogonal projection that does not compute a pseudo-inverse matrix. This is
useful when we want to make several regressions with random data vectors for
simulation purposes.
