function [z]=jacobistep(n,gamma,b,y)
% fa un passo del metodo di Jacobi applicato
% al problema ax=b ove a=In+gamma L
a=faia(n,gamma);
d=diag(a);
dd=diag(d);
cc=a-dd;
di=1./d;
ddi=diag(di);
z=ddi*(b-cc*y);
end
