%Script to test function mygauss.

trials =20

error = zeros(trials,1);

normResidual=zeros(trials,1);

for i=1:trials
 vec = [1:i];
 a = vander(vec);
 one = ones(i,1);
 b = a * one;
 x = mygauss(a,b);
 error(i,1) = norm(x-one);
 residual = a*x- a*one;
 normResidual(i,1) = norm(residual);
end;

error
normResidual