forw_band.m

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

function [b] = forw_band (L, p, b)
[n,n]=size(L);
for j = 1:n
   for i=j+1:min(j+p,n); 
      b(i) = b(i) - L(i,j)*b(j);
   end
end