SUBROUTINE FMMM (QF2,F1C,F2C,PSP,PD)
C *** SEE R H LANDAU, PROGRAM LPOTT, 1981
C     FMMM INTERPOLATES MALCOLM MCMILLAN FORM FACTORS AND FORMS F1C AND
C     F2F FROM THEM (BASED ON APPROX TRITON WAVE FUNCTNS)
C     PSP = S PRIME PROB     PD= DSTATE PROB     QF2= Q2 IN FERM-2
      IMPLICIT  REAL*8   (A-H,O-Z)
      DIMENSION Q2(18), F3456(4,18), Y(4), Q20(27), F0(1,27), Y0(1)
      COMMON /MMM/ F3456
C     N.B. WE MUST HAVE F3456 SAME DIMENSION IN BLOCK DATA AS HERE FOR T
C
      DATA Q2/0.0,0.1,1.1,2.1,3.1,4.1,5.1,6.1,7.1,8.1,9.1,10.1,11.1,12.
     11 ,13.1,14.1,15.1,16.1/
      DATA Q20
     1     /0.,0.3,0.9,1.5,2.1,2.7,3.3,3.9,4.5,5.1,5.7,6.3,6.9,7.5,8.0,
     2  9.0,10.0,11.,12.,13.,14.,15.,16.,17.,18.,19.,20./,F0/
     31.0,0.88124,0.69363,0.55333,0.44569,0.36153,0.29475,0.24114,
     40.19770,0.16221,0.13305,0.10896,0.08897,0.07232,0.605673E-1,
     50.417046E-1,0.276478E-1,0.171613E-1,0.935277E-2,0.35677E-2,
     6-0.679445E-3,-0.375343E-2,-0.592979E-2,-.741903E-2,-0.838265E-2,
     7-0.894523E-2,-0.920333E-2/
      IF (QF2.LE.16.1) GO TO 20
      Q = SQRT(QF2)
      IF ((0.4111*QF2.LT.150.).AND.(1.91*(Q-4.4721)**2.LT.150.))
     1 GO TO 10
C     Q2 TOO LARGE FOR EXP, SET=0
      F1C = 0.
      F2C = 0.
      RETURN
C     Q2 OUTSIDE OF TABLE, USE ANNALTIC FORMS FITTED TO LAST OINTS
 10   Y0(1) = EXP(-0.4111*QF2)-0.94722E-2*EXP(-1.906*(Q-4.4721)**2)
      Y0(1) = Y0(1)*1.23469
      Y(1) = 0.589748*EXP(-0.124043*QF2)
      Y(2) = 0.325014*EXP(-0.125324*QF2)
      Y(3) = (1./19.0048)*(EXP(-0.4667*QF2)-0.25369E-01*EXP(-1.94305*(Q-
     14.1833)**2))
      Y(4) = -0.831500*EXP(-0.587182*QF2)
      GO TO 30
 20   CONTINUE
C     INTERPOLATE ON FORM FACTORS,,N.B. MAY HAVE TO CAHGE THIS FOR MORE
C     ACCURACY AND INTERPOLATE ON LOGS OF FF AND THE EXPEONENTIATE
      CALL LAGRNG (QF2,Q2,Y,F3456,18,4,2,18,4)
      CALL LAGRNG (QF2,Q20,Y0,F0,27,1,2,27,1)
 30   PS = 1.-PSP-PD
      F1C = PS*Y0(1)+0.5*PSP*Y(2)+19.0048*PD*Y(3)
      F2C = SQRT(PSP*PS)*Y(1)+19.0048*PD*Y(4)
C     WRITE(6,900) QF2,Y0(1),Y(1),Y(2),Y(3),Y(4),F1C,F2C
C     900 FORMAT(" ", F10.3,7E15.5)
      RETURN
      END