PROGRAM D12R10
C	Driver for routine EVLMEM
	PARAMETER(N=1000,M=10,NFDT=16)
	DIMENSION DATA(N),COF(M),WK1(N),WK2(N),WKM(M)
	OPEN(5,FILE='SPCTRL.DAT',STATUS='OLD')
	READ(5,*) (DATA(I),I=1,N)
	CLOSE(5)
	CALL MEMCOF(DATA,N,M,PM,COF,WK1,WK2,WKM)
	WRITE(*,*) 'Power spectrum estimate of data in SPCTRL.DAT'
	WRITE(*,'(1X,T6,A,T20,A)') 'f*delta','power'
	DO 11 I=0,NFDT
		FDT=0.5*I/NFDT
		WRITE(*,'(1X,2F12.6)') FDT,EVLMEM(FDT,COF,M,PM)
11	CONTINUE
	END