在R中有符号的ODE求解器吗？(ODE =常微分方程)

Question

问题：在R中有符号ODE求解器吗？(ODE = 常微分方程)

恐怕没有..。但让我从专家那里确认..。

例如，解决：

> (5x-6)^2 y' = 5(5x-6) y - 2

这里:y-未知函数，y‘-它的导数

(用手求解很容易:y= 1/(5 (5x-6) ) + C* (5x-6)，但我想从R中得到这个答案)。

我所知道的：

1)有数值(而非符号)求解者：

我知道有像库(DeSolve)这样的数字代码求解器，参见这里的答案：R语言能找到一阶微分方程的一般解吗？

2)有一些符号包: (但它们似乎不包含ODE求解器)

在R中有符号包，如packages和rSymPy，以及基R中的一些基本符号计算，参见：https://stats.stackexchange.com/questions/4775/symbolic-computation-in-r/4778

3)简要概述R：https://cran.r-project.org/web/views/DifferentialEquations.html的各种微分方程求解器

然而，我找不到求和乐谱的解(()

Answer 1

我在Ryacas**，上玩过，您可以在不需要太多工作的情况下为一些简单的ODE获得象征性的解决方案。不幸的是，YACAS无法为您的示例ODE找到解决方案。然而，根据您正在探索的ODEs，这可能仍然是有用的。如果没有，我很乐意删除这个职位。**

作为初始的简单示例，让我们考虑下面的ODE：y'' + y = 0

1. 加载库 图书馆(Ryacas)；
2. 由于Ryacas只是YACAS的一个接口，所以我们可以使用YACAS的OdeSolve来解决ODE OdeSolve(y‘\’+y == 0)#表达式(C70* exp(x * complex_cartesian(0，-1)) + C74 * exp(x *# complex_cartesian(0，1) 这给出了正确的解决方案const * exp(- ix) + const * exp(+ ix)。
3. 不幸的是，在使用特定示例时，OdeSolve无法找到解决方案： OdeSolve(y‘== (5 * (5 *x- 6) *y- 2)/(5 *x- 6)^2) #表达式(y( 2) - (5 *x- 6) *y(0)-2)/ (5 *x- 6)^2 ) 当我们使用YACAS在线演示时也会发生同样的情况。