import numpy as np
def func(x):
result = x - np.tan(x)
return result
def bisection(f,a,b,tol,maxit):
if a > b:
a, b = b, a
if np.abs(f(a))<=1.0e-16:
return a
if np.abs(f(b))<=1.0e-16:
return b
if f(a)*f(b) > 0:
print(“Zeros point may not exist”)
exit(1)
a1, b1 = a, b
for i in range(maxit):
bc = (a1+b1)/2
print(“it = {},[a,b] = [{},{}],f((a+b)/2) = {}”.format(i,a1,b1,f(bc)))
if np.abs(f(bc))<=1.0e-16:
return bc
elif (f(a1)*f(bc)<0):
b1 = bc
else:
a1 = bc
return bc
def main():
maxit = 100
a, b = 0.01, 60
tol = 1.0e-8
x=bisection(func,a,b,tol,maxit)
print(“x= {}”.format(x))
if name == ‘main’:
main()
3、结果
it = 0,[a,b] = [0.01,60],f((a+b)/2) = 36.20670948365492
it = 1,[a,b] = [0.01,30.005],f((a+b)/2) = 15.850580622920619
it = 2,[a,b] = [0.01,15.0075],f((a+b)/2) = 4.728146144446539
it = 3,[a,b] = [0.01,7.50875],f((a+b)/2) = 3.0488086231383953
it = 4,[a,b] = [0.01,3.759375],f((a+b)/2) = 4.965180853113451
it = 5,[a,b] = [0.01,1.8846874999999998],f((a+b)/2) = -0.44321739582592645
it = 6,[a,b] = [0.9473437499999999,1.8846874999999998],f((a+b)/2) = -4.993061970154832
it = 7,[a,b] = [1.416015625,1.8846874999999998],f((a+b)/2) = 14.19370489381378
it = 8,[a,b] = [1.416015625,1.6503515625],f((a+b)/2) = -25.04101886482726
it = 9,[a,b] = [1.53318359375,1.6503515625],f((a+b)/2) = 49.26910366880788
it = 10,[a,b] = [1.53318359375,1.591767578125],f((a+b)/2) = -118.61635689888206
it = 11,[a,b] = [1.5624755859375,1.591767578125],f((a+b)/2) = 159.67140020221845
it = 12,[a,b] = [1.5624755859375,1.5771215820312499],f((a+b)/2) = -1000.692164743844
it = 13,[a,b] = [1.5697985839843749,1.5771215820312499],f((a+b)/2) = 376.9823019098605
it = 14,[a,b] = [1.5697985839843749,1.5734600830078125],f((a+b)/2) = 1202.041886371078
it = 15,[a,b] = [1.5697985839843749,1.5716293334960936],f((a+b)/2) = -12139.058342277169
it = 16,[a,b] = [1.5707139587402343,1.5716293334960936],f((a+b)/2) = 2665.9689019298944
it = 17,[a,b] = [1.5707139587402343,1.571171646118164],f((a+b)/2) = 6828.644930977909
it = 18,[a,b] = [1.5707139587402343,1.5709428024291991],f((a+b)/2) = 31199.129846181575
it = 19,[a,b] = [1.5707139587402343,1.5708283805847167],f((a+b)/2) = -39748.58767499719
it = 20,[a,b] = [1.5707711696624755,1.5708283805847167],f((a+b)/2) = 289997.1272461388
it = 21,[a,b] = [1.5707711696624755,1.570799775123596],f((a+b)/2) = -92126.95115168976
it = 22,[a,b] = [1.5707854723930357,1.570799775123596],f((a+b)/2) = -270047.07122987835
it = 23,[a,b] = [1.5707926237583159,1.570799775123596],f((a+b)/2) = -7852130.792086863
it = 24,[a,b] = [1.5707961994409558,1.570799775123596],f((a+b)/2) = 602234.3925935033
it = 25,[a,b] = [1.5707961994409558,1.570797987282276],f((a+b)/2) = 1304519.461708586
it = 26,[a,b] = [1.5707961994409558,1.570797093361616],f((a+b)/2) = 3128850.159251721
it = 27,[a,b] = [1.5707961994409558,1.570796646401286],f((a+b)/2) = 10402990.014126034
it = 28,[a,b] = [1.5707961994409558,1.5707964229211209],f((a+b)/2) = -64045667.53043196
it = 29,[a,b] = [1.5707963111810384,1.5707964229211209],f((a+b)/2) = 24840906.13909882
it = 30,[a,b] = [1.5707963111810384,1.5707963670510796],f((a+b)/2) = 81161176.10896431
it = 31,[a,b] = [1.5707963111810384,1.5707963391160589],f((a+b)/2) = -607404991.2282902
it = 32,[a,b] = [1.5707963251485486,1.5707963391160589],f((a+b)/2) = 187356891.7622508
it = 33,[a,b] = [1.5707963251485486,1.5707963321323037],f((a+b)/2) = 541849931.0295086
it = 34,[a,b] = [1.5707963251485486,1.570796328640426],f((a+b)/2) = 10041096317.569733
it = 35,[a,b] = [1.5707963251485486,1.5707963268944873],f((a+b)/2) = -1293027499.757818
it = 36,[a,b] = [1.5707963260215179,1.5707963268944873],f((a+b)/2) = -2968291555.4557796
it = 37,[a,b] = [1.5707963264580025,1.5707963268944873],f((a+b)/2) = -8428021001.515761
it = 38,[a,b] = [1.5707963266762448,1.5707963268944873],f((a+b)/2) = -104924530373.4907
it = 39,[a,b] = [1.570796326785366,1.5707963268944873],f((a+b)/2) = 22207347960.62738
it = 40,[a,b] = [1.570796326785366,1.5707963268399268],f((a+b)/2) = 56338852158.332436
it = 41,[a,b] = [1.570796326785366,1.5707963268126464],f((a+b)/2) = 243336218459.12137
it = 42,[a,b] = [1.570796326785366,1.5707963267990062],f((a+b)/2) = -368927491579.71277
it = 43,[a,b] = [1.570796326792186,1.5707963267990062],f((a+b)/2) = 1429839342339.8782
it = 44,[a,b] = [1.570796326792186,1.570796326795596],f((a+b)/2) = -994550654629.8223
it = 45,[a,b] = [1.5707963267938911,1.570796326795596],f((a+b)/2) = -6533813968064.459
it = 46,[a,b] = [1.5707963267947436,1.570796326795596],f((a+b)/2) = 3662284195633.0815
it = 47,[a,b] = [1.5707963267947436,1.5707963267951697],f((a+b)/2) = 16697052243463.29
it = 48,[a,b] = [1.5707963267947436,1.5707963267949565],f((a+b)/2) = -21417587553013.977
it = 49,[a,b] = [1.57079632679485,1.5707963267949565],f((a+b)/2) = 151512722541384.88
it = 50,[a,b] = [1.57079632679485,1.5707963267949032],f((a+b)/2) = -49887138373250.875
it = 51,[a,b] = [1.5707963267948766,1.5707963267949032],f((a+b)/2) = -148752623989352.16
it = 52,[a,b] = [1.57079632679489,1.5707963267949032],f((a+b)/2) = -1.6331239353195368e+16
最后
Python崛起并且风靡,因为优点多、应用领域广、被大牛们认可。学习 Python 门槛很低,但它的晋级路线很多,通过它你能进入机器学习、数据挖掘、大数据,CS等更加高级的领域。Python可以做网络应用,可以做科学计算,数据分析,可以做网络爬虫,可以做机器学习、自然语言处理、可以写游戏、可以做桌面应用…Python可以做的很多,你需要学好基础,再选择明确的方向。这里给大家分享一份全套的 Python 学习资料,给那些想学习 Python 的小伙伴们一点帮助!
👉Python所有方向的学习路线👈
Python所有方向的技术点做的整理,形成各个领域的知识点汇总,它的用处就在于,你可以按照上面的知识点去找对应的学习资源,保证自己学得较为全面。
👉Python必备开发工具👈
工欲善其事必先利其器。学习Python常用的开发软件都在这里了,给大家节省了很多时间。
👉Python全套学习视频👈
我们在看视频学习的时候,不能光动眼动脑不动手,比较科学的学习方法是在理解之后运用它们,这时候练手项目就很适合了。
👉实战案例👈
学python就与学数学一样,是不能只看书不做题的,直接看步骤和答案会让人误以为自己全都掌握了,但是碰到生题的时候还是会一筹莫展。
因此在学习python的过程中一定要记得多动手写代码,教程只需要看一两遍即可。
👉大厂面试真题👈
我们学习Python必然是为了找到高薪的工作,下面这些面试题是来自阿里、腾讯、字节等一线互联网大厂最新的面试资料,并且有阿里大佬给出了权威的解答,刷完这一套面试资料相信大家都能找到满意的工作。
一个人可以走的很快,但一群人才能走的更远!不论你是正从事IT行业的老鸟或是对IT行业感兴趣的新人,都欢迎加入我们的的圈子(技术交流、学习资源、职场吐槽、大厂内推、面试辅导),让我们一起学习成长!