使用Maxima求解非线性微分方程组：是否有可用工具包？

嘿，好问题！其实Maxima本身就自带处理非线性微分方程组的能力，不需要额外安装工具包——它的内置微分方程求解函数就能搞定这类问题。咱们结合你给出的示例来一步步看怎么操作：

首先，先把你的方程组在Maxima中定义好：

eqn1: A(r)*E(r)/r - 'diff(A(r), r) = 5 + 1/r;
eqn2: r*'diff(E(r)/r) - A(r)^2/r = 0;

第一步：简化第二个方程 #

先处理eqn2，我们可以先展开并化简其中的导数项。对E(r)/r求导后代入方程，再用ratsimp化简：

%eqn2_expanded: r*((diff(E(r),r)*r - E(r))/r^2) - A(r)^2/r = 0;
eqn2_simplified: ratsimp(%eqn2_expanded);

化简后会得到：

'diff(E(r), r) - E(r)/r - A(r)^2/r = 0

整理一下就是关于E(r)的一阶线性微分方程：

'diff(E(r), r) = (E(r) + A(r)^2)/r

第二步：求解E(r)的通解 #

用Maxima的ode2函数求解这个关于E(r)的方程：

sol_E: ode2(diff(E(r), r) = (E(r) + A(r)^2)/r, E(r), r);

你会得到类似这样的通解（%c是积分常数）：

E(r) = r*(%c + integrate(A(r)^2/r^2, r))

第三步：代入第一个方程求解A(r) #

把E(r)的通解代入eqn1，得到只含A(r)的微分方程：

eqn1_substituted: subst(E(r) = rhs(sol_E), eqn1);
eqn1_simplified: ratsimp(eqn1_substituted);

到这一步得到的方程可能会比较复杂，如果符号求解有难度，可以尝试用Maxima的desolve函数进一步尝试符号解；如果符号解不存在或者难以求出，还可以用数值方法求解——比如用龙格-库塔方法的rk函数，只需要给初始条件和求解区间即可。

补充：数值求解的思路 #

如果符号解卡壳了，数值求解是个不错的备选。比如假设初始条件A(1)=1，E(1)=1，求解区间r=1到r=10，可以用：

rk([diff(A,r) = A*E/r - 5 - 1/r, diff(E,r) = (E + A^2)/r], [A,E], [1,1], r, 1, 10, 0.1);

这个命令会输出一系列离散点的数值解，你可以用plot2d可视化结果。

总结一下：Maxima不需要额外工具包，靠内置的ode2、desolve（符号求解）和rk（数值求解）就能处理非线性微分方程组，关键是先简化方程组，逐步消元求解。

内容的提问来源于stack exchange，提问作者Astrolabe