COPYRIGHT (C) 1986 ADAM FRITZ, 133 MAIN ST., AFTON, N.Y. 13730
C
C     PROGRAM: MATGEN
C
C     VERSION: 1.0                            DATE: 04/25/86
C              MICROSOFT FORTRAN
C
C     DESCRIPTION:
C
C        GENERATE A TEST MATRIX USING HILBERT COEFFICIENTS.
C
C     AUTHOR:
C
C        ADAM FRITZ
C        133 MAIN STREET
C        AFTON, NEW YORK 13730
C
      SUBROUTINE MATGEN(A, LDA, N)
C
      INTEGER LDA, N, I, J, NM1, IP1
      REAL A(LDA,1), T, AIJ
C
      IF (N .LT. 1 .OR. N .GT. LDA) STOP
      NM1 = N - 1
      IF (N .LT. 2) GO TO 30
         DO 20 I=1, NM1
            T = 1.E0/FLOAT(2*I-1)
            A(I,I) = T
            IP1 = I+1
            DO 10 J=IP1, N
               AIJ = 1.E0/FLOAT(I+J-1)
               A(I,J) = AIJ
               A(J,I) = AIJ
   10       CONTINUE
   20    CONTINUE
   30 CONTINUE
      T = 1.E0/FLOAT(2*N-1)
      A(N,N) = T
C
      RETURN
      END
C
COPYRIGHT (C) 1986 ADAM FRITZ, 133 MAIN ST., AFTON, N.Y. 13730