function [u,it] = Richardson(A,b,u0,prec,omega,itMax)

u = u0;
r0 = b-A*u;
normr0 = norm(r0);
r = r0;
it = 0;

while norm(r)>prec*normr0 && it<itMax
    u = u + omega*r;
    r = b-A*u;
    it = it+1;
end

