function [u,degu,v,degv] = poldiv_e(a,dega,b,degb,EPS)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% polynomial division:                                    % 
%                                                         %
%    a = b*u + v,    degv<degb                            %
%                                                         %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% M. Sebek, June 1990                                     %
% modified by S. Pejchova, August 1994                    %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% function used: polred_e 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if degb<0
   error('poldiv: division by a zero polynomial');
end

if nargin==4
   EPS=1e-8;
end


u=0; degu= -Inf;
v=a; degv=dega;

while degv-degb>=0
  d=degv-degb;
  c=-v(degv+1)/b(degb+1);
  [v,degv]=polred_e(v,degv,b,degb,c,d,EPS);
  [u,degu]=polred_e(u,degu,[1],0,-c,d,EPS);
end