U=[1,2,3:4,5,6;7,8,9];
B=[5;6;7];

x=[];
n= lenght(B);
x(n) = b(n)/U(n,n);

for i=n-1:-1:1
    somma = 0;
    for k