function [e,dege] = poleneps(A,degA,i,j,EPS) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% e is (i,j) entry of a polynomial matrix A              %
%                                                        %
% e(s) = A(s)(i,j)                                       %
%                                                        %
% clearing by EPS employed                               % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% M. Sebek, June 1990                                    %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% function used: polsize;


[rA,cA] = polsize(A,degA);
e = A(i,j : cA : cA*degA +j);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% clearing via EPS                                       %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

EPS=norm(A,Inf)*EPS;

for i=1:degA
  if abs(e(i))<=EPS
    e(i)=0;
  end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% setting dege                                           %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

dege=degA;

while e(dege+1)==0
  e(dege+1)=[];
  dege=dege-1;
  if dege<0
    dege=-inf;
    e=0;
    return
  end
end