clm function

Constrained least squares