function y=rho0mu(z,mu)
if mu==1
  y=-log(1-1/z);
else
  y=(1/(mu-1))*z^(1-mu)*((1-1/z)^(1-mu)-1);
end
