10	!
	!	COMPLEX LOG
	!
	!	KEYED IN BY TIMOTHY SHANE MCCOY 
	!	AT CENTRAL STATE UNIVERSITY
	!	ON PDP 11/45
	!	VERSION RSTS/E  V06C
	!
20	INPUT C1,C2
30	S1,S2,A1,A2 = 0
40	X1 = ((C1-1) * (C1+1) + C2 * C2)/((C1+1) * (C1+1) + C2 * C2)
50	X2 = (C2 * (C1+1) - C2 * (C1-1))/((C1+1) * (C1+1) + C2 * C2)
55	PRINT X1;X2
60	FOR N = 1 TO 100
70	P1 = 1
80	P2 = 0
90	FOR J = 1 TO 2*N-1
100	T1 = P1 * X1 - P2 * X2
110	T2 = P1 * X2 + P2 * X1
112	P1 = T1
114	P2 = T2
120	NEXT J
130	S1 = S1 + P1/(2*N-1)
140	S2 = S2 + P2/(2*N-1)
150	IF ABS(S1 - A1) > 1E - 6
	THEN 170
160	IF ABS(S2-A2) < 1E - 6
	THEN 190
170	A1 = S1
172	A2 = S2
175	NEXT N
180	PRINT "DID NOT CONVERGE"
185	GOTO 200
190	PRINT "DID CONVERGE"
200	PRINT "SUM ="2 * S1" +"2 * S2"I"
210	PRINT "LOG =".5 * LOG(C1*C1+C2*C2)" +"ATN(C2/C1)"I"
220	GOTO 10
230	END