text:       PROGRAM xks2d2s
  C     driver for routine ks2d2s
        INTEGER NMAX
        PARAMETER(NMAX=1000)
        INTEGER idum,j,jtrial,n1,n2,ntrial
        REAL d,gasdev,prob,shrink,u,v,x1(NMAX),y1(NMAX),x2(NMAX),y2(NMAX)
  1     write(*,*) 'INPUT N1,N2'
        read(*,*,END=999) n1,n2
        if (n1.gt.NMAX) then
          write(*,*) 'n1 too large.'
          goto 1
        endif
        if (n2.gt.NMAX) then
          write(*,*) 'n2 too large.'
          goto 1
        endif
        write(*,*) 'WHAT SHRINKAGE?'
        read(*,*,END=999) shrink
        write(*,*) 'HOW MANY TRIALS?'
        read(*,*,END=999) ntrial
        if (ntrial.gt.NMAX) then
          write(*,*) 'Too many trials.'
          goto 1
        endif
        idum=-287-ntrial-n1-n2
        do 13 jtrial=1,ntrial
          do 11 j=1,n1
            u=gasdev(idum)
            v=gasdev(idum)*shrink
            x1(j)=u+v
            y1(j)=u-v
  11      continue
          do 12 j=1,n2
            u=gasdev(idum)*shrink
            v=gasdev(idum)
            x2(j)=u+v
            y2(j)=u-v
  12      continue
          call ks2d2s(x1,y1,n1,x2,y2,n2,d,prob)
          write(*,'(1x,a7,2f12.6)') 'D,PROB= ',d,prob
  13    continue
        goto 1
  999   write(*,*) 'NORMAL COMPLETION'
        STOP
        END
  