我为Matrix向后计算编写了一些代码,但由于某种原因,它从3变为2,但不会从3变为1。如果我告诉它从3变为1,它将完全跳过指令。
function [x] = backwards(U,b)
[n,n]=size(U);
x(n)=b(n)/U(n,n);
for i=n-1:1
sum=0;
for k=i+1:n
sum=sum+(U(i,k)*x(k));
x(i)=(b(i)-sum)/U(i,i);
end
end
fprintf('solution:');
fprintf('\n');
disp(x');
end