我想解决R中的两个嵌套微分方程。假设x'= x + a给出x（T），y'= y + x + b给出y（0），其中a和b是常数。第一个方程式与第二个方程式无关，因为它具有终止条件，所以可以通过反转时间用deSolve求解。由于第二方程式具有初始条件，因此需要使用第一方程式的解并随着时间的推移而求解。有什么方法可以存储用于x的解的值吗？
我尝试使用X <-out [，“ X”]，其中out是x方程的函数求解的输出，并将X插入y的函数求解中，但是我得到以下消息： / p>
func（）返回的导数的数量必须等于初始条件向量的长度。
我了解问题的本质，但我不知道要解决这个问题。