function x=esvschritt(A, b, x)
% macht einen Schritt des ESV, SEHR naiv programmiert
[n m]=size(A);
for i=1:n
    x(i,1)=(b(i,1)-A(i,1:i-1)*x(1:i-1,1)...
        -A(i,i+1:end)*x(i+1:end,1))/A(i,i);   
end


