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

L = tril(A,-1);
D = diag(diag(A));
DL = D + omega*L;
f = omega * ( DL \ b );
U = triu(A,1);
xk = zeros((N-1)^2,1); err = []; res = [];
for k = 1:400,
   xk = DL \ ( (1-omega)*D*xk - omega*U*xk ) + f;
   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),