anchor.diff:      end do
   
     end subroutine
  -
   !the fourth phase of verlet integration. Calculating the final velocity for the next phase
     subroutine verlet_integration_fourth(object, siz, list, timescale)
       implicit none
  
positive.code:       do i = 1, siz
          call total_acceleration(i, siz, list)
        end do
  
        !Fourth step for all objects
        do i = 1, siz
          call verlet_integration_fourth(i, siz, list, timescale)
  
negative.code:     close(unit=1, status='keep')
  
    end subroutine
  
  !Subroutine for one iteration of verlet_integration for one object the first phase. siz is the amount of objects in list
    subroutine verlet_integration_first(object, siz, list, timescale)
  