c--------------------------------------------------c
c     qmcb.f computes critical values              c
c       for 2-sided mcc (dunnett) inference        c
c     last revision: 03-12-92 by jch               c
c--------------------------------------------------c
      program main
      implicit doubleprecision (a-h,o-z)
c6    implicit real*8 (a-h,o-z)
c     real conf,fact1
      external glv,ghc1
      dimension fact1(100)
      dimension a(100),b(100),c(100)
      dimension ph(48),wh(48),whln(48),zh(48)
      dimension pl(48),wl(48)
      dimension ah(7),bh(12)
      common/list2/a,b,c
      common/list5/dinfnu
      common/list6/ph,wh,zh,whln
      common/list7/pl,wl
      common/list8/zero,zsr,zlog,zln
      common/list9/sr2,pi,srpi,pib2
      common/list10/iin,iout,iouti
      common/list12/ah,bh
c     
      write(6,1010)
 1010 format(' ','input number of treatments')
      read(5,*) k
      km1=k-1
      write(6,1020)
 1020 format(' ','input confidence coefficient')
      read(5,*) conf
      cc=dble(conf)
      write(6,1060)
 1060 format(' ','input array of factor pattern = sqrt(communality)')
      read(5,*) (fact1(j),j=1,km1)
      write(6,1070)
 1070 format(' ','input degrees of freedom')
      read(5,*) nu
      dnu=dble(float(nu))
c
      call qmccv(km1,fact1,1.0d0,nu,1,cc,d)
c
   30     write(6,230) k,conf,nu,d
  230     format(' ','k=',i2,1x,'cc=',f4.3,1x,'nu=',i3,1x,
     *               'd=',f5.3)
c
      stop
      end