%MOD_D  Modification by the same divisor
%
function ab=mod_d(N,m,z,ab0,rho0_mu)
%
% The mu-th component of the vector rho0_mu is
% the input parameter \rho_0^{[\mu]} to Algorithm 1
% in the mu-th modification by the divisor t-z
%
global imod mu 
ab=ab0;
for mu=1:m
  ab0=ab;
  r0=feval(rho0_mu,z,mu);
  if imod==1
    ab=schri4(N,ab0,z,r0);
  else
    ab=schri4(N,ab0,conj(z),r0);
  end
end
