SageMath表达式无法求值：如何获取预期数值结果？

你的代码核心问题是符号变量的绑定逻辑错误：先定义了符号变量u_5和beta，之后给这两个变量赋值数值时，只是将变量名重新指向了数值，但最初构建的表达式y仍然引用的是原始符号变量。执行y.subs(u==u_5)时，仅替换了u为当前u_5的数值，而表达式内的符号变量u_5和beta并未被替换，导致结果保留符号项。

以下是两种可行的修正方案：

方案一：一次性替换所有符号变量 #

直接在subs中传入所有变量的取值字典，确保所有符号都被替换：

var('u_5 u beta')
y = u - u_5 + beta*(exp(u)-u-4)
# 传入所有变量的数值映射
p = y.subs({u: -4.0, u_5: -4.0, beta: 1.0})
# 转换为数值结果
p.n()

方案二：调用表达式时直接传入参数 #

SageMath的表达式对象支持直接调用传入参数，自动替换所有符号变量：

var('u_5 u beta')
y = u - u_5 + beta*(exp(u)-u-4)
# 直接传入所有变量的数值计算
p = y(u=-4.0, u_5=-4.0, beta=1.0)
p.n()

两种方案执行后都会得到数值结果：0.0183156388887342（本质是exp(-4.0)的近似值，因为代入后-u_5 -4.0部分抵消为0）。

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