def METHOD_NAME(dataset):
    fit = Fit()
    result = fit.run(datasets=dataset)
    model = dataset.models[0].spectral_model
    assert_allclose(model.amplitude.error, 1.9e-14, rtol=1e-2)
    model.amplitude.scan_n_values = 3
    model.amplitude.scan_n_sigma = 1
    model.amplitude.interp = "lin"
    profile = fit.stat_profile(
        datasets=dataset,
        parameter="amplitude",
    )
    ts_diff = profile["stat_scan"] - result.total_stat
    assert_allclose(
        model.amplitude.scan_values, [1.97e-13, 2.16e-13, 2.35e-13], rtol=1e-2
    )
    assert_allclose(ts_diff, [110.244116, 0.0, 110.292074], rtol=1e-2, atol=1e-7)
    value = model.parameters["amplitude"].value
    err = model.parameters["amplitude"].error
    model.amplitude.scan_values = np.array([value - err, value, value + err])
    profile = fit.stat_profile(
        datasets=dataset,
        parameter="amplitude",
    )
    ts_diff = profile["stat_scan"] - result.total_stat
    assert_allclose(
        model.amplitude.scan_values, [1.97e-13, 2.16e-13, 2.35e-13], rtol=1e-2
    )
    assert_allclose(ts_diff, [110.244116, 0.0, 110.292074], rtol=1e-2, atol=1e-7)