The CoeffLinear procedure

This procedure returns a matrix $U$ of dimension n $\times$ m+1 with $U(i,j)$ equal to the coefficient of $x_j$ in the equation $i$ while $U(i,{\tt m}+1)$ is the constant term of the equation. In the previous example this procedure will be:

 
void CoeffLinear(MATRIX &U)
{
U(1,1)=1; U(1,2)=0.5; U(1,3)=0.3;
U(2,1)=2; U(1,2)=-0.5; U(1,3)=0.2;
}