涉及牛顿法的寻根问题

时间:2019-02-13 03:11:33

标签: python python-3.x

enter image description here

enter image description here

我附上了方程式的图片,我应该使用牛顿法找到其根源。我的代码给了我很大的扎根,应该接近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)

0 个答案:

没有答案