解二阶微分方程组
这个微分方程组用Maple是怎么解的呀? 很多常微分方程(组)可考虑求数值解.特别是含有参数的常微分方程(组).
以下方法仅供参考,不一定合乎需要.
1.定义一个通用的含有参数的常微分方程(组)求解程序odesolve.
2.对不同的参数赋值,调用odeslove程序求相对应的解.
odesolve程序如下:
restart;
odesolve:=proc(para::set)
local test;
test:=proc(eqns::set,rag::`..`,var::list(list),opt1::set:={},opt2::set:={})
local sol;
sol:=dsolve(eqns,type=numeric,range=rag,op(opt1));
return plots:-odeplot(sol,var,op(opt2));
end proc:
return test(eval(args,para),args);
end proc:
附件是maple例子文件:
太感谢了…这样解很好~我用级数法求解的来着。
dsolve(eqs,{u(t),v(t),w(t)},type=series);
页:
[1]