Function:      fmatrix

Purpose:       Allocate space for a matrix with subscript range
               M[ylo..yhi][xlo..xhi] for floats.
               
Category:      UTILITY

File:          matrix.c

Author:        M. Vogelaar 

Use:           float **fmatrix( int xlo, int ylo, int xhi, int yhi );

                        
               fmatrix : Output pointer to array of pointers to array
                         of floats.
                         e.g. the matrix M[ylo..yhi][xlo..xhi]
               xlo     : First subscript in x
               ylo     : First subscript in y               
               xhi     : Last subscript in x
               yhi     : Last subscript in y
                 
                
Example:       float      **image;
 
               image = fmatrix( blo[0], blo[1], bhi[0], bhi[1] );
               x = 3; y = 0;
               element = M[y][x];
               freefmatrix( blo[0], blo[1] );

Comment:       This routine is NOT callable from FORTRAN.

Notes:         After freeing memory with 'freefmatrix' the pointer to
               the matrix is NOT reset to NULL.
               
Updates:       28 Jul,  2000: VOG, Document created.