求一个微分方程组的解
 

遇到一个微分方程组的求解问题：
我编的程序如下：
syms A B C D E F G y z t
eq1='Dy=A+B*y-C*y^2';
eq2='Dz=z*(-B+2*C*y-2*D)+z*(((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^0.5+((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^(-0.5))-exp(-0.05*t)';
[y,z]=dsolve('eq1','eq2','y(0)＝2000','z(0)=0','t')

可是运行以后，程序报错
??? Error using ==> maple
at offset 2, `}` unexpected

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

不知道啥意思！请教论坛的高手！谢谢！

回复: 求一个微分方程组的解
 

帮你顶一下，我也很关心这个问题!

回复: 求一个微分方程组的解
 

[QUOTE=zoom;38595]遇到一个微分方程组的求解问题：
我编的程序如下：
syms A B C D E F G y z t
eq1='Dy=A+B*y-C*y^2';
eq2='Dz=z*(-B+2*C*y-2*D)+z*(((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^0.5+((E*z*exp(0.05*t)+(E^2*z^2*exp(0.1*t)-G)^0.5)/F)^(-0.5))-exp(-0.05*t)';
[y,z]=dsolve('eq1','eq2','y(0)＝2000','z(0)=0','t')

可是运行以后，程序报错
??? Error using ==> maple
at offset 2, `}` unexpected

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

不知道啥意思！请教论坛的高手！谢谢！[/QUOTE]
可能是eq2中少了一个括号。
太长了，不容易发现。

回复: 求一个微分方程组的解
 

[QUOTE=reader1;38632]可能是eq2中少了一个括号。
太长了，不容易发现。[/QUOTE]

我检查过了，应该没有少括号！但还是谢谢你！