【计算方法】二分法求方程数值解 python实现

#[a,b]为区间,e是误差限,k为计算的轮次,x为得到的数值解
import math
def dichotomize(a, b, e):
	k=0
	while(abs(a-b)/2>=e):
		k=k+1
		x = (a+b)/2
		if func(x)<0:
			a = x
		elif func(x)>0:
			b = x
		else:
			break
	return x,k

#求解的方程
 def func(x):
	y = pow(x, 3) - x -1
	return y

极限学习工程计算,希望我考试能过!
算法自己实现一遍就很开心哈哈。

你可能感兴趣的:(python,算法)