分子模拟课作业

RIJSQ .LT. RCUTSQ ) THENР SR2 = SIGSQ / RIJSQР SR6 = SR2 * SR2 * SR2Р UIJ = SR6 * ( SR6 - 1.0 )Р WIJ = SR6 * ( SR6 - 0.5 )Р U = U + UIJР W = W + WIJР ENDIFР ENDIFР100 CONTINUEР U = 4.0 * UР W = 48.0 * W / 3.0Р RETURNР END SUBROUTINE ENERGYРc *****************************************************************РC **第三个子函数——该函数用于生成0~1的随机数**Р REAL FUNCTION RANF ( DUMMY )РC *******************************************************************Р INTEGER L, C, MР PARAMETER ( L = 1029, C = 221591, M = 1048576 ) Р INTEGER SEEDР REAL*8 DUMMYР SAVE SEEDР DATA SEED / 0 /