e2pg313b.m

X0=[1.00 1.00 1.00 1.00 1.00 2.00 2.00 2.00;
    2.00 2.00 4.00 4.00 6.00 2.00 2.00 4.00;
       0 1.00 0 1.00 2.00 0 1.00 0;
   -2.52 -2.71 -8.34 -8.40 -14.60 -0.62 -0.47 -6.49];
X1=[2.00 3.00 3.00 3.00 3.00 3.00 3.00 3.00;
    6.00 2.00 2.00 2.00 4.00 6.00 6.00 6.00;
             0 0 1.00 2.00 1.00 0 1.00 2.00;
   -12.46 1.36 1.40 1.60 -4.64 -10.34 -10.43 -10.30];
Xd=[X0 X1];
al=zeros(2,15); al=['coeff x1      :';'coeff x2      :'];
b=mregg(Xd([1 2 4],:),1,0,al);