N = 100;
omega = 0.9;
format long
[A,b] = twodlaplace(N);
disp('Konditionszahl von A');
condest(A),
x = A \ b;

D = diag(diag(A)); f = D \ b;
LU = A-D;
disp('Spektralradius von B');
normest(omega*inv(sqrt(D)) * LU * inv(sqrt(D)) + (1-omega)*speye((N-1)^2)),
xk = zeros((N-1)^2,1); err = []; res = [];
for k = 1:1000,
   xk = xk + omega*(D \ ( b-A*xk ) );
   err = [err;norm(x-xk)];
   res = [res;norm(A*xk - b)];
end
semilogy(err,'r-'),
hold on
semilogy(res,'b-'),
hold off
makelarge
err(end) / norm(x),