def METHOD_NAME():
    fitter = LevMarLSQFitter()
    b = BlackBody(3000 * u.K, scale=5e-17, output_units=u.Jy / u.sr)
    wav = np.array([0.5, 5, 10]) * u.micron
    fnu = np.array([1, 10, 5]) * u.Jy / u.sr
    b_fit = fitter(b, wav, fnu, maxiter=1000)
    assert_quantity_allclose(b_fit.temperature, 2840.7438355865065 * u.K)
    assert_quantity_allclose(b_fit.scale, 5.803783292762381e-17)