%RHO01_MU
%
function y=rho01(z,mu)
f0='%8.0f %4.0f %20.12e\n'; 
%for Examples 1 and 2
y=(-1)^(mu-1)*rho0mu(z,mu);
%for Example 3
%global imod m dig
%digits(dig)
%if imod==1
%  y=(-1)^(mu-1)*rho0mu(z,mu);
%else
%  c=parfrac(dig,m,z);
%  cmax=max(subs(abs(c(mu,:))))
%  sn=vpa(0);
%  for nu=1:mu
%    sn=sn+c(mu,nu)*rho0mu(conj(z),nu);
%  end
%  abssn=abs(subs(sn))
%  sk=vpa(0);
%  for nu=mu+1:mu+m
%    sk=sk+c(mu,nu)*rho0mu(z,nu-mu);
%  end
%  abssk=abs(subs(sk)) 
%  y=(-1)^(m+mu-1)*(sn+sk);
%end
%for Example 4
%z0=1/sqrt(1-z);
%if mu==1
%  y=z0*log((1-z0)/(1+z0));
%else
%  R=(1/2)*log((1+z0)/(1-z0));
%  for k=1:mu-1
%    R=(1-1/(2*k))*R+z0/(2*k*(1-z0^2)^k);
%  end
%  y=-2*z0^(2*mu-1)*R;
%end
%for Examples 5 and 6
%global imod m dig
%digits(dig);
%s=symengine;
%if imod==1
%  y=(-1)^mu*z^(1-mu)*exp(-z)*vpa(feval(s,'Ei',vpa(mu),-z));
%else
%  c=parfrac(dig,m,z);
%  cmax=max(subs(abs(c(mu,:))))
%  sn=vpa(0);
%  for nu=1:mu
%    sn=sn+c(mu,nu)*(conj(z))^(1-nu)*vpa(feval(s,'Ei',vpa(nu),-conj(z)));
%  end
%  abssn=abs(subs(sn))
%  sk=vpa(0);
%  for nu=mu+1:mu+m
%    sk=sk+c(mu,nu)*z^(1-nu+mu)*vpa(feval(s,'Ei',vpa(nu-mu),-z));
%  end
%  abssk=abs(subs(sk))
%  y=(-1)^(mu+m)*(exp(-conj(z))*sn+exp(-z)*sk);
%end
