SUBROUTINE UNIF ( IX, U )
        INTEGER K1, IX
        REAL U
        K1 = IX / 127773
        IX = 16807 * ( IX - K1 * 127773 ) - K1 * 2836
        IF ( IX .LT. 0 ) IX = IX + 2147483647
        U = IX * 4.656612875e-10
        END

        PROGRAM Simulation
        REAL U
        INTEGER IX, I
        WRITE (*,*) ' Initialisez le generateur'
        READ (*,*) IX
        DO 10 I = 1, 20
        CALL UNIF ( IX, U )
        WRITE (*,*) IX, U
10      CONTINUE
        END