如何使用Python (SciPy)解出一定数量的时间步长的ODE？

import numpy as np
from scipy.integrate import solve_ivp

def ode_func(t, y):
    # 定义ODE的方程
    dydt = ...
    return dydt

t0 = 0  # 初始时间
t_end = 10  # 结束时间
num_steps = 100  # 时间步长数量
t_span = np.linspace(t0, t_end, num_steps)  # 时间跨度
y0 = ...  # 初始条件

sol = solve_ivp(ode_func, [t0, t_end], y0, t_eval=t_span)

t = sol.t  # 时间数组
y = sol.y  # 解数组