非线性非齐次常微分方程组的求解

wrx98

y"+2ay'+by=cx^2
x"+2dx'+ex=fxy+gcos(t)
其中x和y是时间t的函数，a,b,c,d,e,f,g都是常数，我用四阶龙格库塔法进行求解，首先将方程化为一阶的四个方程组，进行求解，但是由于系数a,b,c,d,e,f,数值很大，所以迭代几步就出现Nan,a,b,c,d,e,f在数值10^6次以上，
x'(0)=0
y'(0)=0
x(0)=0
y(0)=0
各位大哥大姐，如果有什么好的处理方法，可以给我介绍一下的，非常感激了哦，真的
我的联系电话是13858243461
毕业前最痛苦的日子，希望得到你的帮助

[ 本帖最后由 eight 于 2007-6-11 11:24 编辑 ]

sssssxxxxx921

那估计是病态方程组，龙格库塔法是不适合于求解病态方程的你用ode15s求解算算看看