In reply to Rick Wicklin.

Thanks a lot for your immediate response. I prefer to do the integration analytically, because it usually makes the computation much faster and more precise. It seems that the partial derivative of the CDF w.r.t. rho is equal to the PDF: