现在要算一个内弹道过程，其中有六个微分方程和一个代数方程，它们都互相关联，若用龙格库塔（matlab）来计算，出现了这点疑问：
[latex]$$ \frac{dZ}{dt}=ap^n $$[/latex]   ，其中a是一个常数，p是t的函数，p的计算表达式里面不含微分，通过初值就可以直接求出当前的p（t）的值，我的问题在于p不是微分方程故不可能用龙格库塔来迭代计算，这里的Z需要用龙格库塔，那么就要求Z的k1，k2，k3，k4，在dZ/dt=f的表达式里面，p是一个t的函数，那k1=f（xi，yi），k2=f（xi，yi+k1/2）......这里面的x，y都是自变量，那dz/dt的p也是自变量吗？怎么代入呢....不知道怎么写，求助 返回小木虫查看更多