C
C       ESGL.FT
C       CALCULATE E USING AN INFINITE SERIES - DOUBLE PRECISION
C       CONTINUE ITERATING UNTIL THE CALCULATED VALUE DOES NOT CHANGE

        DOUBLE PRECISION E, E0

        I = 0
        E = 1.0

10      E0 = E
        I = I + 1 
        D = 1.0
        DO 20 J = 1, I
        D = D * FLOAT( J )
20      CONTINUE
        E = E + ( 1.0 / D )
        WRITE (4,30) I, E
30      FORMAT ( ' I = ',I3,' E = ', D36.30 )
        IF ( E .NE. E0 ) GOTO 10

40      WRITE (4,50)
50      FORMAT ( ' ALL DONE ' )
        STOP
        END