C FOCAL15 V6A
01.10 A "RADIUS= ",R,!,"ORIGIN   X= ",X,"     Y= ",Y,!,"NO. OF"
01.11 A " POINTS= ",N,!!
01.15 S X1=X-R;S I=(4*R)/N; S NN=N/2; S RR=R*R
01.20 F J=1,NN; S X(J)=X1; S Y(J)=FSQT(FABS(RR-(X1-X)^2))+Y;S X1=X1+I
01.25 S X1=R+X
01.28 T X1,!
01.30 F J=NN+1,N;S X(J)=X1;S Y(J)=-FSQT(FABS(RR-(X1-X)^2))+Y;S X1=X1-I
01.35 D 10
01.40 GOTO 11.1
02.05 S MAX=R-X
02.10 T "TRANSFORM CONSTANT  (MAX = ",MAX,")  ";A LL,!;S L=LL^2
02.12 S ZL=FSIN(Y/(X+LL))/FCOS(Y/(X+LL)); T "ANGLE OF 0 LIFT = ",ZL*180/3.14,!
02.20 F J=1,N;S Z=X(J)^2+Y(J)^2;S X(J)=(L*X(J))/Z+X(J);S Y(J)=(-L*Y(J))/Z+Y(J)
02.30 D 10
05.10 S MAX=-10000; S MIN=100000
05.15 F J=1,N; D 8
05.20 A "SIZE OF PLOT ",SP,!
05.25 S XSC=(SP)/(MAX-MIN)
05.28 S YSC=(SP*(2/5))/(MAX-MIN)
05.29 S MM=-(XSC*MIN)
05.30 F J=1,N; S Y(J)=FITR(Y(J)*XSC+MM)+5;S X(J)=FITR(X(J)*YSC+MM)+5
05.35 S M=-1; F J=1,SP; S M=M+1; D 6
06.05 S C=0;S A(2)=0
06.10 F K=1,N; D 7
06.20 IF (C-1) 6.9,6.8,6.3
06.30 IF (A(1)-A(2)) 6.4,6.5,6.5
06.40 S T=A(1);S A(1)=A(2);S A(2)=T
06.50 F JJ=0,A(2); T " "
06.55 T "*"
06.80 F JJ=0,A(1)-A(2); T " "
06.85 T "*"
06.90 T !
07.10 IF (M-X(K))7.3,7.5,7.3
07.30 RETURN
07.50 S C=C+1; S A(C)=Y(K); RETURN
08.10 IF (MAX-X(J)) 8.15,8.2,8.2
08.15 S MAX=X(J)
08.20 IF (MAX-Y(J)) 8.25,8.3,8.3
08.25 S MAX=Y(J)
08.30 IF (Y(J)-MIN) 8.35,8.4,8.4
08.35 S MIN=Y(J)
08.40 IF (X[J]-MAX) 8.45,8.5,8.5
08.45 S MIN=X[J]
08.50 R
10.05 T "PLOT READY ",!!
10.10 F J=1,N; S XX=FDPY(2048*X(J),2048*Y(J))
10.20 GOTO 10.1
11.10 S Z=0