function [x1,iter] = template_im_1(a,b,alpha,x0,tol,maxiter)
n=length(b);
x1=zeros(n,1);
b=alpha*b;
a=alpha*a;
a=speye(n)-a;
for iter=1:maxiter
    x1=a*x0+b;
    err=norm(x1-x0, inf)/norm(x1, inf);
    if(err<=tol)
        disp('convergenza entro il numero di iterazioni');
        break
    end
    x0=x1;
end
end