使用Maxima与abs_integrate计算多项式绝对值积分遇阻求助

abs_integrate包对三次及以上多项式的绝对值积分存在处理局限：这类多项式在[-1,1]区间内的实根求解逻辑超出了包的自动处理能力，而绝对值积分的核心是在根的位置拆分区间、去除绝对值符号，包无法完成这一步时就会直接返回未求值的积分表达式。低次多项式根的数量少且易求解，因此能正常计算。

• 手动拆分积分区间，通过先找根再分段积分的方式实现计算：
  1. 先求解多项式在[-1,1]内的所有实根，以此拆分积分区间；
  2. 在每个子区间内判断多项式的正负，去除绝对值后分别积分再求和。
     示例代码：

  load("abs_integrate");
  normaU(p) := block(
      /* 获取区间内的实根并排序 */
      real_roots = realroots(p, x, -1, 1),
      sorted_roots = sort(real_roots),
      intervals = cons(-1, append(sorted_roots, 1)),
      total = 0,
      /* 遍历每个子区间计算积分 */
      for i from 1 to length(intervals)-1 do (
          a = intervals[i],
          b = intervals[i+1],
          mid = (a + b)/2,
          /* 用区间中点判断多项式符号 */
          sign_val = sign(subst(x=mid, p)),
          total = total + integrate(sign_val * p, x, a, b)
      ),
      total
  );
  normaU(x^3 - 1/4);
  

• 用realroots替代solve获取实根：realroots(p, x, -1, 1)会直接返回[-1,1]内的实根列表，避免solve返回复数根的冗余处理，效率更高。
• 验证结果：可以用数值积分函数quad_qags(abs(x^3 - 1/4), x, -1, 1)得到数值结果，与符号积分结果对比确保正确性。

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