clear all
n=input('dimensione della matrice='); 
a1=2*eye(n);
 
a1(1,2:n)=0.2; 
a1(2:n,1)=0.2 
[l1,u1,p1]=lu(a1) 



a2=2*eye(n);
a2(n,1:n-1)=0.2;
a2(1:n-1,n)=0.2
[l2,u2,p2]=lu(a2)