你能给出一个自适应步长scipy.integrate.LSODA函数的简单例子吗？

import numpy as np
from scipy.integrate import LSODA

def func(t, y, p):
    # 定义常微分方程组的函数
    dydt = np.zeros_like(y)
    dydt[0] = -p * y[0]  # 示例方程为 dy/dt = -p * y
    return dydt

# 设置初始条件和参数
y0 = np.array([1.0])
p = 0.5

# 创建 LSODA 求解器对象
solver = LSODA(func, t0=0.0, y0=y0, args=(p,))

# 设置求解器的参数
solver.set_initial_step(0.1)  # 初始步长
solver.set_f_params(p)  # 传递额外的参数给函数 func

# 求解常微分方程组
t_values = np.linspace(0, 10, 100)  # 时间范围
y_values = np.zeros((len(t_values), len(y0)))  # 存储解的数组
y_values[0] = y0

for i in range(1, len(t_values)):
    solver.integrate(t_values[i])
    y_values[i] = solver.y

# 打印结果
for t, y in zip(t_values, y_values):
    print(f"t = {t:.2f}, y = {y[0]:.4f}")