# Hitung Z score & probability
    # Sunu Pradana
    # Samarinda, Indonesia, Agustus 2014
    #		
    # Creative Commons Attribution-ShareAlike 4.0
    # 
    # Usage (contoh penggunaan): 
    # mungkin(1, 0, 1, -2, 2)
    # mungkin(5, 200, 31, 187.4, 221.2)
    #
    # sd1 = Population standard deviation
    # xbar1 = Population mean
    # n1 = Jumlah sample (sample size)
    # value1 = left value
    # value2 = right value
    # Untuk nilai ekstrim (for extrme value), gunakan minimal: 
    # kiri (left)   = -1E99
    # kanan (right) =  1E99


    mungkin=function(sd1,xbar1,n1,value1,value2)
    {
        es1 = sd1/(sqrt(n1))
        z1 = (value1-xbar1)/es1
        z2 = (value2-xbar1)/es1
        kem=pnorm(z2)-pnorm(z1)
    
        cat("\n")
        cat( sprintf("%08s % 20.4f
       \r%08s % 20.4f
       \r%08s % 20.4g
       \r\n%08s % 20.4E
       \r%08s % 20.4E
       \n" ,
        "sd1:",sd1,"xbar1:",xbar1,"n1:",n1,"value1:",value1,"value2:",value2), fill = TRUE);
        
       #cat(" \n ");
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Jumlah sample",                   ":", n1,     "[trunc2:]", round(n1,2)),     fill = TRUE )
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Population mean",                 ":", xbar1,  "[trunc2:]", round(xbar1,2)),  fill = TRUE )
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Population standard deviation",   ":", sd1,    "[trunc2:]", round(sd1,2)),    fill = TRUE )    
        cat("\n")
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai batas kiri (left)",         ":", value1, "[trunc2:]", round(value1,2)), fill = TRUE )
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai batas kanan (right)",       ":", value2, "[trunc2:]", round(value2,2)), fill = TRUE )
        cat("\n")
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Standard Error",                  ":", es1,    "[trunc2:]", round(es1,2)),    fill = TRUE )
        cat("\n")
        cat( sprintf("%-31s %1s % 10.4g %10s % 5.2E ","Nilai Z kritis kiri (left)",      ":", z1,     "[trunc2:]", round(z1,2)),     fill = TRUE )
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2E ","Nilai Z kritis kanan (right)",    ":", z2,     "[trunc2:]", round(z2,2)),     fill = TRUE )
        cat("\n")
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Nilai kemungkinan (probability)", ":", kem,    "[trunc2:]", round(kem,2)),    fill = TRUE )
        cat( sprintf("%-31s %1s % 10.4E %10s % 5.2f ","Nilai kemungkinan dalam persen",  ":", kem*100,"[trunc2:]", round(kem*100,2)),fill = TRUE )
        cat("\n")
    }