Newton Method to Calculate Root Python

What is the Newton Method?

The Newton-Raphson method, also known as Newton's method, is an iterative numerical technique used to find successively better approximations to the roots of a real-valued function. It's particularly useful when analytical solutions are difficult or impossible to find.

The method works by:

1. Starting with an initial guess for the root
2. Using the function and its derivative to estimate a better approximation
3. Repeating the process until the approximation is sufficiently accurate

The Newton-Raphson method typically converges quadratically when close to the root, making it very efficient for many problems.

Python Implementation

Here's a Python function that implements the Newton-Raphson method:

def newton_raphson(f, df, x0, tol=1e-6, max_iter=100):
    """
    Newton-Raphson method to find a root of f(x) = 0.

    Parameters:
    f: function - the function for which to find a root
    df: function - the derivative of f
    x0: float - initial guess
    tol: float - tolerance for stopping criterion
    max_iter: int - maximum number of iterations

    Returns:
    float - the approximate root
    """
    x = x0
    for i in range(max_iter):
        fx = f(x)
        if abs(fx) < tol:
            return x
        dfx = df(x)
        if dfx == 0:
            raise ValueError("Derivative is zero. No solution found.")
        x = x - fx/dfx
    raise ValueError("Maximum iterations reached. No solution found.")

To use this function, you need to define your function and its derivative. For example, to find the square root of 2:

def f(x):
    return x**2 - 2

def df(x):
    return 2*x

root = newton_raphson(f, df, x0=1.5)
print(f"Approximate root: {root}")

This implementation includes:

• Error handling for division by zero in the derivative
• A maximum iteration limit to prevent infinite loops
• A tolerance parameter to control the accuracy of the result

Example Calculation

Let's find the square root of 2 using the Newton-Raphson method.

After just 2 iterations, we've found that √2 ≈ 1.4142, which is accurate to 4 decimal places.