anchor.diff:        if(kmi.lt.1) kmi=1
         kma=idint(qmax*dfloat(nbins))+1
         kkmi=idint(rn*qmin)
  -      if(kkmi.lt.minbin) kkmi=minbin
  +      if(kkmi.lt.minbinpts) kkmi=minbinpts
         bpow=0.0d0
   c
   c     The following variables are defined for the extension
  
positive.code: c
  c
        subroutine eebls(n,t,x,u,v,nf,fmin,df,nbins,qmin,qmax,
       &     p,bper,bpow,depth,qtran,in1,in2)
  c
  c------------------------------------------------------------------------
  
negative.code: c     -- Dimensions of arrays {y(i)} and {ibi(i)} MUST be greater than
  c     or equal to  *nbinsmax*.
  c     -- The lowest number of points allowed in a single bin is equal
  c     to   MAX(minbinpts,qmin*N),  where   *qmin*  is the minimum transit
  c     length/trial period,   *N*  is the total number of data points,
  c     *minbinpts*  is the preset minimum number of the data points per
  c     bin.
  c
  c========================================================================
  