e2 <- function(x) {
    out <- 0*x
    for(i in 1:length(x))
    out[i] <-integrate(function(y) qnorm(y),lower=0,upper=x[i])$value
    out }
integrate(e2,lower=0, upper=a)$value


above is my code , when a is small , say a<0.45  the result is right .
however , when a>0.5
the result is incorrect .  why ?    thank you .

	[[alternative HTML version deleted]]


