正態分佈函式值

《概率論及其應用（卷1•第3版） 》第 7 章提到正態分佈函式：

隨後在表 7.1中對於 0.00 ≤ x ≤ 3.29 給出了正態分佈函式值：

我使用以下 PARI/GP 程式來生成這張表格中的 330 個值：

f(x)=0.5+intnum(y=0,x,exp(-0.5*y^2))/sqrt(2*Pi)
x=0;while(x<3.3,printf("%.2f:%.4f\n",x,f(x));x+=0.01)

發現其中一些值是錯誤的，如下所示：

x: 正確值 表格中的值 差
0.04: 0.5160 0.5159 1
0.68: 0.7517 0.7518 1
0.71: 0.7611 0.7612 1
0.88: 0.8106 0.8016 90
0.99: 0.8389 0.8380 9
1.13: 0.8708 0.8718 10
1.33: 0.9082 0.9083 1
1.34: 0.9099 9.9099 90000
1.58: 0.9429 0.9430 1
1.63: 0.9484 0.9485 1
1.75: 0.9599 0.9509 90
1.97: 0.9756 0.9758 2
1.98: 0.9761 0.9762 1
2.00: 0.9772 0.9773 1
2.21: 0.9864 0.9865 1
2.86: 0.9979 0.9989 10
2.87: 0.9979 0.9980 1
2.92: 0.9982 0.9983 1
2.93: 0.9983 0.9984 1
3.00: 0.9987 0.9986 1
3.05: 0.9989 0.9988 1
3.08: 0.9990 0.9989 1
3.22: 0.9994 0.9993 1
3.27: 0.9995 0.9994 1