感谢！我学习一下~

ivf自带有个erf函数，变换后可以用于求标准正态分布函数。

[Fortran] 纯文本查看 复制代码

program test
real(kind=8) a,x
x=1.0d20
a= 0.5*erf(sqrt(2.0)*x) + 0.5
print*,a

x=0.0d20
a= 0.5*erf(sqrt(2.0)*x) + 0.5
print*,a

x=-1.0d20
a= 0.5*erf(sqrt(2.0)*x) + 0.5
print*,a
end 

感谢版主指点，我尝试去做一下

这是一个简单的函数，没必要用到 MKL

但由于 Φ(x) 是一个超越函数，所以数值计算上只能用级数求和来达到。可以从 x 反向求和，一直到足够精度为止（不可能到负无穷）

http://fcode.cn/code_prof-33-1.html
这是你想要的吗？