*     G11BCF Example Program Text
*     Mark 17 Release. NAG Copyright 1995.
*     .. Parameters ..
      INTEGER          NIN, NOUT
      PARAMETER        (NIN=5,NOUT=6)
      INTEGER          MMAX, LTMAX
      PARAMETER        (MMAX=5,LTMAX=54)
*     .. Local Scalars ..
      INTEGER          I, IFAIL, J, K, MAXST, MCELLS, MDIM, NCELLS,
     +                 NCOL, NDIM, NROW
      CHARACTER        STAT
*     .. Local Arrays ..
      DOUBLE PRECISION AUXT(LTMAX), STABLE(LTMAX), TABLE(LTMAX),
     +                 WK(LTMAX)
      INTEGER          IDIM(MMAX), ISDIM(MMAX), IWK(3*MMAX),
     +                 MLEVEL(MMAX)
*     .. External Subroutines ..
      EXTERNAL         G11BCF
*     .. Executable Statements ..
      WRITE (NOUT,*) 'G11BCF Example Program Results'
*     Skip heading in data file
      READ (NIN,*)
      READ (NIN,*) STAT, NCELLS, NDIM
      IF (NCELLS.LE.LTMAX .AND. NDIM.LT.MMAX) THEN
         READ (NIN,*) (TABLE(I),I=1,NCELLS)
         READ (NIN,*) (IDIM(J),J=1,NDIM)
         READ (NIN,*) (ISDIM(J),J=1,NDIM)
         MAXST = LTMAX
         IFAIL = 0
*
         CALL G11BCF(STAT,TABLE,NCELLS,NDIM,IDIM,ISDIM,STABLE,MAXST,
     +               MCELLS,MDIM,MLEVEL,AUXT,IWK,WK,IFAIL)
*
         WRITE (NOUT,*)
         WRITE (NOUT,*) ' Marginal Table'
         WRITE (NOUT,*)
         NCOL = MLEVEL(MDIM)
         NROW = MCELLS/NCOL
         K = 1
         DO 20 I = 1, NROW
            WRITE (NOUT,99999) (STABLE(J),J=K,K+NCOL-1)
            K = K + NCOL
   20    CONTINUE
      END IF
*
99999 FORMAT (10F8.2)
      END
