Program for bisection method in fortran 77