anchor.diff:       return
     end if
   
  +  call initialize_dirichlet(counts, nc)
     call compute_multiplicities(counts, nc)
   
     estimate = hdir(alpha)
  
positive.code:      return
    end if
  
    call compute_multiplicities(counts, nc)
  
    call compute_integration_range()
  
negative.code: 
  subroutine nsb(n,counts,nc,estimate,err_estimate)
    use iso_fortran_env
    use dirichlet_mod, only: initialize_dirichlet, compute_multiplicities, dirichlet_finalize
    use nsb_mod, only: hnsb
    use nsb_mod, only: compute_integration_range
    implicit none
  