function [A,degA,DEGA,LCA] = polmatce(A,degA,EPS)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Elementwise finds degrees (DEGA), leading coefficients (LCA)%
% and clears A via EPS. If possible, decreases degA.          %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% M. Sebek, September 1990                                    %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% functions used: poleneps, polsize, polseten;

[rA,cA] = polsize(A,degA); 
DEGA = zeros(rA,cA); LCA = zeros(rA,cA);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% creating DEGA and LCA and clearing A                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for i=1:rA
    for j=1:cA
      [a,dega]=poleneps(A,degA,i,j,EPS);
      A=polseten(A,degA,a,dega,i,j);
      DEGA(i,j)=dega;
      if dega<0
	LCA(i,j)=0;
      else
	LCA(i,j)=a(dega+1);
      end
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% decreasing degA                                             %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

m=max(max(DEGA));
if m<0
  A=zeros(rA,cA);
  degA=0;
else
  A(:,cA*(m+1)+1:cA*(degA+1))=[];
  degA=m;
end