Newton-Raphson Code

logicle_nr.py

def newton_method(f, df, x0, err, W):
    """
    Estimate p using Newton-Raphson method.
    
    Output: root approximation given error criterion (err)
    """
    delta = err+1
    x = x0
    while delta > err:
        x1 = x - f(x, W) / df(x)
        delta = abs(x1 - x)
        x = x1
    return x

def newton_method2(f, df, x0, err, W, S, T, M, p):
    """
    Estimate X using Newton-Raphson method.
    
    Output: root approximation given error criterion (err)
    """
    delta = err+1
    x = x0
    while delta > err:
        x1 = x - f(S, x, p, W, M, T) / df(x, p, W, M, T)
        delta = abs(x1 - x)
        x = x1
    return x