anchor.diff:        do while (dist > cutoff)
         !do while (i < 10)
             !! Reflection first
  -          write(*,*) i
  +          !write(*,*) i
             minsave = min
  +          minval = value(min)
             xsave = pos(1, min)
             ysave = pos(2, min)
             call reflect(pos, min, dim, alpha)
  
positive.code: subroutine downhill
        implicit none
        integer, parameter         :: dim = 4
        integer :: min, max, i, minsave, maxsave
        real*8, dimension(2, dim)  :: pos
        real*8, dimension(dim - 1)     :: value
        real*8  :: x, y, alpha = 0.5, beta = 0.5, gamma = 0.5, dist, cutoff = 0.00001
        real*8  :: xsave, ysave
  
        interface
  
negative.code:               call expand(pos, max, dim, gamma)
                call findval(pos, value, dim)
                call minmax(value, min, max)
  
                !! Setting things back, if expansion is worse than 
                !! standard reflection
  