Foldnorminv-func
x = foldnorminv(p,mu,sigma) computes the deviate x, associated with the given lower tail probability, p, of the folded normal distribution, with distribution parameters mu and sigma.