new_contents: ##########################TESTS ON MULTIPLE ORBITS#############################
  import numpy
  from galpy import potential
  
  # Test that integrating Orbits agrees with integrating multiple Orbit instances
  def test_integration_1d():
      from galpy.orbit import Orbit, Orbits
      times= numpy.linspace(0.,10.,1001)
      orbits_list= [Orbit([1.,0.1]),Orbit([0.1,1.]),Orbit([-0.2,0.3])]
      orbits= Orbits(orbits_list)
      # Integrate as Orbits
      orbits.integrate(times,
                       potential.toVerticalPotential(potential.MWPotential2014,1.))
      # Integrate as multiple Orbits
      for o in orbits_list:
          o.integrate(times,
                      potential.toVerticalPotential(potential.MWPotential2014,1.))
      # Compare
      for ii in range(len(orbits)):
          assert numpy.amax(numpy.fabs(orbits_list[ii].x(times)-orbits.x(times)[ii])) < 1e-10, 'Integration of multiple orbits as Orbits does not agree with integrating multiple orbits'
          assert numpy.amax(numpy.fabs(orbits_list[ii].vx(times)-orbits.vx(times)[ii])) < 1e-10, 'Integration of multiple orbits as Orbits does not agree with integrating multiple orbits'
      return None
      
  
subject: Test that 1D orbit integration agrees between Orbit and Orbits