다음과 같은 서브루틴을 사용 했는데요
Subroutine potential(A,B,C,D,n,ef,dgi,bi,lam,ga)
dimension A(3000),B(3000),C(3000),D(3000),ef(3000),dgi(3000)
Do 970 i=1,n
A(i)=-bi*lam*ef(i)
B(i)=-(2**(3.d0/2.d0))*bi*ga*lam*dgi(i)
C(i)=-bi*lam*(ef(i)-2*ga*dgi(i))
D(i)=-(2**(3.d0/2.d0))*bi*ga*lam*dgi(i)
970 co