fortran中有没有贝塞尔函数可以调用啊？？？

fortran中有没有贝塞尔函数可以调用啊？？？ help中搜到的以下内容是啥意思，大侠们帮帮忙啊，我是初学者。
Module: USE DFPORT
Syntax result = DBESJ0 (posvalu)
result = DBESJ1 (posvalu)
result = DBESJN (n, posvalu)
result = DBESY0 (posvalu)
result = DBESY1 (posvalu)
result = DBESYN (n, posvalu)

posvalue
(Input) REAL(8). Independent variable for a Bessel function. Must be greater than or equal to zero.

n
(Input) Integer. Specifies the order of the selected Bessel function computation.

如果没有，有没有大侠们编好的计算贝塞尔函数的源程序呢？

双精度计算第一类零阶贝塞尔函数：

FUNCTION bessj0(x)
REAL bessj0,x
REAL ax,xx,z
DOUBLE PRECISION p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,&
r1,r2,r3,r4,r5,r6,s1,s2,s3,s4,s5,s6,y
SAVE p1,p2,p3,p4,p5,q1,q2,q3,q4,q5,r1,r2,r3,r4,r5,r6,&
s1,s2,s3,s4,s5,s6
DATA p1,p2,p3,p4,p5/1.d0,-.1098628627d-2,.2734510407d-4,&
-.2073370639d-5,.2093887211d-6/, q1,q2,q3,q4,q5&
/-.1562499995d-1,.1430488765d-3,-.6911147651d-5,&
.7621095161d-6,-.934945152d-7/
DATA r1,r2,r3,r4,r5,r6/57568490574.d0,-13362590354.d0,&
651619640.7d0,-11214424.18d0,77392.33017d0,&
-184.9052456d0/,s1,s2,s3,s4,s5,s6/57568490411.d0,&
1029532985.d0,9494680.718d0,59272.64853d0,&
267.8532712d0,1.d0/
if(abs(x)<8.) then
y=x**2
bessj0=(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6)))))/&
(s1+y*(s2+y*(s3+y*(s4+y*(s5+y*s6)))))
else
ax=abs(x)
z=8./ax
y=z**2
xx=ax-.785398164
bessj0=sqrt(.636619772/ax)*(cos(xx)*(p1+y*(p2+y*&
(p3+y*(p4+y*p5))))-z*sin(xx)*(q1+y*(q2+y*&
(q3+y*(q4+y*q5)))))
endif
END FUNCTION bessj0

besj0, besj1, besjn, besy0, besy1, besyn 这些都是。具体请看编译器说明书（帮助文档）

你找到的帮助信息就可以帮助你了。你没有必要再百度提问。 仔细看你的帮助信息，按照它的样子使用贝塞尔函数就可以。