Fortran77怎么编译高斯误差函数erf(x)? 50

最近在利用fortran77写abaqus子程序时，碰到要用高斯误差函数erf(x),但fortran77内部函数库里没有这个函数，有人碰到过类似的问题吗？需要怎么解决？我在网上找了高斯误差函数的源代码，编译有错，目前还没走通。。。