SUBROUTINE INIT(X, LDX, M, N)
C
C     PURPOSE:
C
C     The subroutine INIT initializes an M by N matrix X with the M by N
C     identity matrix, characterized by unit diagonal entries and zero
C     off-diagonal elements.
C
C     ARGUMENT LIST:
C
C     X - DOUBLE PRECISION array of DIMENSION (LDX,N)
C         On return, X contains the M by N identity matrix.
C     LDX - INTEGER
C         LDX is the leading dimension of the array X (LDX >= M).
C     M - INTEGER
C         M is the number of rows of the matrix X.
C     N - INTEGER
C         N is the number of columns of the matrix X.
C
C     CONTRIBUTOR: S. Van Huffel, (ESAT Laboratory, KU Leuven).
C
C     REVISIONS: 1988, February 15.
C
C     .. Scalar Arguments ..
      INTEGER LDX, M, N
C     .. Array Arguments ..
      DOUBLE PRECISION X(LDX,*)
C     .. Local Scalars ..
      INTEGER I, J
C     .. Executable Statements ..
C
      DO 20 J = 1, N
         DO 10 I = 1, M
            X(I,J) = 0.0D0
   10    CONTINUE
         X(J,J) = 1.0D0
   20 CONTINUE
      RETURN
C *** Last line of INIT ***********************************************
      END