xtrnlu.f

Go to the documentation of this file.

      SUBROUTINE xtrnlu(X,M,NB,UG)
*
*
*       Harmonic test potential.
*       ------------------------
*
      IMPLICIT REAL*8 (a-h,m,o-z)
      REAL*8 x(1),m(1)
*
*
      ug=0.0
      l=0
      DO i=1,nb
      DO k=1,3
      l=l+1
      ug=ug-5.0*x(l)**2*m(i)
      END DO
      END DO
*
      RETURN
      END