标签:dfun return 非线性 线割 print num plt fun 迭代法
from matplotlib import pyplot as plt
import numpy as np
def fun(x):
return x**3-2*x-5
def dfun(x):
return 3*x*x-2
def newton(fun,dfun,a,b,eps):
err = 1
x = b
k = 0
lada = 1
x_r = []
x_r.append(x)
x = x - fun(x) / dfun(x)
k = k + 1
while err > eps:
x_r.append(x)
x = x - (fun(x)/fun(x_r[-2]))*(x-x_r[-2])/(fun(x)/fun(x_r[-2])-1)
err = fun(x)
k = k + 1
print('牛顿两点线割法的迭代次数为{:d}次'.format(k))
print('x的迭代过程:')
num = 0
for i in x_r:
print('{:.7f}'.format(i), end=' ')
num = num+1
if num%5 == 0:
print(' ')
print(' ')
return x
x = newton(fun, dfun, 0, 10, 0.0000001)
print('牛顿两点线割迭代法的求解结果为{:.7}'.format(x))
x1 = np.linspace(-10, 10, 1000)
y1 = fun(x1)
plt.plot(x1, y1)
plt.show()
标签:dfun,return,非线性,线割,print,num,plt,fun,迭代法 来源: https://blog.csdn.net/weixin_40653652/article/details/112133072
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。