C LAMP.FOR - TEST OF VT11-F COMMON/DFILE/IBUF(4000) DIMENSION X(10),Y(10),Z(10),I3(7),IDBUF(1000) C IT = # TICKS BETWEEN PICTURES IT=180 C PARAMETERS FOR RANDU IX=IFIX(SECNDS(0.)/1000.) IY=IFIX(SECNDS(0.)/1000.) C SCREEN PARAMETERS GT40: G4=768, G5=250 ; GT44: G4=500, G5=0 G4=1000. G5=300. C N = NUMBER OF PTS. (# LINE SEGMENTS + 1) N=5 C *************** CALL INIT(4000) C ****************** I3(1)=1 I3(2)=5 I3(3)=4 I3(4)=3 I3(5)=7 I3(6)=2 I3(7)=6 C *************** P1=3.1415926/4. S9=SIN(P1/10.) C9=COS(P1/10.) F=3.2E-04 C HERE IS THE BIG LOOP 160 S1=0. DO 230 I=1,N CALL FUDGE(W,IX,IY) X(I)=W CALL FUDGE(W1,IX,IY) Y(I)=W1 CALL FUDGE(W2,IX,IY) Z(I)=W2 S=W**2 +W1**2 +W2**2 IF (S.GE.S1) S1=S 230 CONTINUE R=SQRT(250000./S1) DO 270 I=1,N X(I)=X(I)*R Y(I)=Y(I)*R 270 Z(I)=Z(I)*R C C ROTATE ABOUT AXIS C DO 650 K=1,7 II=I3(K) I8=II/4 I7=II/2 - I8*2 I6=II - I8*4 - I7*2 IZ=I8 + I7 + I6 Z1 = SQRT (FLOAT ( IZ)) S= S9/ Z1 IF (IZ . GE. 3) S=S+F C= 1. - (1. - C9) / FLOAT (IZ) C C PLOT IT NOW C CALL INIT CALL WINDW (-G4, -G4, G4, G4) DO 610 L=1,80 DO 548 M=1,N IF (I8.EQ.0) GO TO 532 V=Y(M) Y(M)= C*V-S*X(M) X(M) = S*V + C*X(M) 532 IF (I7.EQ.0) GO TO 540 V=Z(M) Z(M) = C*V - S*X(M) X(M) = S*V + C*X(M) 540 IF (I6.EQ.0) GO TO 548 V=Y(M) Y(M) = C*V - S*Z(M) Z(M) = S*V + C*Z(M) 548 CONTINUE Z1=X(1) Z2=Y(1)-G5 CALL APNT(Z1,Z2,0,-5) DO 600 I=2,N Z1=X(I)-X(I-1) Z2=Y(I)-Y(I-1) 600 CALL VECT(Z1,Z2) 610 CONTINUE C CALL TIME(IT) C630 CALL TIMR(IE) C IF (IE.NE.0) GO TO 630 630 CALL WAIT(300.,0.) 650 CONTINUE CALL CMPRS GO TO 160 END SUBROUTINE FUDGE(X,IX,IY) CALL RANDU(IX,IY,Z) IX=IY X=600.*Z-200. IF(X.LE.70) X=X-140. RETURN END