[laosam280]

好了，今天可以解答这个问题了，不过为了简便，我没有对每个初值修改方程，回头你自己去修改吧。
结果和代码都在附件里面。顺便说一声，因为你的是非线性方程，步长最好自己控制，否则计算时间很长，我用了50个点，就算了61.64seconds