fork download
copy 
  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. X2 = min(a/2.0,b/2.0)
  8.  
  9. DO WHILE ( ABS(X2-X1) > d)
  10.   h = ABS(X2 - X1)/3
  11.   F1 = (a - 2*(X1+h))* (b - 2*(X1+h))*(X1+h)
  12.   F2 = (a-2*(X2-h))*(b-2*(X2-h))*(X2-h)
  13.  IF ( F1 <= F2 ) THEN
  14.   X1=X1 + h 
  15.  ELSE 
  16.   X2=X2 - h
  17. END IF
  18. Xmax = (X1 + X2)/2
  19. PRINT *, '体積最大を得るｘ(暫定)は',Xmax
  20. PRINT *, 'その時の体積は', (a-2*(Xmax))*(b-2*(Xmax))*(Xmax)
  21. END DO
  22. PRINT *, '体積最大を得るx(最終結果)は',Xmax
  23. PRINT *, '体積最大は', (a-2*(Xmax))*(b-2*(Xmax))*(Xmax)
  24.  
  25. END PROGRAM ex5_9
  26.  
Success #stdin #stdout 0s 5296KB
comments (?)
stdin
copy 
100,175
stdout
copy 
 体積最大を得るx(最終結果)は   0.00000000E+00
 体積最大は   0.00000000E+00
https://ideone.com/1R0XuL
language:
Fortran (gfortran 8.3)
created:
3 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!