SUBROUTINE MATA

PURPOSE
   PREMULTIPLY A MATRIX BY ITS TRANSPOSE TO FORM A
   SYMMETRIC MATRIX

USAGE
   CALL MATA(A,R,N,M,MS)

DESCRIPTION OF PARAMETERS
   A  - NAME OF INPUT MATRIX
   R  - NAME OF OUTPUT MATRIX
   N  - NUMBER OF ROWS IN A
   M  - NUMBER OF COLUMNS IN A. ALSO NUMBER OF ROWS AND
	NUMBER OF COLUMNS OF R.
   MS  - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A
	  0 - GENERAL
	  1 - SYMMETRIC
	  2 - DIAGONAL

REMARKS
   MATRIX R CANNOT BE IN THE SAME LOCATION AS MATRIX A
   MATRIX R IS ALWAYS A SYMMETRIC MATRIX WITH A STORAGE MODE=1

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
   LOC

METHOD
   CALCULATION OF (A TRANSPOSE A) RESULTS IN A SYMMETRIC MATRIX
   REGARDLESS OF THE STORAGE MODE OF THE INPUT MATRIX. THE
   ELEMENTS OF MATRIX A ARE NOT CHANGED.