我附上了方程式的图片,我应该使用牛顿法找到其根源。我的代码给了我很大的扎根,应该接近7。
from math import *
def newton(x):
fx= (((3*10**-5)*exp((log1p(10**3/(3*10**-5)))*(1-exp(-.12*x))))-1)
dfx= (.12*(3*10**-5)*exp((log1p(10**3/(3*10**-5)))*(1-exp(-.12*x))-.12*x))
tolerance= .000001
x0=100
xnew = x0
xold= x0-1
iterations =0
while (abs(xnew - xold) > tolerance or iterations>1000):
xold= xnew
xnew= xnew- (fx/dfx)
iterations= iterations +1
return abs(xnew)