corp.f

来自「Fortran科学计算程序,计算原子间力的大小,能级分布」· F 代码 · 共 12 行

**==corp.spg  processed by SPAG 4.52O  at 10:32 on 18 Oct 1996
      FUNCTION CORP(R, Rho)
      IMPLICIT NONE
      INCLUDE 'potential.inc'
      DOUBLE PRECISION sig3, ri3, R, CORP, Rho
 
      sig3 = SIG2*SQRT(SIG2)
      ri3 = sig3/(R*R*R)
      CORP = 4*PI*EPS4*(Rho**2)*sig3*(2*ri3*ri3*ri3/9-ri3/3)
      RETURN
      END