anchor.diff:     real(kind = xp) :: E         ! Total energy
      integer :: iter = 0          ! Number of iterations ran
      integer :: dump_freq = 10    ! Frequency at which the system is sampled
  -   integer :: maxiter = 10000   ! Maximum number of iterations
  +   integer :: maxiter = 3       ! Maximum number of iterations
   
  +   real :: start_time, total_time
  +
  +   !---------------------------------------------
  +   ! Start time
  +   !---------------------------------------------
  +   call cpu_time(start_time)
   
      !---------------------------------------------
      ! Read number of points
  
positive.code: 
        call integrate(v, a, dt/2)  ! Compute v(t+dt/2)
        call integrate(r, v, dt)    ! Compute r(t+dt)
        call compute_force(m, r, a) ! Compute a(t+dt)
        call integrate(v, a, dt/2)  ! Compute v(t+dt)
  
        t = t + dt
  
negative.code:    deallocate(v)
     deallocate(a)
  
     !---------------------------------------------
     ! Elapsed time
     !---------------------------------------------
     call cpu_time(total_time)
     total_time = total_time - start_time
     print*, 'Elapsed time:', total_time
  
  end program n_body
  