testlinsolve.m

Solve linear system and test it.

function s=textlinsolve(x,y)
s=0; z=0;
while(z<1000)
    A=randn(x,y);
    b=randn(x,1);
    if (norm(A*(lin_solve(A,b))-b)<norm(A*(A\b)-b)) s=s+1;
    end
    z=z+1;
end
s/z