Program for bisection method in fortran