anchor.diff:  c Renamed to use alpha and beta arguments for compatibility with 
   c random.beta.
   
  -cf2py double precision dimension(nx),intent(in) :: x
  -cf2py double precision dimension(na),intent(in) :: alpha
  -cf2py double precision dimension(nb),intent(in) :: beta
  +cf2py real dimension(nx),intent(in) :: x
  +cf2py real dimension(na),intent(in) :: alpha
  +cf2py real dimension(nb),intent(in) :: beta
   cf2py integer intent(hide),depend(x) :: nx=len(x)
   cf2py integer intent(hide),depend(alpha),check(na==1 || na==len(x)) :: na=len(alpha)
   cf2py integer intent(hide),depend(beta),check(nb==1 || nb==len(x)) :: nb=len(beta)
  -cf2py double precision intent(out) :: like
  +cf2py real intent(out) :: like
         IMPLICIT NONE
         INTEGER i,nx,na,nb
  -      DOUBLE PRECISION like
  -      DOUBLE PRECISION x(nx),alpha(na),beta(nb), atmp, btmp
  -      DOUBLE PRECISION gammln
  +      REAL like
  +      REAL x(nx),alpha(na),beta(nb), atmp, btmp
  +      REAL gammln
   
         atmp = alpha(1)
         btmp = beta(1)
  
positive.code: c Combinations of x red balls
          if (nd .NE. 1) draws_tmp = draws(i)
          if (ns .NE. 1) s_tmp = success(i)
          if (nt .NE. 1) t_tmp = total(i) 
          like = like + combinationln(t_tmp-s_tmp, x(i))
          like = like + combinationln(s_tmp,draws_tmp-x(i))
          like = like - combinationln(t_tmp, draws_tmp)
        enddo
        return
        END
  
negative.code: c Return the uniform likelihood of x.
  c CREATED 12/06 DH
  
  cf2py real dimension(n), intent(in) :: x
  cf2py real dimension(nlower), intent(in) :: lower
  cf2py real dimension(nupper), intent(in) :: upper 
  cf2py integer intent(hide), depend(x) :: n=len(x)
  cf2py integer intent(hide), depend(lower) :: nlower=len(lower)
  cf2py integer intent(hide), depend(upper) :: nupper=len(upper)
  cf2py real intent(out) :: like
  
          IMPLICIT NONE
          
          INTEGER n, nlower, nupper, i
          REAL x(n), lower(nlower), upper(nupper)
          REAL like, low, high
                  
          low = lower(1)
          high = upper(1)       
  