標準正規分布において，上側確率が p に対するパーセント点を与える関数です。

# 正規分布のパーセント点を求める。

function pxg(p,     y, a, b, c, d, e, f, g, h, i, j, k, l, m, n, o)
{
   a =  0.161614603871299e-15
   b = -0.178837229464158e-13
   c =  0.888943322341436e-12
   d = -0.261195052190738e-10
   e =  0.500533655324426e-09
   f = -0.646627744010232e-08
   g =  0.553751002988654e-07
   h = -0.283629174340630e-06
   i =  0.539877747512783e-06
   j =  0.122127251950965e-05
   k =  0.154574365345128e-04
   l = -0.234768842092152e-03
   m = -0.830686207721635e-03
   n =  0.370684581838722e-01
   o =  1.57079636931376
   y = -log(4*p*(1-p))
   return sqrt(((((((((((((((y*a+b)*y+c)*y+d)*y+e)*y+f)*y+g)*y+h)*y+i)*y \
          +j)*y+k)*y+l)*y+m)*y+n)*y+o)*y)
}