downsan.m

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

%求解下三角方程组的通用子函数，downsan.m
function  b=downsan(L,b)
n=length(L(:,1));
for j=1:n-1
    b(j)=b(j)/L(j,j);
    b(j+1:n)=b(j+1:n)-b(j)*L(j+1:n,j);
end
    b(n)=b(n)/L(n,n);