请求提供R中使用rim包调用Maxima的行内代码示例

我们在建模课程中使用Maxima，希望让学生在R中通过rim包以行内代码（非RMarkdown）调用Maxima完成符号运算，同时在R中完成绘图等数值计算工作。目前测试时遇到调用问题：比如计算矩阵行列式时，直接用R中存储的变量名无法运行，必须使用Maxima输出行号（如%o12）才能执行。下面是解决该问题的方法及更多完整可行的交互示例：

核心问题说明 #

maxima.get()函数是获取Maxima的运算结果并返回给R，但不会在Maxima的环境中定义对应变量。如果要在后续Maxima命令中使用变量名，需要先在Maxima环境中显式定义变量，而不是仅将结果赋值给R变量。

完整示例代码 #

1. 基础符号运算：解方程与积分 #

library(rim)

# 解符号方程：逻辑增长模型的平衡点
maxima.get("solve((1 + r*(1 - x/K))*x = x, x);")

# 计算不定积分，并在R数据框中生成数值结果
myint <- maxima.get("integrate(1/(1-x), x);")
# 构造R数据框（避开x=1的奇点）
df <- data.frame(x = seq(0, 0.9, by = 0.1))
# 将符号积分结果在R环境中求值
maxima.eval(myint, code = TRUE, envir = df)

2. 矩阵符号运算（解决行列式调用问题） #

要在Maxima中使用变量名，需先在Maxima环境中定义变量：

# 第一步：在Maxima环境中定义矩阵变量mymat
maxima.get("mymat: matrix([0, 1, a], [1, 0, 1], [1, 1, 0]);")

# 第二步：直接使用变量名计算行列式，此时可正常运行
maxima.get("determinant(mymat);")

# 扩展：计算矩阵的逆
maxima.get("invert(mymat);")

# 扩展：符号矩阵乘法
maxima.get("mymat . transpose(mymat);")

3. 符号函数结合R数值绘图 #

将Maxima得到的符号函数转为R可执行的函数，再进行绘图：

# 定义符号函数：逻辑增长模型的导数
maxima.get("dlogistic: diff((1 + r*(1 - x/K))*x, x);")

# 将符号导数转为R函数
dlogistic_fun <- maxima.eval(dlogistic, code = TRUE, function = TRUE)

# 设定参数值，绘制函数曲线
r <- 0.5
K <- 10
x_seq <- seq(0, K, by = 0.1)
plot(x_seq, dlogistic_fun(x_seq, r = r, K = K), type = "l", 
     xlab = "x", ylab = "dx/dt", main = "逻辑增长模型的导数曲线")

4. 多步符号运算与数值验证 #

# 第一步：求解微分方程的符号解
maxima.get("ode2(diff(x,t) = r*x*(1 - x/K), x, t);")

# 第二步：将符号解转为R函数
logistic_sol <- maxima.eval("%", code = TRUE, function = TRUE)

# 第三步：设定初始条件与参数，生成数值结果
t_seq <- seq(0, 20, by = 0.5)
x0 <- 1
r <- 0.3
K <- 10
x_vals <- logistic_sol(t_seq, x0 = x0, r = r, K = K)

# 第四步：在R中绘图
plot(t_seq, x_vals, type = "l", xlab = "时间t", ylab = "种群数量x", 
     main = "逻辑增长模型的种群动态")

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