R kemungkinan normal

mungkin=function(sd1,n1,xbar1,value1,value2)
{
    es1 = sd1/(sqrt(n1))
    z1 = (value1-xbar1)/es1

    cat("Jumlah sample: ", n1); cat("   trunc 2: ", round(n1,2))
    cat("\n")

    cat("Population mean: ", xbar1); cat("   trunc 2: ", round(xbar1,2))
    cat("\n")

    cat("Population standard deviation: ", sd1); cat("   trunc 2: ", round(sd1,2))
    cat("\n")

    cat("Nilai batas kiri: ", value1); cat("   trunc 2: ", round(value1,2))
    cat("\n")

    cat("Nilai batas kanan: ", value2); cat("   trunc 2: ", round(value2,2))
    cat("\n");cat("\n")

    cat("Standard Error: ", es1); cat("   trunc 2: ", round(es1,2))
    cat("\n")
    cat("Nilai Z kritis kiri: ", z1); cat("   trunc 2: ", round(z1,2))
    cat("\n")
    z2 = (value2-xbar1)/es1
    cat("Nilai Z kritis kanan: ", z2); cat("   trunc 2: ", round(z2,2))
    cat("\n");cat("\n")
    kemungkinan=pnorm(z2)-pnorm(z1)
    cat("Nilai kemungkinan: ", kemungkinan); cat("   trunc 2: ", round(kemungkinan,2))
    cat("\n")
    cat("Nilai kemungkinan dalam persen: ", 100*kemungkinan);
    cat("   trunc 2: ", round(kemungkinan*100,2))
    cat("\n")
}