如何为每次迭代更新 ODE 求解器中的初始条件？

Question

我希望将初始值（pCH3OH_in、pH2O_in、pCO_in、pH2_in、pCO2_in、t）替换为 pCH3OH、pH2O、pCO、pH2、pCO2 的最后值，即 pCH3OH[-1]、pH2O[-1]等， ODES系统求解后得到。 我想提一下代码运行正常，但我想进行此调整。

Xi_feed = [0, 0, 0.110005, 0.782425, 0.10757]
press = np.linspace(5e06, 15e06, 6)
temp = np.linspace(473.15, 623.15, 7)

for i, p in enumerate(press):
    for j, t in enumerate(temp):  
        pCH3OH_in = Xi_feed[0]*p
        pH2O_in = Xi_feed[1]*p
        pCO_in = Xi_feed[2]*p
        pH2_in = Xi_feed[3]*p
        pCO2_in = Xi_feed[4]*p

        L = np.linspace(0, length[0])

        sol = solve_ivp(func, [L[0], L[-1]], np.array([pCH3OH_in, pH2O_in, pCO_in, pH2_in, 
        pCO2_in, t], dtype=object), t_eval = L) # call the function
        pCH3OH, pH2O, pCO, pH2, pCO2, T = sol.y