C
C See eq. 9.18 of Kris Heyde, "The Nuclear Shell Model", Springer-Verlag,
C 1994.
C
       FUNCTION VN0105(N1,L1,LAMB,N2,L2)
C      VN0105(N1,L1,LAMB,N2,L2) =(N1,L1/RHO**LAMB/N2,L2)
C        WHERE /N,L) IS THE RADIAL PART OF THE WAVE FUNCTION
C           ASSOCIATED WITH THE SYMMETRICAL HARMONICAL
C          OSCILLATOR POTENTIAL AND
C        WHERE RHO IS THE REDUCED R-COORDINATE
 1    NUL=0
      IF(N2-N1)2,3,3
2      NA=N1
      N1=N2
      N2=NA
      NA=L1
      L1=L2
      L2=NA
 3    S=0.
       NN1=2*N1
       NN2=2*N2
       DO 9 K=1,N1,1
      KK=2*K
       T1=VN0107(NN1)/VN0107(KK)
      M1=2*(N1-K+1)
      M2=2*(N2-N1+K)
       T2=VN0107(NN2)/(VN0107(M1)*VN0107(M2))
       M3=2*(N1-K+1)+L1+L2+LAMB+1
       T3=VN0107(M3)
      T4=1.
      T5=1.
      IF(K-1)31,5,31
31     AL=L1-L2-LAMB
      KA=K-2
       DO 4 I=NUL,KA,1
      T4=T4*(AL/2.+FLOAT(I))
4      CONTINUE
5      IF((K.EQ.1).AND.(N1.EQ.N2))GOTO 7
      BL=L2-L1-LAMB
      KB=K-2+N2-N1
       DO 6 J=NUL,KB,1
      T5=T5*(BL/2.+FLOAT(J))
6      CONTINUE
7      S=S+(T1*T2*T3*T4*T5)*10.**(FLOAT(LAMB)/2.
     1+FLOAT(-2*K+N1-N2+2))
 9    CONTINUE
11     N5=2*(N1+L1)+1
       N6=2*(N2+L2)+1
10     SQ=VN0107(NN1)*VN0107(NN2)*VN0107(N5)*VN0107(N6)
       SQ1=SQRT(SQ)
12     VN0105=(S*FLOAT((-1)**(N1+N2)))/SQ1
      RETURN
       END
      FUNCTION VN0107(N)
C     VN0107(N)=GAMMA(N/2)/(1O**(N/2-1))
C        WHERE GAMMA IS THE WELL KNOWN GAMMA-FUNCTION AND
C        WHERE N IS AN INTEGER
C     FCT(M)=(M-1)!/(10**(M-1))
C       WHERE M IS A POSITIVE INTEGER
      COMMON/VN06FC/FCT(40)
10     I=0
  1   K=N/2
      IF(K*2-N)3,2,12
  2   IF(N)15,15,21
 21   GA=FCT(K)
      IF(I)90,90,13
 3    K=N/2
      IF(K)31,6,31
31     K1=2*K
       GA=(FCT(K1)*SQRT(31.4159))/(FCT(K)*2.**(2*K-1))
      IF(I)90,90,13
  6   GA=SQRT(31.4159)
       IF(I)90,90,13
 12   I=1
      N=N*(-1)+2
       GOTO 3
13     G=GA
       N=N*(-1)+2
       FA=(-1)**((1-N)/2)
       GA=(FA*31.4159)/G
      GOTO 90
90    VN0107=GA
15     RETURN
       END