File: EX1POW.FT of Tape: Various/ETH/eth11-3

(Source file text) 

---

	DIMENSION PSPECT(32),IMAG(32)
	INTEGER REAL(32),SCALE,ERROR
	DATA N,I1,I2,T,A,B,C,D/32,0,0,0.,100.,200.,.5,2./
	DATA PI/3.141593/,IMAG/32*0/

	DT=PI/8.
	CALL RANDU(I1,I2,X)
	DO 1 I=1,N
	REAL(I)=A*SIN(C*T)+B*COS(D*T)+(-1.**I)*15.*RAN(I1,I2)
    1	T=T+DT
	TYPE 1000
	TYPE 1001,REAL
	TYPE 1002,IMAG

	CALL FFT(ERROR,N,REAL,IMAG,0,SCALE)
	IF(ERROR.NE.0) PRINT 3000,ERROR

	CALL POWRSP(N,REAL,IMAG,PSPECT)
	TYPE 2000
	TYPE 1001,REAL
	TYPE 1002,IMAG
	TYPE 2001,PSPECT

	STOP

 1000	FORMAT(' INPUT DATA FOR FFT')
 1001	FORMAT(/' ***REAL PART***',/,(4I16))
 1002	FORMAT(' ***IMAGINARY PART***',/,(4I16))
 2000	FORMAT(//' RESULTS OF THE FFT')
 2001	FORMAT(//' POWER SPECTRUM RESULTS',/,(4F16.1))
 3000	FORMAT(///'ERROR CODE FROM FFT = ',I6,///)
	END

---