def METHOD_NAME(solvers):
    testsolver, scipysolver, ode = solvers
    state = testsolver.initialize(ode)
    time_scipy = scipysolver.t
    state_scipy = scipysolver.y
    np.testing.assert_allclose(state.t, time_scipy, atol=1e-13, rtol=1e-13)
    np.testing.assert_allclose(state.rv.mean[0], state_scipy[0], atol=1e-13, rtol=1e-13)