C	PROGRAM SPHERE
C	DISPLAY ROTATING SPHERE OF DOTS
	COMMON C,S
	DIMENSION C(96),S(96)
C
	PIO4=ATAN(1)
	PIR=PIO4/12.
	ZC=40.
	R=10.
	ZCR=ZC/R
	ZS=20.
	S=1.
C*** ZC IS DISTANCE FROM EYE TO CENTER OF SPHERE
C*** R IS RADIUS OF THE SPHERE
C*** ZS IS DISTANCE TO SCREEN
C*** S IS LAST SCALE FACTOR
C*** SET UP TABLES OF SIN AND COS
	DO 2 I=1,96
	RI=FLOAT(I)*PIR
	C(I)=COS(RI)
	S(I)=SIN(RI)
2	CONTINUE
C
	DO 25 IP=1,NP
	IA=IRAN(96)
	IB=IRAN(96)
	ID=IRAN(64)
	RD=S*R*FLOAT(ID)/64.
C
	CIA=C(IA)
	CIAZS=CIA*ZS
	Y1=S(IA)*ZS
	DO 24 IF=1,48
	MF=LAND(IF/8+2,7)*8
	ISA=IREM(0)*256
	IAT=ISA+I
	IBR=IF+IF
C*** STORE POINT(IA,IB+IBR,ID) AT [MF,IAT]
	IB=IREM((IB+IBR)/96)+1
	X1=CIAZS*S(IB)
C	Y1=SIAZS
	Z1=CIA*C(IB)
	ZM=RD/(ZCR-Z1)
C
	IX=X1*ZM
	IY=Y1*ZM