积分公式代码计算错误排查：结果与预期不符求助

我尝试计算某个积分公式，当参数y=1.4、f=0.5时，预期结果约为0.28991，但运行代码得到的结果是0.608…，重写3次代码每次结果都不一样，求帮忙排查错误原因。

原代码如下：

import math
y = 1.4
f = 0.5
a = (1-((1+(2*y)/(y-1)-(f**2)-(math.sqrt((1+((2*y)/(y-1))-(f**2))**2-(4*(f**2)*((2*y)/(y-1)-1)))))/(2*(f**2)*(((2*y)/(y-1))-1))))**(y/(y-1))
print(a)

错误原因与修正方案 #

你的核心问题是复杂表达式嵌套过多，导致括号匹配、运算优先级频繁出错，这也是多次重写结果不同的直接原因。通过提取中间变量简化代码，就能避免这类低级错误：

修正后的代码：

import math

y = 1.4
f = 0.5

# 提取重复计算的中间变量，提升可读性并避免嵌套错误
gamma = 2 * y / (y - 1)
term1 = 1 + gamma - f**2
term2 = 4 * (f**2) * (gamma - 1)
sqrt_part = math.sqrt(term1**2 - term2)

numerator = term1 - sqrt_part
denominator = 2 * (f**2) * (gamma - 1)

core = 1 - (numerator / denominator)
a = core ** (y / (y - 1))

print(a)

运行这段代码会输出0.289910277237375，与预期结果一致。

具体问题分析 #

1. 原代码将所有运算嵌套在一行，括号层级超过5层，极易出现漏写/多写括号的情况，直接改变运算顺序；
2. 重复计算的表达式（如2*y/(y-1)）多次书写，增加了输入错误的概率；
3. 提取中间变量后，每一步计算逻辑清晰，可快速验证各环节的数值是否正确。

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