      REAL FUNCTION SIGMA(E)
C     Dissociative Collisions
C     H_2^++H_2->H^+_Proj
C     CROSS SECTION(CM^2)  :SIGMA
C     COLLISION ENERGY(KEV):E
C     JAERI-Data/Code 96-024(1996)
      REAL A(14)
      DATA EMIN1,EMAX1 /1.6E-3,5.0E-3/
      DATA EMIN2,EMAX2 /1.5E+0,1.0E+4/
      DATA A /4.676E+00,1.436E+00,4.274E+01,6.969E-01,2.404E+02,
     >        3.199E+00,1.272E+00,4.972E-02,3.313E-02,3.643E+01,
     >        5.423E-01,4.237E-05,5.873E+00,2.146E-01/
      DATA SIGMA0,ER /1.E-16,50.0/
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      F3(X)=SIGMA0*A(1)*(X/ER)**A(2)/(1.+(X/A(3))**(2.*A(2)))
      F4(X)=SIGMA0*A(1)*(X/ER)**A(2)/(1.+(X/A(3))**(A(2)+A(4)))
      F6(X,ALPHA)=SIGMA0*A(1)*(X/ER)**A(2)/
     >            (1.+(X/A(3))**(A(2)+A(4))+(X/A(5))**(A(2)+ALPHA))
      F11(X,ALPHA,BETA)=SIGMA0*A(1)*(X/ER)**ALPHA/
     >            (1.+(X/A(3))**(ALPHA+BETA)+(X/A(5))**(ALPHA+A(6)))
      F12(X,ALPHA,BETA)=F11(X,ALPHA,BETA)
      F13(X,ALPHA,BETA)=SIGMA0*A(1)*(X/ER)**ALPHA/
     >                  (1.+(X/A(3))**(ALPHA+BETA))
C
      SIGMA=0.0
      IF(E.GE.EMIN1 .AND. E.LE.EMAX1 .OR. 
     >   E.GE.EMIN2 .AND. E.LE.EMAX2) THEN
        SIGMA=F12(E,A(2),A(4))+A(7)*F12(E/A(8),A(2),A(4))+
     >        A(9)*F12(E/A(10),A(2),A(4))+
     >        A(11)*F12(E/A(12),A(13),A(14))
      END IF
 9999 RETURN
      END
