/* Copyright (C) 1986 Adam Fritz, 133 Main St., Afton, N.Y. 13730 */

MATGEN ( A, lda, n )
   float A[][ldx] ;
   int lda, n ;

     /* Program:                                            */
     /*                                                     */
     /*    MATGEN                                           */
     /*                                                     */
     /* Version:                       Date:                */
     /*                                                     */
     /*    C                              04/25/86          */
     /*                                                     */
     /* Description:                                        */
     /*                                                     */
     /*    Generate a test matrix using Hilbert co-         */
     /*    efficients.  The matrix is symmetric and         */
     /*    poorly conditioned.                              */
     /*                                                     */
     /* Author:                                             */
     /*                                                     */
     /*    Adam Fritz                                       */
     /*    133 Main Street                                  */
     /*    Afton, New York 13730                            */

{
   int i, j, ip1 ;
                                /* Validate Leading Dimension */
   if ( lda > 0 )
   {
                                /* Validate Order */
      if ((n > 1) && (n <= lda))
      {
                                /* Form Matrix */
         for ( i=0; i<n; i++ )
         {
            A[i][i] = 1.0/((float) 2*i+1) ;
            ip1 = i + 1 ;
            for ( j=ip1; j<n; j++ )
               A[i][j] = A[j][i] = 1.0/((float) i+j+1) ;
         }
      }
      else
      {
         printf("Error: Invalid MATRIX Order, %d\n", n) ;
         bios(0) ;
      }
   }
   else
   {
      printf("Error: Invalid LEADING DIMENSION, %d\n", lda) ;
      bios(0) ;
   }
}

/* Copyright (C) 1986 Adam Fritz, 133 Main St., Afton, N.Y. 13730 */