fork download
copy 
  1. PROGRAM ex14
  2. REAL::h, FX1, FX2, X1, X2, X3, X4, Xmax
  3. READ *,a, b
  4.  X1=0
  5.  X2=MIN(a, b)/2
  6. DO WHILE(X2-X1>0.01)
  7.  FX1=(a-2*X3)*(b-2*X3)*X3
  8.  FX2=(a-2*X4)*(b-2*X4)*X4
  9.  h=(X2-X1)/3
  10.  X3=X1+h
  11.  X4=X2-h
  12. IF(FX1<FX2) THEN
  13.  X1=X3
  14. ELSE
  15.  X2=X4
  16. END IF
  17. Xmax=(X1+X2)/2
  18. PRINT *, Xmax, (a-2*Xmax)*(b-2*Xmax)*Xmax
  19. END DO
  20.  
  21. END PROGRAM ex14
  22.  
Success #stdin #stdout 0.01s 5324KB
comments (?)
stdin
copy 
100 567
stdout
copy 
   16.6666679       592963.000    
   11.1111126       470795.625    
   14.8148165       560221.562    
   17.2839527       602140.938    
   18.9300442       622433.250    
   20.0274372       632622.500    
   20.7590332       637950.000    
   21.2467651       640855.875    
   21.5719185       642507.688    
   21.7886868       643482.625    
   21.9332008       644076.562    
   22.0295410       644447.688    
   22.0937691       644684.000    
   22.1365891       644836.750    
   22.1651344       644936.375    
   22.1841640       645001.750    
   22.1968498       645045.000    
   22.2053089       645073.562    
   22.2109489       645092.562    
   22.2147064       645105.125    
   22.2172127       645113.500    
   22.2188835       645119.188
https://ideone.com/37bFUi
language:
Fortran (gfortran 8.3)
created:
7 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!