2楼: Originally posted by nagami at 2014-11-13 10:56:57
你的初值，使指数项在开始时变得很大，这样方程其实就是刚性的。
换刚性求解器+打靶法，应该可以

Variable c=[0,1],a=[47,], b=[23,], a'=[,0], b'=[,0];
Plot a,b;
ODEFunction a''=0.1883*a*(3.38/(b+1.4)-3.42*exp(-6.4*b)+1);
            b''=0.1321*a*(3.38/(b+1.4)-3.42*exp(-6.4*b)+1);

4楼: Originally posted by dingd at 2014-11-13 14:34:14
1stOpt求解：

Variable c=,a=, b=, a'=, b'=;
Plot a,b;
ODEFunction a''=0.1883*a*(3.38/(b+1.4)-3.42*exp(-6.4*b)+1);
            b''=0.1321*a*(3.38/(b+1.4)-3.42*exp(-6.4*b)+1);


常微分方程(边值 ...

3楼: Originally posted by lockyli at 2014-11-13 14:18:15
果然是大神，一眼看出问题了，我又看了下我的微分方程，方程①和②中指数想应该为3.42*exp（-6.4b），你能帮我再看下，十分感谢了！...

6楼: Originally posted by nagami at 2014-11-13 17:18:48
计算的结果与3楼这个仁兄的一样，a，b其实画出来是几乎一样的，方程的是形式几乎相同

11.jpg
...

7楼: Originally posted by nagami at 2014-11-13 17:21:35
“这位仁兄”，打字太快没留意，表示尊敬我得更正下...

8楼: Originally posted by lockyli at 2014-11-13 19:47:35
兄弟，你能告诉下你用的软件名以及程序代码吗？因为我还需要解决一些类似的问题，也就是方程组的参数会改变一些，小弟先在此谢谢了。...