Python近似圆周率时无输出与math domain error问题求助

1. return位置导致无输出的原因 #

把return G放在print语句上方时，函数执行到return就会直接终止并返回值，后续的print代码根本不会被执行，所以控制台看不到任何输出。必须把print放在return前面，才能让打印逻辑正常运行。

2. math domain error的根源 #

你的循环逻辑存在边界漏洞：

• 循环条件是while x2 <= 1，但每次循环的第一步就会将x2更新为x1 + d
• 当循环到最后一轮时，x1会被赋值为1，此时进入循环后，x2 = 1 + d，这个值必然大于1
• 把大于1的x传入halfCircle函数后，1 - pow(x,2)会变成负数，而math.sqrt()无法处理负数输入，直接抛出math domain error

修复后的代码 #

同时修复了字符串拼接的错误（浮点数不能直接与字符串用+拼接），优化了循环逻辑：

import math

def halfCircle(x):
    # 提前判断x的范围，避免sqrt处理负数报错
    if abs(x) > 1:
        return 0.0
    # 因为x²<=1，所以1-x²非负，不需要abs
    y = math.sqrt(1 - pow(x, 2))
    return y

def areaUnderGraph(N):
    x1 = -1.0
    G = 0.0
    d = 2 / N
    
    # 改用固定次数的for循环，避免浮点数边界判断的坑
    for _ in range(N):
        x2 = x1 + d
        y1 = halfCircle(x1)
        y2 = halfCircle(x2)
        
        # 取区间两端较小的y值计算矩形面积（下和近似）
        a = d * min(y1, y2)
        G += a
        x1 = x2
    
    # 用str()转换浮点数，或者用f-string更简洁
    print(f"The Area under the graph is {G}")
    # 当前计算的是半圆面积（π/2），乘以2得到π的近似值
    print(f"Approximation of π: {G * 2}")
    return G

areaUnderGraph(100)

额外说明 #

• 原代码的while循环容易因浮点数精度问题出现边界判断错误，改用固定N次的for循环更可靠，毕竟我们本来就是要把区间分成N块
• 你当前计算的是下和（用每个区间两端较小的y值计算矩形面积），后续可以实现上和（取区间两端较大的y值），取两者的平均值后乘以2，近似结果会更接近真实的π值

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