forw_band.m

various matlab programs to slove various problems

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