      PROGRAMME MAIN
C     From - Wed Aug 28 18:25:35 2002
C     Spaltenzaehlung in den folgenden zwei  Zeilen
C2345678901234567890123456789012345678901234567890123456789012345678901234567890
C        |10       |20       |30       |40       |50       |60         |72
C ********************************************************************
C     Bildschirm- Ein- und Ausgabe 
C     Iterative Annaeherung an die Kreiszahl pi
C     pi/4 = 1/1 -1/3 +1/5 -1/7 ... (konvergiert sehr langsam)
C
C     Autor:  Wolfgang Loebnitz  (Wed Aug 28 18:25:35 2002)
C ********************************************************************
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      PI=3.1415
      PIEXAKT=4*ATAN(1.0)
      WRITE(*,100)
      READ(*,110) I
      PV = 0.D0
      V  = -1.D0
      II = 2*I-1
      DO 120 J=1,II,2
      V  = -V
      Q  = V/J
      PV = PV+Q
C     WRITE(*,290) J, Q, PV
  120 CONTINUE	
      PI = PV*4 
      FEHLER=PIEXAKT-PI
      WRITE(*,300) I
      WRITE(*,310) PI
      WRITE(*,310) PIEXAKT
      WRITE(*,310) FEHLER
 9999 CONTINUE		
  100 FORMAT("Anzahl der Iterationen zur Berechnung von pi?")
  110 FORMAT(I8)
  290 FORMAT(I8,2X,F12.5,2X,F12.5)
  300 FORMAT("Kreiszahl pi nach",2X,I8,1X," Iterationen:")
  310 FORMAT(1X,F18.15)
C ---------------------------------------------------------------------
      STOP
      END
