[求助]解常微分方程！！

[hitzhjtopku]

我的程序为：
a=0.03755;
g=9.8;
m=0.4387;
j1=0.25*m*a*a;
j2=j1;
j3=0.5*m*a*a;
x=pi/6;
eq1='(j1+m*a*a)*D2x+(j3+m*a*a)*z*sin(x)*Dy+m*g*a*cos(x)=j1*Dy*Dy*sin(x)*cos(x)';
eq2='j1*D(Dy*sin(x))+j1*Dx*Dy*cos(x)-j3*z*Dx=0';
eq3='(j3+m*a*a)*Dz=m*a*a*Dy*Dx*sin(x)';
[x,y,z]=dsolve('eq1','eq2','eq3')


运行完后提示为：??? Error using ==> maple
at offset 2, `}` unexpected

Error in ==> dsolve at 145
var_set = maple([var_set ' intersect ' var_set]);

Error in ==> chungun1 at 11
[x,y,z]=dsolve('eq1','eq2','eq3')


为什么啊？我调一下午了。。。。。。。。:ft:

[hitzhjtopku]

:lovely: 贵人帮忙啊

[walehenry]

dsolve是线性常微分方程的求解器，你这不是这种类型，所以出错
你可以把微分方程组编为一个函数，用0de45来求解

[hitzhjtopku]

你可以把微分方程组编为一个函数，用0de45来求解???
用龙阁酷达法？
大侠！帮忙示范下！

[hitzhjtopku]

是不是要用差分的方法解决非线性微分方程

[hitzhjtopku]

高手帮忙啊

[hitzhjtopku]

求助！！！！！！