function [ynew,iter] = my_gs_28_04(x,v,z,b,tol,maxiter) %Jacobi method for arrow matrices n=length(x); yold=zeros(n,1); ynew=zeros(n,1); iter=0; for k=1:n if(x(k)==0) disp('metodo non applicabile'); return end end err=inf; while(err>tol & iter<=maxiter) %calcolare ynew da yold ynew(1)=(b(1)-v'*yold(2:n))/x(1); for k=2:n ynew(k)=(b(k)-z(k-1)*ynew(1))/x(k); end err=norm(ynew-yold,inf)/norm(ynew,inf); iter=iter+1; yold=ynew; end end