%%
%%A = lap2D(10,10);

A = fd3d(20,10,1,0,0,0,0);
n = size(A,1);
%%
b = A*ones(n,1);
%%-------------------- omega = 1.56
x = randn(n,1);
om = 1.6;
[xout, nrm1] = sorRelax (A, x, b, 200, om,1.e-06);
%%-------------------- omega = 1.0
x = randn(n,1);
om = 1.0;
[xout, nrm2] = sorRelax (A, x, b, 200, om,1.e-06);
%%-------------------- plot
semilogy([1:length(nrm1)], nrm1,'r--')
hold
semilogy([1:length(nrm2)], nrm2,'b')
legend('om=1.6','om=1.0')
