function [l,u]=GE(A)
r=rank(A);
u=A;
l=eye(r);
   for p=1:(r-1)
       for m=(p+1):r
           l(m,p)=u(m,p)/u(p,p);
           for n=1:r
               y=u(m,n)-l(m,p)*u(p,n);
               u(m,n)=y;
           end
       end
  end