% utilizza template2 per calcolare in sequenza le norme-1 delle colonne di inv(U)
% e ne prende il massimo, cioè norm(U,1)

alpha=-1+sqrt(eps);
n=100;
maxnorm=0;
for k=1:n
    ek=zeros(n,1);
    ek(k)=1;
    vk=template2(alpha,ek);
    nm=norm(vk,1);
    if nm>maxnorm
        maxnorm=nm;
    end
end
maxnorm
