def METHOD_NAME(lhaaso1):
    m = lhaaso1[0].spectral_model("KM2A")
    dnde, dnde_err = m.evaluate_error(50 * u.TeV)
    assert dnde.unit == "cm-2 s-1 TeV-1"
    assert_allclose(dnde.value, 0.33e-16, rtol=1e-3)
    assert_allclose(dnde_err.value, 0.05e-16, rtol=1e-3)
    m = lhaaso1[11].spectral_model("WCDA")
    dnde, dnde_err = m.evaluate_error(3 * u.TeV)
    assert dnde.unit == "cm-2 s-1 TeV-1"
    assert_allclose(dnde.value, 0.69e-13, rtol=1e-3)
    assert_allclose(dnde_err.value, 0.16e-13, rtol=1e-3)