for k=1:m,
        if k ~= i,
                      ratio=a(k,t)/a(i,t);
                      a(k,:)=a(k,:)-a(i,:)*ratio;
        end
end
a(i,:)=a(i,:)/a(i,t);

