def test_Bss_2fast_approxiamtion():
t, yout, A_ssB, A_ssB_2fast, ydot = integrate_rd(
1.0, 1e-10, 220e5, 217)
assert np.allclose(yout[:, 0, 0], A_ssB_2fast)
def test_Bss_approxiamtion():
t, yout, A_ssB, A_ssB_2fast, ydot = integrate_rd(
1.0, 1e-3, 220, 217)
assert np.allclose(yout[:, 0, 0], A_ssB)
_assert_mp(t, ydot, yout)