SUBROUTINE POLY(X,PF,CONS,FA,MPPF,NCONS,NX,MP)
      IMPLICIT REAL*8 (A-H,O-Z)
      PARAMETER (IT=63)
      DIMENSION MPPF(IT+1,*)
      DIMENSION X(NX),PF(*),CONS(*)
      E=X(1)
      Z=X(2)
      ALE=LOG(E)
      ALZ=LOG(Z)
      PF(1)=1
      PF(2)=ALZ
      PF(3)=ALE
      
      PF(4)=ALZ*ALZ
      PF(5)=ALZ*ALE
      PF(6)=ALE*ALE
      
      PF(7)=ALZ*ALZ*ALZ
      PF(8)=ALZ*ALZ*ALE
      PF(9)=ALZ*ALE*ALE
      PF(10)=ALE*ALE*ALE
      
      PF(11)=PF(7)*ALZ
      PF(12)=PF(8)*ALZ
      PF(13)=PF(9)*ALZ
      PF(14)=PF(10)*ALZ
      PF(15)=PF(10)*ALE
      
      PF(16)=PF(11)*ALZ
      PF(17)=PF(12)*ALZ
      PF(18)=PF(13)*ALZ
      PF(19)=PF(14)*ALZ
      PF(20)=PF(15)*ALZ
      PF(21)=PF(15)*ALE
      
      FA=0
      DO I=1,NCONS
        FA=FA+CONS(I)*PF(I)
      ENDDO  
      FA=EXP(FA)
      DO I=1,NCONS
        PF(I)=PF(I)*FA
      ENDDO  
      RETURN
      END