Files
unisa_analisi_numerica_2014…/liner_system_resolver.m

17 lines
240 B
Matlab

U=[2,2,4;0,-7,-11;0,0,2];
b=[5;-8;-2];
x=[];
n= length(b);
x(n) = b(n)/U(n,n);
for i=n-1:-1:1
somma = 0;
for k=i+1:n
somma=somma+U(i,k) * x(k);
end
x(i) = (b(i) - somma)/ U(i,i);
end
disp ('La soluzione è ')
x