qnorm(matProbabilit​ies,dblMean,matSigm​as,boolUseApproxima​tion)

This function uses the normal inverse distribution function to compute the number of standard deviations that probability values are away from the mean, assuming a Gaussian distribution; i.e., it transforms known probabilities to standard deviation units (see also: http://en.wikipedia.org/wiki/Probit )
Syntax: matSigmaOut = qnorm(matProbabilities,dblMean,matSigmas,boolUseApproximation)
By Jorrit Montijn, 08-01-15 (University of Amsterdam), based on original code by Anders Holtsberg, 13-05-94. Improvements include:
- Approximation of inverse error function using scaled logarithmic function for probability values lower than "eps" (~8sd)
- Corrected potential input check failure for matrices with more than 2 dimensions
- Comments added and code made more readable

N.B.: for output values of >~8 sd's the returned values are approximations and are not directly computed using erfinv(). The logarithmic approximator function is relatively precise, but caution is advised. If you do not wish to use this approximation, supply a fourth input argument set to 'false'; the function will then return 'Inf' where the direct computation fails

Any suggestions on improvements on the approximator function of the inverse error function are welcome