nongauss.m

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

%用不选主元的gauss消去法对矩阵进行LU分解，nongauss.m　
function [L,U]=nongauss(A)
L=[];
U=[];
n=length(A(:,1));
for k=1:n-1
    A(k+1:n,k)=A(k+1:n,k)/A(k,k);
    A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);
end
L=tril(A);
U=triu(A);
for i=1:n
    L(i,i)=1;
end