Cancelled effects of previous commit due to several numerical troubles
Updated Psi function. In the general case, the r and r.int.0 functions should be evaluated on samples.
Using mapply, this becomes possible.
Vectorize takes more time because of several checks, internal objects build
This should be applied to Phi
Applied mapply to the evaluation of integrals in Phi