anchor.diff:  implicit none
   real*8, parameter :: inf = 1000000000.d0
   contains
  +
  +function isElIn(item,list) result(is)
  +  integer, intent(in) :: list(:), item
  +  integer :: i
  +  logical :: is
  +  is=.false.
  +  do i=1,size(list)
  +    if (list(i) == item) then
  +      is=.true.
  +    endif
  +  end do
  +
  +end function isElIn
   
   subroutine shuffle(a)
     ! http://rosettacode.org/wiki/Knuth_shuffle#Fortran
  
positive.code:   integer :: len, start,finish , i, temp
    len = size(item,1)
    call pick_range(2,len,start,finish)
    print *, start, finish
    print "(10I2)", (i, i=1,10)
    ! Inverse given range of genes
    do i=start,(finish+start-1)/2
      temp = item(i)
  
negative.code:   genNo = size(population,1)
    do i=1,popNo
      population(:,i) = (/ (j, j=1,genNo) /)
      call shuffle(population(2:genNo,i))
    enddo!i
  end subroutine growPopulation
  
  