upsan.m

来自「1、 写出此两点边值问题的积分形式； 2、 构造有限元方程； 3、 编程计算」· M 代码 · 共 8 行

%求解上三角方程组的通用子函数，upsan.m
function y=upsan(U,y)
n=length(U(:,1));
for j=n:-1:2
    y(j)=y(j)/U(j,j);
    y(1:j-1)=y(1:j-1)-y(j)*U(1:j-1,j);
end
y(1)=y(1)/U(1,1);