关于Python解方程

牛顿迭代源码：
```python
def d(fu, dx=1e-5):
    return lambda x: (fu(x + dx) - fu(x)) / dx  # 定义式求导

def solve(fu, prec=1e-6, maxit=10, l=-1e9, r=1e9, dx=1e-5):
    """

:param fu: 你要求解的函数
    :param prec: 精度
    :param maxit: 最大迭代次数
    :param l: 二分左边界
    :param r: 二分右边界
    :return: 方程的根
    """
    monotony = fu(r) > fu(l)  # 单调性
    while r - l > 0.5: 
        mid = (l + r) / 2
        if fu(mid) > 0:
            if monotony:
                r = mid
            else:
                l = mid
        else:
            if monotony:
                l = mid
            else:
                r = mid
    x = (l + r) / 2  # 猜测值
    times = 0
    while abs(fu(x)) > prec:
        if times > maxit:
            raise RuntimeError("Max iteration exceed!")  # 最大迭代次数超过限制
        times += 1
        x -= fu(x) / d(fu,dx)(x)  # 迭代主体
    return x

if __name__ == '__main__':
    print(solve(lambda x: x * x - 3, prec=1e-15))

```
这样求得$\sqrt{3} = 1.7320508075688772$(虽然求得是负根)
如果想求正根可以用`solve(lambda x: x * x - 3, prec=1e-15, l=0)`（锁定最小值为0，当然就是求正根了） _强调_

Programmer-1048576/mathprolib

内容风险标识

评论 (29)

Programmer-1048576/mathprolib .gitee-modal { width: 500px !important; }

内容风险标识