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 *, '体積最大を得るx(最終結果)は',Xmax
  24. PRINT *, '体積最大は', (a-2*(Xmax))*(b-2*(Xmax))*(Xmax)
  25.  
  26. END PROGRAM ex5_9
  27.  
Success #stdin #stdout 0s 5316KB
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    
 体積最大を得るx(最終結果)は   20.4899788    
 体積最大は   162073.047
https://ideone.com/ScOOjU
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!