Program for bisection method in fortran code do loop