INORMALX.OMS
Script File:
# Description:
# Demonstrates that inormal is the inverse of cnormal.
#
clear
# grid of values between -1 and +1
p = seq(20) / 10d0 - 1.
# format for print out
format double "f25.7"
# evaluate normal distribution on the grid
cp = cnormal(p)
# evaluate inverse on cp
icp = inormal(cp)
# print heading with fields right adjusted
print align("p;inormal(cnormal(p))", ";", [25, 26], "right")
write("screen", [p, icp])
Output:
p inormal(cnormal(p))
-0.9000000 -0.9000000
-0.8000000 -0.8000000
-0.7000000 -0.7000000
-0.6000000 -0.6000000
-0.5000000 -0.5000000
-0.4000000 -0.4000000
-0.3000000 -0.3000000
-0.2000000 -0.2000000
-0.1000000 -0.1000000
0.0000000 -0.0000000
0.1000000 0.1000000
0.2000000 0.2000000
0.3000000 0.3000000
0.4000000 0.4000000
0.5000000 0.5000000
0.6000000 0.6000000
0.7000000 0.7000000
0.8000000 0.8000000
0.9000000 0.9000000
1.0000000 1.0000000
