fork download
copy 
  1. PROGRAM ex5_7
  2.  
  3.   IMPLICIT NONE
  4.   REAL :: r = 0.0E0, x, d = 1.0E-5
  5.  
  6.   DO
  7.     PRINT *, 'Input a positive real number: '
  8.     READ *, r
  9.     IF ( r > 0.0E0 ) EXIT
  10.   END DO
  11.  
  12.   x = r                             ! 初期値は0以外なら何でもよい
  13.  
  14.   DO WHILE ( ABS( x ** 2 - r ) / r > d )
  15.     x = 0.5E0 * ( x + r / x )
  16.     PRINT '( 1X, E12.5 )', x
  17.   END DO
  18.  
  19. END PROGRAM ex5_7
Success #stdin #stdout 0.01s 5288KB
comments (?)
stdin
copy 
144
stdout
copy 
 Input a positive real number: 
  0.72500E+02
  0.37243E+02
  0.20555E+02
  0.13780E+02
  0.12115E+02
  0.12001E+02
  0.12000E+02
https://ideone.com/PjiceM
language:
Fortran (gfortran 8.3)
created:
10 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!