假设a>0，利用迭代法求方程的解：x²-a=0，给出代码

当你需要使用迭代法求解像 \( x^2 - a = 0 \) 这样的方程时，牛顿法（Newton-Raphson method）是一个常用的数值方法。假设我们知道函数 \( f(x) = x^2 - a \)，零点就是我们需要找到的解。以下是使用Python实现的一个简单版本：

def newton_method(a, initial_guess=1, tolerance=1e-6, max_iterations=100):
    def func_and_derivative(x):
        return x**2 - a, 2 * x

    x = initial_guess
    for _ in range(max_iterations):
        fx, dx = func_and_derivative(x)
        if abs(fx) < tolerance:
            break
        x -= fx / dx
    else:
        print("警告：达到最大迭代次数，结果可能不准确")
        
    return x

# 示例：求解当a=5时的解
solution = newton_method(5)
print(f"解大约是 {solution}")