使用Regula Falsi（试位法）求函数根时遇ZeroDivisionError求助

核心问题1：测试函数与初始区间不满足试位法前提 #

试位法（Regula Falsi）的生效条件是：区间[a,b]内函数连续，且两端点函数值异号（f(a)*f(b)<0），以此保证区间内存在实根。但你测试的函数f(x)=x²+2x+3：

• 判别式Δ=2²-4×1×3=-8<0，无实根；
• 在区间(-2,3)上，f(-2)=3>0，f(3)=18>0，两端点函数值同号，完全不满足试位法的初始要求。

核心问题2：ZeroDivisionError的触发原因 #

按你给出的代码和测试用例，第一次迭代时分母f(b)-f(a)=18-3=15≠0，理论上不会触发除零错误。出现错误的大概率场景是：

• 测试函数定义错误：比如误将Python幂运算符**写成按位异或^（^不是幂运算），或函数表达式写错，导致某一步中f(a)与f(b)值相等；
• 初始区间传入错误：比如传入的区间端点相同（如(-2,-2)），直接导致f(b)-f(a)=0。

修复建议 #

1. 更换合法测试用例：选择有实根且满足端点异号的函数和区间，例如测试f(x)=x²-2，区间取(1,2)（f(1)=-1<0，f(2)=2>0）。
2. 修正函数定义语法：确保幂运算使用**，例如：

   f = lambda x: x**2 + 2*x + 3  # 语法正确，但无实根，仅作示例
   

3. 添加输入校验与代码优化：在函数中加入初始条件检查，拦截无效输入，同时减少重复计算：

   def regulafalsi(f, itv0, its=20):
       a = itv0[0]
       b = itv0[1]
       
       # 初始条件校验
       fa = f(a)
       fb = f(b)
       if fa * fb >= 0:
           raise ValueError("初始区间两端点函数值必须异号")
       if a == b:
           raise ValueError("区间两端点不能相同")
       
       for i in range(its):
           # 计算试位法的c值
           c = (fb * a - fa * b) / (fb - fa)
           fc = f(c)
           
           # 更新区间端点
           if fc * fa > 0:
               a = c
               fa = fc
           else:
               b = c
               fb = fc
       return c
   

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