1. 原理
2. 代码
def newton_solve(N):
x0 = N / 2
x1 = x0 - (x0 ** 2 - N) / (2 * x0)
max_iter = 1000
iter = 0
while iter < max_iter and abs(x0 - x1) > 1e-5:
x0 = x1
x1 = x0 - (x0 ** 2 - N) / (2 * x0)
iter += 1
return x1
res = newton_solve(2)
print(res)