牛顿法解方程python

牛顿法是一种用于解方程的迭代方法，在Python中可以使用以下代码来实现牛顿法求解方程：

import numpy as np

def f(x):
    y = x**3 - x - 1  # 方程的表达式
    return y

def g(x):
    y = 3*x**2 - 1  # 方程的导函数
    return y

def newton_method():
    x0 = 1.5  # 初值
    e = 10**(-9)  # 误差要求
    L = 0  # 初始化迭代次数

    while abs(f(x0)) > e:  # 采用残差来判断
        x1 = x0 - f(x0) / g(x0)  # 迭代公式, x(n+1) = x(n) - f(x(n)) / f'(x(n))
        x0 = x1
        L += 1  # 统计迭代次数

    return x1

solution = newton_method()
print("方程的解为:", solution)

你可以根据需要修改方程的表达式和初值，运行以上代码即可得到方程的解。