Apply fixed‑point iteration (#587)

Apply fixed‑point iteration (#587)

Rewrite equation as x=g(x)x = g(x) and iterate xn+1=g(xn)x_{n+1} = g(x_n).