Code (filename hello.f):
PROGRAM BISECTION DOUBLE PRECISION A, B, C, FA, FB, FC, TOL EXTERNAL F TOL = 1.0D-6 A = 0.0D0 B = 2.0D0 IF (F(A)*F(B) .GE. 0.0) STOP 'No sign change' DO 10 I = 1, 100 C = (A + B) / 2.0D0 FC = F(C) IF (FC .EQ. 0.0 .OR. (B-A)/2.0 .LT. TOL) THEN WRITE(*,*) 'Root: ', C STOP ENDIF FA = F(A) IF (FA*FC .LT. 0.0) THEN B = C ELSE A = C ENDIF 10 CONTINUE END DOUBLE PRECISION FUNCTION F(X) DOUBLE PRECISION X F = X**3 - X**2 - 2.0D0 RETURN END fortran 77 compiler for windows 10 64-bit free download
By using the modern, open-source GFortran compiler, you breathe new life into legacy code while maintaining the speed and stability of a 64-bit Windows 10 system. Code (filename hello