解二阶微分方程组

maplev12

这个微分方程组用Maple是怎么解的呀？

maplelab

很多常微分方程(组)可考虑求数值解.
特别是含有参数的常微分方程(组).

以下方法仅供参考,不一定合乎需要.

1.  定义一个通用的含有参数的常微分方程(组)求解程序odesolve.

2.  对不同的参数赋值,调用odeslove程序求相对应的解.

odesolve程序如下:

1. restart;
   
2. odesolve:=proc(para::set)
   
3. local test;
   
4. test:=proc(eqns::set,rag::`..`,var::list(list),opt1::set:={},opt2::set:={})
   
5. local sol;
   
6. sol:=dsolve(eqns,type=numeric,range=rag,op(opt1));
   
7. return plots:-odeplot(sol,var,op(opt2));
   
8. end proc:
   
9. return test(eval(args[2],para),args[3..-1]);
   
10. end proc:

复制代码

附件是maple例子文件:

maplev12

太感谢了…这样解很好～我用级数法求解的来着。

1. dsolve(eqs,{u(t),v(t),w(t)},type=series);

复制代码