python牛顿迭代公式_牛顿迭代法Python实现

例1：给定方程$$f(x) = x^2 + sin x - 1 = 0$$，判别该方程有几个实根，并用牛顿法求出方程所有实根，精确到$$10^{-4}$$.

解：利用画图法观察根的所在区间为(-2,-1)和(0,1)，其中蓝色为$$y=f(x)$$的曲线，橘黄色的直线是$$y=0$$.



画图代码：import numpy as np

import matplotlib.pyplot as plt

def fun(x):

return x**2 + np.sin(x) - 1

def plotCurves():

x = np.linspace(-3, 3, 100)

y = fun(x)

y1 = 0.*x

plt.plot(x,y)

plt.plot(x,y1)

def main():

plotCurves()

if __name__ == '__main__':

main()

另外一种观察方法是，将原来的函数分解为$$y=sin x$$和$$y=1-x^2$$两个函数，然后观察它们的交点，如下图所示：import numpy as np