anchor.diff:  
         do i = 1,dim -1
             call random_number(pos(1,i))
  -          pos(1,i) = (pos(1,i) - 0.5)* 10
  +          pos(1,i) = (pos(1,i) - 0.5)
             call random_number(pos(2,i))
  -          pos(2,i) = (pos(2,i) - 0.5) * 10
  +          pos(2,i) = (pos(2,i) - 0.5)
         end do
         
   end subroutine
  
positive.code:       do while (dist > cutoff)
        !do while (i < 10)
            !! Reflection first
            write(*,*) i
            minsave = min
            xsave = pos(1, min)
            ysave = pos(2, min)
            call reflect(pos, min, dim, alpha)
  
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
  