solution = solve_ivp(
    exponential_decay,
    t_span=time_span,
    y0=[INITIAL_VALUE],
    t_eval=sample_times,
    method="RK45",
    rtol=1e-8,
    atol=1e-10,
)