7楼: Originally posted by 文刀星 at 2016-01-13 09:54:11
提供一个思路参考：将19式两边对t求导，可以解出dy(19)/dt。于是将dY/dt作为微分方程组，其中Y=。剩下的y'作为参数，它们都是可以由Y表示的，假设令剩下的y'=Y1，则Y1=f(Y)。
在ode的目标方程中可以这么写：
func ...

2楼: Originally posted by zho123855 at 2016-01-10 16:14:46
1,尝试能否消元获得普通微分方程组，然后使用ode求解。
2,如果不能或者比较麻烦的话，直接尝试差分求解。