initial_guess = (2.5, 1.0, 0.3)
bounds = (0.0, [5.0, 5.0, 2.0])

params, covariance = curve_fit(
    decay_model,
    x,
    y,
    p0=initial_guess,
    bounds=bounds,
)