REAL DATA(18,16)
        REAL ARRAY1(16,16), ARRAY2(16,16), VECT1(16), VECT2(16)
        real vect3(18), vect4(18)
C
	OPEN(UNIT=21,STATUS='OLD',FILE='spectr.dat')
C
	N = 18
	M = 16
	DO I = 1, N
	   READ(21,100)(DATA(I,J),J=1,M)
  100	   FORMAT(8F7.1)
 	ENDDO
C	CALL OUTMAT(N,M,DATA)
C
	METHOD = 3
	IPRINT = 3
	CALL PCA(N,M,DATA,METHOD,IPRINT,ARRAY1,VECT1,VECT2,
     X              vect3,vect4,ARRAY2,IERR)
             
	IF (IERR.NE.0) GOTO 9000
C
	GOTO 9900
 9000	WRITE (6,*) ' ABNORMAL END: IERR =', IERR
 9900	CONTINUE
	END