Question: Inverse CDF computations failing or returning complex values

with(Statistics):
Phiinverse := (x) -> Quantile(Normal(0, 1), x) :

evalf(Phiinverse(1-1/2^24));

and the result is 5.691523310-2.147372008*I

Now, I can get a real-valued solution to this by using the equivalent
functionality from stats instead of Statistics.

with(stats):
PhiinvClassic:= (x) -> statevalf[icdf, normald[0, 1]](x):

evalf(PhiinvClassic(1-1/2^24));

Please Wait...