@@ -0,0 +1,327 @@
+import pytest%0Afrom numpy.testing import assert_allclose%0Aimport numpy as np%0Afrom scipy.integrate import RK23, RK45%0A%0A%0A@pytest.mark.parametrize(%22solver%22, %5BRK23, RK45%5D)%0Adef test_coefficient_properties(solver):%0A    assert_allclose(np.sum(solver.B), 1, rtol=1e-15)%0A    assert_allclose(np.sum(solver.A, axis=1), solver.C, rtol=1e-15)%0A

Add a test for Runge-Kutta coefficient properties