test(tests/test_simulate.py): add more cases for some Michaelis-Mente…

…n model tests (#618)

This adds more "extreme" edge cases for both the simple and the "consuming" Michaelis-Menten model tests

tests/test_simulate.py

@@ -116,9 +116,9 @@ def test_get_y_conservation_in_equilibria() -> None:
         (cat0, sub0, keq, kcat, method)
         for cat0 in [0.3, 0.4]
         for sub0 in [0.01, 0.02]
-        for keq in [10.0, 100.0]
-        for kcat in [1e10, 1e11]
-        for method in ("RK23", "LSODA")
+        for keq in [1.0, 10.0, 100.0]
+        for kcat in [1e-1, 1e1, 1e10, 1e11, 1e13]
+        for method in ("RK23", "LSODA", "Radau", "BDF")
     ],
 )
 def test_simple_michaelis_menten(
@@ -164,9 +164,9 @@ def test_simple_michaelis_menten(
         (cat0, sub0, keq, kcat, method)
         for cat0 in [0.3, 0.4]
         for sub0 in [0.01, 0.02]
-        for keq in [10.0, 100.0]
-        for kcat in [1e10, 1e11]
-        for method in ("RK23", "LSODA")
+        for keq in [1.0, 10.0, 100.0]
+        for kcat in [1e-1, 1e1, 1e10, 1e11, 1e13]
+        for method in ("RK23", "LSODA", "Radau", "BDF")
     ],
 )
 def test_consuming_michaelis_menten(