COMMON IGX1, IGX2
	DIMENSION IGX1(1600), IGX2(1600)
	DIMENSION IXS(480), IYS(480)
	CALL GCLEAR
	CALL GSTART
	DO 80 N = 1,601,5
	R = FLOAT(N-1) / 100.0
	IX = IFIX(R+100.5) - 100
	X = R - FLOAT(IX)
	IF (IREM(IX/2)) 10,20,10
10	X = -X
20	JX1 = IFIX(180.0*X + 90.5)
	JY1 = 99
	DO 70 I = 1,N
	THETA = 2.0*3.141593*FLOAT(I) / FLOAT(N)
	X = R * COS(THETA)
	Y = R * SIN(THETA)
	IX = IFIX(X+100.5) - 100
	IY = IFIX(Y+100.5) - 100
	X = X - FLOAT(IX)
	Y = Y - FLOAT(IY)
	IF (IREM(IX/2)) 30,40,30
30	X = -X
40	IF (IREM(IY/2)) 50,60,50
50	Y = -Y
60	JX2 = IFIX(180.0*X + 90.5)
	JY2 = IFIX(198.0*Y + 99.5)
	IXS(I) = JX2
	IYS(I) = JY2
70	CONTINUE
	CALL GCLEAR
	CALL GPLOT(0,0,0)
	CALL GPLOT(180,0,1)
	CALL GPLOT(180,198,1)
	CALL GPLOT(0,198,1)
	CALL GPLOT(0,0,1)
	CALL GPLOT(JX1,JY1,0)
	DO 75 I = 1,N
	CALL GPLOT(IXS(I),IYS(I),1)
75	CONTINUE
80	CONTINUE
	CALL GSTOP
	STOP
	END
OAT(N-1) / 100.0
	IX = IFIX(R+100.5) - 100
	X = R - FLOAT(IX)
	IF (IREM(IX/2)) 10,20,10
10	X = -X
20	JX1 = IFIX(180.0*X + 90.5)
	JY1 = 99
	DO 70 I = 1,N
	THETA = 2.0*3.141593*FLOAT(I) / FLOAT(N)
	X = R * COS(THETA)
	Y = R * SIN(THETA