fork download
  1. PROGRAM ex5_9
  2.  
  3. IMPLICIT NONE
  4. REAL :: a, b, X1=0, X2, d, h, F1, F2, Xmax
  5. READ *, a, b
  6.  
  7.  
  8. X2 = min(a, b)/2
  9.  
  10. DO WHILE ( ABS(X2-X1) > d)
  11. h = (X2 - X1)/3
  12. F1 = (a - 2*(X1+h))* (b - 2*(X1+h))*(X1+h)
  13. F2 = (a-2*(X2+h))*(b-2*(X2+h))*(X2+h)
  14. IF ( F1 < F2 ) THEN
  15. X1=X1 + h
  16. ELSE
  17. X2=X2 - h
  18. END IF
  19. Xmax = (X1 + X2)/2
  20. PRINT *, " Xmax, (a-2*(Xmax+h))*(b-2*(Xmax+h))*(Xmax+h) "
  21. END DO
  22.  
  23. END PROGRAM ex5_9
  24.  
Success #stdin #stdout 0s 5288KB
stdin
100,175
stdout
Standard output is empty