test `sol.converged` after `sol = fsolve`. check: https://people.sc.fsu.edu/~jburkardt/f_src/fsolve/fsolve.f90