SUBROUTINE WPLOT(I,J,K) WRITE(4,1) I, J, K 1 FORMAT(3A2) RETURN END X(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 CALL WPLOT(0,0,0) CALL WPLOT(180,0,1) CALL WPLOT(180,198,1) CALL WPLOT(0,198,1) CALL WPLOT(0,0,1) CALL WPLOT(JX1,JY1,0) DO 70 I = 1,N THETA = 2.0*3.141593*FLOAT(I) / FLOA