有没有一种在Python中递归调用ODE变量的方法？

import numpy as np
from scipy.integrate import solve_ivp

# 定义ODE
def ode_func(t, y):
    return -y

# 定义递归函数
def recursive_ode_solver(t_span, y0, t_eval):
    def recursive_helper(t_span, y0, t_eval, current_t):
        if len(t_eval) == 0:
            return []
        t_eval_next = [t_eval[0]]
        sol = solve_ivp(ode_func, t_span, y0, t_eval=t_eval_next)
        result = sol.y.T.tolist()
        return result + recursive_helper(t_span, sol.y[:, -1], t_eval[1:], current_t + sol.t[-1])
    
    return recursive_helper(t_span, y0, t_eval, t_span[0])

# 初始条件
y0 = [1.0]

# 时间范围和评估点
t_span = (0, 5)
t_eval = np.linspace(t_span[0], t_span[1], 10)

# 递归求解ODE
result = recursive_ode_solver(t_span, y0, t_eval)
print(result)