標準正規分布において，標準化得点以上の値をとる確率を与える関数です。

# 正規分布の上側確率

function gxp(x,    is, i, c, p, y, z, pi2)
{
   pi2 = 0.398942280401432677940
   
   is = -1
   y = abs(x)
   c = y*y
   p = 0.0
   z = exp(-c*0.5)*pi2
   if (y < 2.5) {
       for (i = 20; i > 0; i--) {
           p = i*c/(i*2+1+is*p)
           is = -is
       }
       p = 0.5-z*y/(1.0-p)
   }
   else {
       for (i = 20; i > 0; i--) {
           p = i/(y+p)
       }
       p = z/(y+p)
   }
   return (x < 0.0) ? 1.0-p : p
}