fork download
copy 
  1. PROGRAM ex13_1
  2. IMPLICIT NONE
  3. INTEGER :: a, b
  4. REAL :: X1, X2, Xmax, eps, h, x,F1, F2
  5. READ *, a, b
  6. X1 = 0.0
  7. X2 = MIN(a, b) / 2.0
  8. eps = 0.01
  9. h = (X2 - X1) / 3.0
  10. Xmax = (X1 + X2) / 2.0
  11. DO
  12. IF (ABS(X2 - X1) < eps) EXIT
  13. F1 = (a - 2.0*(X1+h))*(b - 2.0*(X1+h))*(X1+h)
  14. F2 = (a - 2.0*(X2-h))*(b - 2.0*(X2-h))*(X2-h)
  15. IF (F1 < F2) THEN
  16. X1 = X1 + h
  17. ELSE
  18. X2 = X2 - h
  19. END IF
  20. h = (X2 - X1) / 3.0
  21. Xmax = (X1 + X2) / 2.0
  22. PRINT *, Xmax, (a - 2.0*Xmax)*(b - 2.0*Xmax)*Xmax
  23. END DO
  24. END PROGRAM ex13_1
Success #stdin #stdout 0s 5320KB
comments (?)
stdin
copy 
100,136
stdout
copy 
   16.6666679       114074.070    
   22.2222233       113031.547    
   18.5185204       115388.922    
   20.9876556       114503.477    
   19.3415642       115414.328    
   18.2441711       115305.906    
   18.9757671       115443.953    
   19.4634972       115390.125    
   19.1383438       115438.789    
   18.9215755       115442.805    
   19.0660896       115442.672    
   18.9697456       115443.898    
   19.0339737       115443.578    
   18.9911575       115444.023    
   19.0197029       115443.828    
   19.0006714       115444.000    
   18.9879837       115444.008    
   18.9964428       115444.008    
   18.9908028       115444.016    
   18.9945641       115444.016    
   18.9920578       115444.023    
   18.9903870       115444.016
https://ideone.com/xAmkMp
language:
Fortran (gfortran 8.3)
created:
2 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!