fork download
copy 
  1. PROGRAM ex5_9
  2.  
  3.  IMPLICIT NONE
  4.  REAL :: a, b, X1=0, X2, d, h, F1, F2, Xmax
  5.  d=0.0001
  6.   READ *, a, b
  7.  
  8. X2 = min(a/2.0,b/2.0)
  9.  
  10. DO WHILE ( ABS(X2-X1) > d)
  11.   h = ABS(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
  21. PRINT *,  (a-2*(Xmax))*(b-2*(Xmax))*(Xmax)
  22. END DO
  23. PRINT *, Xmax
  24. PRINT *,  (a-2*(Xmax))*(b-2*(Xmax))*(Xmax)
  25.  
  26. END PROGRAM ex5_9
  27.  
Success #stdin #stdout 0s 5320KB
comments (?)
stdin
copy 
100,175
stdout
copy 
   16.6666679    
   157407.406    
   22.2222233    
   161179.703    
   18.5185204    
   160862.156    
   20.9876556    
   161997.766    
   19.3415642    
   161666.922    
   20.4389591    
   162072.297    
   21.1705551    
   161932.828    
   20.6828232    
   162061.578    
   20.3576698    
   162067.844    
   20.5744400    
   162070.797    
   20.4299278    
   162071.984    
   20.5262680    
   162072.609    
   20.4620399    
   162072.828    
   20.5048599    
   162072.953    
   20.4763145    
   162072.984    
   20.4953442    
   162073.016    
   20.4826584    
   162073.031    
   20.4911156    
   162073.031    
   20.4854774    
   162073.047    
   20.4892349    
   162073.047    
   20.4867287    
   162073.047    
   20.4883995    
   162073.031    
   20.4895134    
   162073.031    
   20.4902573    
   162073.047    
   20.4897614    
   162073.047    
   20.4900913    
   162073.047    
   20.4898720    
   162073.047    
   20.4897232    
   162073.047    
   20.4898224    
   162073.047    
   20.4898872    
   162073.047    
   20.4899292    
   162073.031    
   20.4899597    
   162073.047    
   20.4899788    
   162073.047    
   20.4899788    
   162073.047
https://ideone.com/yYQ3iw
language:
Fortran (gfortran 8.3)
created:
16 hours 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!