axis([0,2,0,1])
hold on
xx=[0:0.01:2];
for mu=0.1:0.1:0.9
  yy=[];
  for x=xx
    yy=[yy Rho(x,mu)];
  end
  plot(xx,yy,'b')
end