Scratch中计算∫ₐᵇxⁿdx定积分的方法是否可靠？

背景 #

出于兴趣在Scratch中进行编程练习，尝试实现形如$\int_a^b x^n dx$（$n$为常数）的幂函数定积分计算。已编写自定义幂运算积木和积分计算积木，测试案例$\int_3^4a2da$（理论结果约12.333）时结果符合预期，但不确定该实现的可靠性，希望得到优化方案。

现有实现代码 #

自定义幂运算积木 #

define pow(base)(power)
set [result v] to [1]
set [base v] to (base)
set [power v] to (round (power))
set [done v] to [0]
if <(power) < [0]> then
set [base v] to [((1) / (base))]
set [power v] to [-1(() *(round (power))]
end
repeat until <[done v] = [1]>
if <(([power v]) mod (2)) = [0]> then
set [result v] to (([result v]) * ([base v]))
change [power v] by (-1)
end
if <(([power v]) mod (2)) > [0]> then
set [done v] to [1]
end
if <<not <<[done v] = [1]>>>
set [base v] to (([base v]) * ([base v]))
set [power v] to (([power v]) / (2))
end
end

自定义积分计算积木 #

define area from x to y of the form a^zda(x)(y)(z)
set [x v] to (x)
set [y v] to (y)
set [z v] to (z)
pow (x) ((z)+(1))
set [re2 v] to (result)
pow (y) ((z)+(1))
set [re3 v] to (result)
set [re4 v] to (((re3)-(re2))/(3))

现有实现的可靠性问题 #

1. 积分公式硬编码错误：积分积木中将分母固定为3，这仅适配$n=2$的情况，完全不符合通用幂函数积分公式$\frac{1}{n+1}(b^{n+1}-a{n+1})$，一旦$n≠2$，计算结果必然错误。
2. 幂函数实现缺陷：
   • 强制对幂次取整，无法处理非整数幂（如$n=0.5$、$n=-0.5$等场景），直接限制了适用范围。
   • 快速幂逻辑错误：当前循环逻辑无法正确计算大部分整数幂（例如$3^3$、$42$都会得到错误结果），仅测试案例巧合得到正确值。
   • 负幂处理存在语法错误：set [power v] to [-1(() *(round (power))]括号不匹配，运行时会报错，无法正确计算负幂场景。
3. 未覆盖特殊情况：当$n=-1$时，积分公式变为$\ln|b| - \ln|a|$，现有代码未处理该情况，会出现除以0的错误。

优化方案 #

1. 修复积分公式硬编码问题 #

将积分公式中的分母替换为$n+1$，并添加$n=-1$的特殊情况处理：

define area from x to y of the form a^zda(x)(y)(z)
set [x v] to (x)
set [y v] to (y)
set [z v] to (z)
# 处理n=-1的特殊积分情况
if <(z) = [-1]> then
  # 检查输入合法性：对数的真数必须为正
  if <(x) > [0]> and <(y) > [0]> then
    set [re4 v] to (([ln v] of (y)) - ([ln v] of (x)))
  else
    set [re4 v] to [输入错误：n=-1时，上下限必须为正数]
  end
else
  # 通用幂函数积分公式
  pow (x) ((z)+(1))
  set [re2 v] to (result)
  pow (y) ((z)+(1))
  set [re3 v] to (result)
  set [re4 v] to (((re3)-(re2))/((z)+(1)))
end

2. 重构幂运算积木 #

方案一：直接使用Scratch内置幂函数

Scratch自带[power v] of (base)积木，支持整数、非整数及负幂运算，比自行实现更可靠，可直接替换自定义pow积木：

# 替代自定义pow积木，直接用内置功能
set [re2 v] to ([power v] of (x) raised to ((z)+(1)))
set [re3 v] to ([power v] of (y) raised to ((z)+(1)))

方案二：修正自定义快速幂（仅支持整数幂）

如果是为了练习快速幂实现，修正逻辑为标准快速幂流程，并修复负幂处理的语法错误：

define pow(base)(power)
set [result v] to [1]
set [current_base v] to (base)
set [current_power v] to (power)
# 处理负幂
if <(current_power) < [0]> then
  set [current_base v] to ((1)/(current_base))
  set [current_power v] to ((-1)*(current_power))
end
# 标准快速幂循环
repeat until <(current_power) = [0]>
  if <((current_power) mod (2)) = [1]> then
    set [result v] to ((result)*(current_base))
  end
  set [current_base v] to ((current_base)*(current_base))
  set [current_power v] to (floor((current_power)/(2)))
end

3. 增加输入合法性检查 #

添加对非法输入的判断，避免运行错误：

• 当$n$为负数且上下限包含0时，提示“输入错误：负幂运算中上下限不能为0”。
• 当$n=-1$且上下限存在非正数时，提示“输入错误：对数运算的真数必须为正数”。

4. 多案例验证 #

通过以下案例验证优化后的代码正确性：

• $\int_1^2x1dx$：理论结果1.5
• $\int_0^1x{0.5}dx$：理论结果约0.6667
• $\int_2^4x{-1}dx$：理论结果约0.6931

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