gfsolver Hybrid Newton-bisection solver in Fortran It finds the zeros of a function using the basic Newton method and switches to hybrid Newton-bisection if the max number of iterations is exceeded.