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

(Source file text) 

---

	DIMENSION IR(32),IM(32),P(32),A(32)
	DATA N,M/0,0/
	NEXT(V)=(RAN(N,M)*1000.)-500.

	DO 10 I=1,32
	IR(I)=NEXT(X)
   10	IM(I)=NEXT(X)
	TYPE 1002,IR
	TYPE 1003,IM

	DO 11 I=1,32
	X0=IR(I)
	X1=IM(I)
	P(I)=ATAN2(X1,X0)
   11	A(I)=SQRT(X1*X1+X0*X0)
	TYPE 1006,P
	TYPE 1007,A

	CALL PHAMPL(32,IR,IM,P,A)
	TYPE 1004,P
	TYPE 1005,A

 1002	FORMAT(' RANDOM INPUT DATA'/' REAL PART',/,(4I15))
 1003	FORMAT(/,'    IMAGINARY PART',/,(4I15))
 1004	FORMAT(//,'    PHASE ANGLES FROM SUBROUTINE',/,(4F15.5))
 1005	FORMAT(/,' AMPLITUDES FROM SUBROUTINE',/,(4F15.5))
 1006	FORMAT(//,' PHASE ANGLES FROM DIRECT CALCULATIONS',/,(4F15.5))
 1007	FORMAT(/,' AMPLITUDES FROM DIRECT CALCULATIONS',/,(4F15.5))
	STOP
	END

---