about hankel function

xj2070 · 发表于 2005-8-28 16:56

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

1由于国内常见的fortran程序集合中没有直接关于hankel 函数的，且IMSL 函数中也只有 bessel 函数，希望下面的小程序对大家有一点帮助。 2。是基于bessel 的IMSL 函数的 3。subroutine hankel(order,class,x,res) ! Function: Compute the hankel function values ! ! Note: ! ! This is a very useful subroutine,for that in IMSL ! the hankel function is not provided directly. ! The bessel function in IMsl LIBARY is not very good ! when CLASS=3 WHICH REPRENTS THE FIRST KIND BESSEL FUNCTION ! when class=4 WHICH REPRENTS THE second KIND BESSEL FUNCTION ! ! Input: ! ! order-------- the order number describing the hankel's function ! class-------- the hankel fuction's classfication ! x ---------the variable (real) ! ! Output: ! res --------- the results (complex) ! ! Implicit None Integer order,N !the order number requried for computation Integer class !the kind number Real x,CBS1(order+1),CBS2(order+1) !argument x complex res,i 
!--------------------------------------------------------------------
 i=(0,1) ! the false root N=order+1 ! in order compute the order-th results
 Call BSJNS(x,N,CBS1) Call BSYS(0.0,x,N,CBS2) 
!--------------------------------------------------------------------
 If(class==1)Then
 res=CBS1(N)+i*(CBS2(N)) !the first kind hankel function Else if(class==2) then res=CBS1(N)-i*(CBS2(N)) !the second kind hankel function Else if(class==3) then !the REPRENTS THE FIRST KIND BESSEL FUNCTION res=CBS1(N)+i*0.0
 Else if(class==4) then ! the REPRENTS THE Second KIND BESSEL FUNCTION res=CBS2(N)+i*0.0 Else write(*,*) 'the order num is illegal' stop End if End subroutine

christy · 发表于 2005-8-28 20:16

谢谢搂主，希望搂住以后能提供更多的好东西

账号		自动登录	找回密码
密码			我要加入

[Fortran] about hankel function

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（xj2070）about hankel function