c main program for minimization of function in cumom.f
c
      parameter (n = 2, lwork = n*(n+10))
      double precision  param(n), x(n), f, work(lwork), pa(35), 
     c   pb(35), alpha(3)
c
c initial parameter estimates
c
      param(1) = 0.25
      param(2) = 0.25
c
      do 10 i=1,1
        call cumcmom(param,f)
        write(*,*) f
 10   continue
c
C
C MINIMIZE FUNCTION
C
c     CALL UNCMND (N, X0, CALCF, X, F, IERROR, WORK, LWORK)
C
C PRINT RESULTS
C
c      WRITE (*,*) 'UNCMND RESULTS'
c      IF (IERROR .NE. 0) WRITE (*,*) ' ERROR CODE =', IERROR
c      WRITE (*,*) ' F(X*) =', F
c      WRITE (*,*) ' X* ='
c      WRITE (*,800) (X(I), I = 1,N)
C
C
      STOP
800   FORMAT (5F12.6)
      END