forward_col.m

来自「various matlab programs to slove various」· M 代码 · 共 12 行

function [b]=forward_col(L,b)
[n]=mat_square(L);
l = min(diag(abs(L))); 
if l == 0
  disp('Error: the matrix is singular'); b = [];
  break
end
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);