for k = 1:20
    k
    sizes(k) = k^3;
    [A,b]=gendata3(k^3,1);
    x = A\b;
    condition(k) = condest(A);
    norm1(k) = norm(x,1);
    norm2(k) = norm(x);
    resid(k) = norm(A*x-b);
end;