anchor.diff:  end subroutine
   
   !! Finds the minima and maxima
  +!! We need a better way to search through the initial elements in our array!
   subroutine minmax(value, min, max)
         implicit none
         real*8, intent(in)   :: value(:)
  
positive.code: !! minimum currently set to 0.5, should be found via Downhill Simplex!
  subroutine findval(pos, value, dim)
        implicit none
        real*8, intent(in), dimension(:,:):: pos
        real*8, intent(out), dimension(:) :: value
        real*8              :: sourcex = 0.5, sourcey = 0.5
        integer, intent(in) :: dim
        integer             :: i
  
        do i = 1, dim - 2
            value = sqrt((pos(1,i) - sourcex) * (pos(1,i) - sourcex) &
                    + (pos(2,i) - sourcey) * (pos(2,i) - sourcey))
        end do
  
        do i = dim - 2, dim
            value = 0
        end do
  
  end subroutine
  
negative.code:           call random_number(pos(1,i))
            call random_number(pos(2,i))
        end do
        
  end subroutine
  
  !! finds centroid position 
  pure subroutine centroid(pos, min, max, dim)
        implicit none
        real*8, intent(inout)  :: pos(:,:)
  