basicilu.m

various matlab programs to slove various problems

function [a] = basicILU(a)
[n,n]=size(a);
for k=1:n-1
   for i=k+1:n,
      if a(i,k) ~= 0
        a(i,k) = a(i,k) / a(k,k);
        for j=k+1:n
           if a(i,j) ~= 0
             a(i,j) = a(i,j) -a(i,k)*a(k,j);
           end
         end
      end
   end
end
return