import numpy as np
from scipy.optimize import root_scalar

def equation(x):
    return np.cos(x) - x

solution = root_scalar(equation, bracket=[0, 1], method="brentq", xtol=1e-12)

print(f"converged: {solution.converged}")
print(f"flag: {solution.flag}")
print(f"root: {solution.root:.12f}")
print(f"residual: {equation(solution.root):.3e}")
print(f"iterations: {solution.iterations}")
print(f"function calls: {solution.function_calls}")