fork download
copy 
  1. PROGRAM BISECTION_1_1
  2. IMPLICIT NONE
  3. REAL(8) X,A,AI,B,BI,C,CI
  4. REAL(8) :: EPS=1.0D-10
  5. REAL(8) :: e=2.71828182845904
  6. INTEGER I
  7.  
  8. A=10.0D0
  9. B=-2.0D0
  10. C=(A+B)/2
  11.  
  12. WRITE(*,*) 0,C
  13.  
  14. DO I=1,100
  15.  AI=A
  16.  BI=B
  17.  CI=(AI+BI)/2
  18.   IF(F(AI)*F(CI)<=0) THEN
  19.   B=CI
  20.   ELSE
  21.   A=CI
  22.  ENDIF
  23.  WRITE(*,*) I,CI
  24.  IF(ABS(B-A)<EPS) EXIT
  25. ENDDO
  26. STOP
  27.  
  28. CONTAINS
  29.  
  30. REAL(8) FUNCTION F(X)
  31.   REAL(8),INTENT(IN):: X
  32.   F=e**X-1.0D0
  33.   RETURN
  34. END FUNCTION F
  35.  
  36. END PROGRAM BISECTION_1_1
Success #stdin #stdout 0.01s 5428KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
           0   4.0000000000000000     
           1   4.0000000000000000     
           2   1.0000000000000000     
           3 -0.50000000000000000     
           4  0.25000000000000000     
           5 -0.12500000000000000     
           6   6.2500000000000000E-002
           7  -3.1250000000000000E-002
           8   1.5625000000000000E-002
           9  -7.8125000000000000E-003
          10   3.9062500000000000E-003
          11  -1.9531250000000000E-003
          12   9.7656250000000000E-004
          13  -4.8828125000000000E-004
          14   2.4414062500000000E-004
          15  -1.2207031250000000E-004
          16   6.1035156250000000E-005
          17  -3.0517578125000000E-005
          18   1.5258789062500000E-005
          19  -7.6293945312500000E-006
          20   3.8146972656250000E-006
          21  -1.9073486328125000E-006
          22   9.5367431640625000E-007
          23  -4.7683715820312500E-007
          24   2.3841857910156250E-007
          25  -1.1920928955078125E-007
          26   5.9604644775390625E-008
          27  -2.9802322387695312E-008
          28   1.4901161193847656E-008
          29  -7.4505805969238281E-009
          30   3.7252902984619141E-009
          31  -1.8626451492309570E-009
          32   9.3132257461547852E-010
          33  -4.6566128730773926E-010
          34   2.3283064365386963E-010
          35  -1.1641532182693481E-010
          36   5.8207660913467407E-011
          37  -2.9103830456733704E-011
https://ideone.com/e1mH7t
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!