SUBROUTINE MFUN

PURPOSE
   APPLY A FUNCTION TO EACH ELEMENT OF A MATRIX TO FORM A
   RESULTANT MATRIX

USAGE
   CALL MFUN (A,F,R,N,M,MS)
   AN EXTERNAL STATEMENT MUST PRECEDE CALL STATEMENT IN ORDER
   TO IDENTIFY PARAMETER F AS THE NAME OF A FUNCTION

DESCRIPTION OF PARAMETERS
   A - NAME OF INPUT MATRIX
   F - NAME OF FORTRAN-FURNISHED OR USER FUNCTION SUBPROGRAM
   R - NAME OF OUTPUT MATRIX
   N - NUMBER OF ROWS IN MATRIX A AND R
   M - NUMBER OF COLUMNS IN MATRIX A AND R
   MS  - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A (AND R)
	  0 - GENERAL
	  1 - SYMMETRIC
	  2 - DIAGONAL

REMARKS
   PRECISION IS DEPENDENT UPON PRECISION OF FUNCTION USED

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   LOC

METHOD
   FUNCTION F IS APPLIED TO EACH ELEMENT OF MATRIX A
   TO FORM MATRIX R