fork download
copy 
  1. PROGRAM BISECTION_1_2
  2. IMPLICIT NONE
  3. REAL(8) X,A,AI,B,BI,C,CI
  4. REAL(8) :: EPS=1.0D-10
  5. INTEGER I
  6.  
  7. A=18.0D0
  8. B=1.0D0
  9. C=(A+B)/2
  10.  
  11. WRITE(*,*) 0,C
  12.  
  13. DO I=1,100
  14.  AI=A
  15.  BI=B
  16.  CI=(AI+BI)/2
  17.   IF(F(AI)*F(CI)<=0) THEN
  18.   B=CI
  19.   ELSE
  20.   A=CI
  21.  ENDIF
  22.  WRITE(*,*) I,CI
  23.  IF(ABS(B-A)<EPS) EXIT
  24. ENDDO
  25. STOP
  26.  
  27. CONTAINS
  28.  
  29. REAL(8) FUNCTION F(X)
  30.   REAL(8),INTENT(IN):: X
  31.   F=X*log(X)-1.0D0
  32.   RETURN
  33. END FUNCTION F
  34.  
  35. END PROGRAM BISECTION_1_2
Success #stdin #stdout 0.01s 5432KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
           0   9.5000000000000000     
           1   9.5000000000000000     
           2   5.2500000000000000     
           3   3.1250000000000000     
           4   2.0625000000000000     
           5   1.5312500000000000     
           6   1.7968750000000000     
           7   1.6640625000000000     
           8   1.7304687500000000     
           9   1.7636718750000000     
          10   1.7470703125000000     
          11   1.7553710937500000     
          12   1.7595214843750000     
          13   1.7615966796875000     
          14   1.7626342773437500     
          15   1.7631530761718750     
          16   1.7634124755859375     
          17   1.7632827758789062     
          18   1.7632179260253906     
          19   1.7632503509521484     
          20   1.7632341384887695     
          21   1.7632260322570801     
          22   1.7632219791412354     
          23   1.7632240056991577     
          24   1.7632229924201965     
          25   1.7632224857807159     
          26   1.7632227391004562     
          27   1.7632228657603264     
          28   1.7632228024303913     
          29   1.7632228340953588     
          30   1.7632228499278426     
          31   1.7632228420116007     
          32   1.7632228380534798     
          33   1.7632228360744193     
          34   1.7632228350848891     
          35   1.7632228345901240     
          36   1.7632228343427414     
          37   1.7632228344664327     
          38   1.7632228344045870
https://ideone.com/u1Nm7t
language:
Fortran (gfortran 8.3)
created:
2 years 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!