C-PS/8 FOCAL, 1971

01.01 C     FACTOR BY STEVE POULSEN & JOHN HUISMAN
01.02 C     OMSI SOFTWARE DEVELOPMENT GROUP
01.03 C     LAST CHANGE: 9/11/71
01.04 
01.05 C THE PRIME FACTORS OF POSITIVE INTEGERS WILL BE
01.06 C OUTPUT.  INPUT ZERO, THEN ANOTHER NUMBER AND
01.07 C THE FIRST PRIME NOT LESS THAN THAT NUMBER WILL
01.08 C BE PRINTED.
01.09 
01.10 A !! "NUMBER TO BE FACTORED? ",N;I (N)1.1,2.1,1.2
01.20 S C=0;S L=0;S S=FSQT(N);S X=2; G 1.4
01.30 I (2-X)1.35;S X=1
01.35 S X=X+2
01.36 I (N-X)1.8
01.40 I (N/X-FITR(N/X))1.6,1.5,1.6
01.50 S N=N/X;S C=C+1;I (N-1)1.6,1.6;G 1.4
01.60 I (C) 1.7,1.7;T %FITR(FLOG(X)/FLOG(10)+1),X
01.61 I (C-1)1.62,1.62;T " ^",%FITR(FLOG(C)/FLOG(10)+1),C
01.62 I (N-1)1.1,1.1;T " *"
01.65 S C=0;S L=L+1;I (L-3)1.7
01.70 I (X-S)1.3
01.80 T %FITR(FLOG(N)/FLOG(10)+1),N;G

02.10 A !"NUMBER? ",N
02.20 S C=FITR(FSQT(N));S A=1
02.30 S A=A+1;I (A-C)2.4,2.4,2.5
02.40 I ((FITR(N/A)-(N/A)))2.3;S N=N+1;G 2.2
02.50 T %6,!"FIRST PRIME IS",N,!!!
02.60 G 1.1