请教一段Fortran代码

cuijxin

subroutine F(X,Y,D,N0)
real(8)::Y,D,X,A,GAMMA,B
integer::N0
A=N0/2.
B=2./N0
if(N0<=2) then
      D=-2**A*GAMMA(A)*EXP(0.5*Y)**(1-A)
else
      D=-A*2**A*GAMMA(A)*EXP(0.5*Y**B)
end if
X=X
end subroutine F

这个函数里的D=-2**A*GAMMA(A)*EXP(0.5*Y)**(1-A)是什么意思啊？

mofangbao

你应该是不明白那个GAMMA函数是啥意思吧，其实百度一下就有了，γ函数，好像gfortran自带的，那个real后面的GAMMA定义是你自己加的？这个函数的具体介绍看这里：
http://zh.wikipedia.org/zh/%CE%93%E5%87%BD%E6%95%B0
上面表达式中的其他东西应该很简单的吧，不懂的自己翻fortran教材了