1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
function [x] = back(U,Z) x=zeros(1,n); for i = n:-1:1 sumk=0 for j=i+1:n sumk=sumk+U(i,j)*x(j); end; x(i)=(Z(i)-sumk)/U(i,i); end; endfunction