def METHOD_NAME(self, kwargs, expected):
    """
    Test calculated values based on Figure 1 of Bellan 2012
    (DOI: https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2012JA017856).
    """
    # theta and k values need to be single valued for this test to function
    # correctly
    cs = cs_(kwargs["T_e"], kwargs["T_i"], kwargs["ion"])
    va = va_(kwargs["B"], kwargs["n_i"], ion=kwargs["ion"])
    wci = wc_(kwargs["B"], kwargs["ion"])
    beta = (cs / va).value ** 2
    if not np.isclose(beta, 0.4, atol=1e-4):
        pytest.fail(
            f"The Bellan 2012 paper requires a 'beta' value of 0.4 and the test "
            f"parameters yielded {beta:.6f}."
        )
    Lambda = (kwargs["k"] * va / wci).value ** 2
    if not np.isclose(Lambda, 0.4, atol=1e-4):
        pytest.fail(
            f"The Bellan 2012 paper requires a 'Lambda' value of 0.4 and the test "
            f"parameters yielded {Lambda:.6f}."
        )
    ws = two_fluid(**kwargs)
    for mode, val in ws.items():
        norm = (np.absolute(val) / (kwargs["k"] * va)).value ** 2
        assert np.isclose(norm, expected[mode])