def METHOD_NAME(self):
    self.sim.integrator = "ias15"
    jupyr = 11.86*2.*math.pi
    self.sim.dt = jupyr*1e-7
    e0 = self.sim.energy()
    self.assertNotEqual(e0,0.)
    self.sim.integrate(self.sim.dt*1.001)
    self.sim.dt = jupyr
    self.sim.integrate(1e3*jupyr)
    e1 = self.sim.energy()
    self.assertLess(math.fabs((e0-e1)/e1),1e-14)