fork download 
  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
Success #stdin #stdout 0s 5320KB
comments (?)
stdin 
100 167
stdout 
   16.6666679       148518.516    
   11.1111126       125116.594    
   14.8148165       143211.906    
   17.2839527       149770.812    
   18.9300442       151908.734    
   20.0274372       152403.625    
   20.7590332       152338.641    
   21.2467651       152125.297    
   20.9216118       152282.453    
   20.7048416       152354.016    
   20.5603294       152386.844    
   20.4639854       152402.062    
   20.3997574       152409.234    
   20.3569412       152412.703    
   20.3283958       152414.406    
   20.3093643       152415.297    
   20.2966766       152415.766    
   20.2882195       152416.016    
   20.2825813       152416.172    
   20.2788200       152416.266    
   20.2763138       152416.312    
   20.2746429       152416.359
https://ideone.com/9AHnFq
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!