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

(Source file text) 

---

	DIMENSION AMP(32),PHASE(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 PHAMPL(N,REAL,IMAG,PHASE,AMP)
	TYPE 2000
	TYPE 1001,REAL
	TYPE 1002,IMAG
	TYPE 2001,PHASE
	TYPE 2002,AMP

	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(/' ***PHASE ANGLES***',/,(4F16.5))
 2002	FORMAT(' ***AMPLITUDES***',/,(4F16.5))
 3000	FORMAT(///'ERROR CODE FROM FFT = ',I6,///)
	END

---