def METHOD_NAME():
    mag = m._saturn_magnitude(9.014989659493, 8.03160470546889, 0.1055, -26.224864126755417, -26.332275658328648, True)
    assert abs(-0.552 - mag) < 0.0005
    mag = m._saturn_magnitude(9.438629261423, 8.47601935508925, 1.8569, -8.016972372600742, -7.458535528992841, False)
    assert abs(0.567 - mag) < 0.0005
    mag = m._saturn_magnitude(9.026035315474, 10.1321497654765, 169.8958, 29.097094546361053, -26.672903610681395, False)
    assert abs(5.206 - mag) < 0.0005
    mag = m._saturn_magnitude(9.026035315474, 10.1321497654765, 169.8958, 29.097094546361053, -26.672903610681395, True)
    assert np.isnan(mag)
    args = [
        A[9.014989659493, 9.438629261423, 9.026035315474, 9.026035315474],
        A[8.03160470546889, 8.47601935508925, 10.1321497654765, 10.1321497654765],
        A[0.1055, 1.8569, 169.8958, 169.8958],
        A[-26.224864126755417, -8.016972372600742, 29.097094546361053, 29.097094546361053],
        A[-26.332275658328648, -7.458535528992841, -26.672903610681395, -26.672903610681395],
        A[True, False, False, True],
    ]
    magnitudes = m._saturn_magnitude(*args)
    expected = [-0.552, 0.567, 5.206, nan]
    np.allclose(magnitudes, expected, 0, 0.0005, equal_nan=True)