@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-8K TV5A  @1969

01.02 T !!!!!!!!!!,"        LSQ STERN-VOLMER:LEAST SQUARES TREATMENT"
01.04 T " OF THE GENERAL",!,"                         STERN-VOLMER"
01.06 T " EQUATION",!,"               (UNIVERSITY OF TOLEDO,CHEMISTRY,"
01.07 T "J.E.GANO)"
01.08 T !!,"                R0/RN=(1+C*QN)(1+A*QN)/(1+A*B*QN)"
01.09 T !,"              A=KQT(T)   B=0(S)/(0(S)+0(T))   C=KQT(S)
01.10 T !!,"INSTRUCTIONS:",!,"TURN ON X,Y RECORDER BEFORE STARTING"
01.11 T " IF USING PLOTTING OPTION"
01.12 T !,"ORDINATE GOES FROM 0.8 TO MAXIMUM VALUE OF R0/RN PLUS 10%"
01.14 T !,"ABSCISSA GOES FROM ZERO .......................QN PLUS 10%"
01.15 T !,"ANSWER QUESTIONS WITH YES OR NO"
01.20 S N=1;S XR=0;S YT=0

02.01 T !,"TYPE IN ESTIMATED VALUES OF KQT(T),0(S)/(0(S)+0(T)),KQT(S) AND"
02.03 T " THEN",!,"R0,RN,QN,ETC AND TERMINATE WITH -1,-1.",!
02.05 A A0,B0,C0,R0
02.07 A RN,QN
02.09 I (RN)2.19;I (XR-QN)2.27
02.11 S P(N)=RN
02.13 S R(N)=R0/RN;S Q(N)=QN
02.14 I (YT-R(N))2.29
02.15 S N=N+1
02.17 G 2.07
02.19 S END=N-1
02.21 S MX=XR;S MY=YT
02.23 G 3.12
02.27 S XR=QN*1.1;G 2.11
02.29 S YT=R(N)*1.1;G 2.15

03.02 T !!!!!!!!," RN       R0/RN    QN",!
03.04 S N=1
03.06 T %6.4,P(N)," ",R(N)," ",Q(N),!
03.08 S N=N+1
03.10 I (END-N)4.60,3.06,3.06
03.12 T !!,"RESIDUALS A0       B0       C0"
03.14 S HOLDNO=0;S ITNO=12
03.15 S N=1; S ASUM=0;S BSUM=0;S CSUM=0;S DSUM=0;S ESUM=0;S GSUM=0
03.16 S YASUM=0;S YBSUM=0;S YCSUM=0;S SQRSUM=0;S AB=-10;S BC=-20;S AC=-30

04.01 S DNOM=(1+A0*B0*Q(N));S Y0=(1+A0*Q(N))*(1+C0*Q(N))/DNOM
04.03 S ADY=Q(N)*(1-B0)*(1+C0*Q(N))/DNOM^2
04.05 S BDY=-A0*Q(N)*(1+A0*Q(N))*(1+C0*Q(N))/DNOM^2
04.06 S CDY=Q(N)*(1+A0*Q(N))/DNOM
04.07 S ASUM=ASUM+ADY^2;S BSUM=BSUM+BDY^2;S CSUM=CSUM+CDY^2
04.09 S DSUM=DSUM+ADY*BDY;S ESUM=ESUM+BDY*CDY;S GSUM=GSUM+ADY*CDY
04.13 S YASUM=YASUM+ADY*(R(N)-Y0);S YBSUM=YBSUM+BDY*(R(N)-Y0)
04.15 S YCSUM=YCSUM+CDY*(R(N)-Y0)
04.17 S SQRSUM=SQRSUM+(R(N)-Y0)^2
04.19 S N=N+1
04.21 I (END-N)4.23,4.01,4.01
04.23 S UBC=BSUM*CSUM-ESUM^2;S VAC=ASUM*CSUM-GSUM^2
04.25 S WAB=ASUM*BSUM-DSUM^2;S XFE=GSUM*ESUM-DSUM*CSUM
04.27 S YDE=DSUM*ESUM-BSUM*GSUM;S ZDF=DSUM*GSUM-ASUM*ESUM
04.30 S DETVAL=ASUM*UBC+2*DSUM*ESUM*GSUM-CSUM*(DSUM^2)-BSUM*(GSUM^2)
04.32 S XXX=SQRSUM/(END-3)*DETVAL;S ARANGE=FSQT(XXX*UBC)
04.34 S BRANGE=FSQT(XXX*VAC);S CRANGE=FSQT(XXX*WAB)
04.35 I (ARANGE*BRANGE)4.36,4.80
04.36 S P(AB)=XXX*XFE/ARANGE*BRANGE;S P(BC)=XXX*ZDF/BRANGE*CRANGE
04.37 S P(AC)=XXX*YDE/ARANGE*CRANGE
04.38 S ERROR=SQRSUM-HOLDNO
04.40 T %6.4,!,SQRSUM," ",A0," ",B0," ",C0
04.41 S ITNO=ITNO-1;I (ITNO-1)4.55
04.42 I (ERROR)4.44,4.55,4.44
04.44 S A0=A0+(YASUM*UBC+YBSUM*XFE+YCSUM*YDE)/DETVAL
04.46 S B0=B0+(YASUM*XFE+YBSUM*VAC+YCSUM*ZDF)/DETVAL
04.47 S C0=C0+(YASUM*YDE+YBSUM*ZDF+YCSUM*WAB)/DETVAL
04.48 S HOLDNO=SQRSUM
04.50 G 3.15
04.55 G 3.02
04.60 T !,"KQT(T)=",A0," (SIGMA=",ARANGE,")",!!
04.62 T "0S/(0S+0T)=",B0," (SIGMA=",BRANGE,")"
04.63 T !!,"KQT(S)=",C0," (SIGMA=",CRANGE,")"
04.64 T !!,"CORRELATION COEFFICIENTS",!,"  P(AB)=",P(AB),!,"  P(BC)="
04.66 T P(BC),!,"  P(AC)=",P(AC)
04.70 A !!!!!!!!!,"PLOT EXPERIMENTAL?",REPLY
04.72 I (REPLY-155)4.74,10.26;I (0.25E+21-REPLY)4.74,5.10,4.74
04.74 T !,"YOU GOOFED. TRY AGAIN!";G 4.70
04.80 T !!,"DIVISION BY ZERO REQUESTED. IGNORE P(  ) VALUES";G 4.38

05.10 DO 31.0
05.11 A "MANUAL SCALING?",RE
05.12 I (RE-155)5.13,5.15;I (.25E+21-RE)5.13,5.14
05.13 T "WRONG ANSWER",!;G 5.11
05.14 DO 31.30;G 5.18
05.15 SET XL =0;SET YB=0;S XR=MX;S YT=MY
05.18 DO 31.4
05.20 SET N=1;SET IR=R(END)/500;SET IQ=Q(END)/500
05.22 SET X=Q(N)-IQ*4;SET Y=R(N)-IR*6
05.25 DO 31.50
05.26 DO 31.90

06.10 SET X=Q(N)+IQ*4
06.13 DO 31.5
06.15 SET Y=R(N)+IR*6
06.17 DO 31.5
06.20 SET X=Q(N)-IQ*4
06.21 DO 31.5
06.25 SET Y=R(N)-IR*6
06.27 DO 31.5
06.30 DO 31.80
06.32 S N=N+1
06.34 I (END-N)8.1,5.22,5.22

08.10 T !,"PLOT THEORETICAL?"
08.20 A PLOT
08.40 I (PL0T-155)8.42,10.26,8.42
08.42 I (PLOT-(.25E+21))8.45,9.1,8.45
08.45 T "CANNOT YOU SPEAK ENGLISH? TRY AGAIN",!
08.47 G 8.20

09.10 S X=0;S Y=(1+A0*X)*(1+C0*X)/(1+A0*B0*X)
09.40 DO 31.5; DO 31.9
09.60 SET X=X+IQ;SET Y=(1+A0*X)*(1+C0*X)/(1+A0*B0*X)
09.80 I (X-Q(END))10.1,10.2,10.2

10.10 DO 31.5
10.15 GOTO 9.6
10.20 D0 31.8
10.26 T !,"TERMINATE?"
10.27 A TERM
10.28 I (TERM-155)10.29,10.3,10.29
10.29 I (TERM-(.25E+21))10.45,10.4,10.45
10.30 T !!!!!!!!!!
10.32 G 1.20
10.40 Q
10.45 T "YOU CANNOT SPEAK ENGLISH, CAN YOU? TRY AGAIN.",!
10.47 G 10.27

31.01 C FOCAL PLOT ROUTINES
31.10 S SI=1
31.11 S X=11;S Y=-245;D 31.16;I (SI-16)31.12,31.17,31.12
31.12 S X=511;D 31.16;I (SI)31.17
31.13 S Y=255;D 31.16;I (SI)31.17
31.14 S X=11;D 31.16;I (SI)31.17
31.15 G 31.11
31.16 S H(2040)=FDIS(X,Y);A SI
31.17 R ;Q
31.30 D 31.34;D 31.35;D 31.40
31.31 Q
31.34 A "MIN. X VALUE "XL,!,"MAX. X VALUE "XR,!
31.35 A "MIN. Y VALUE "YB,!,"MAX. Y VALUE "YT,!
31.40 S CX(2040)=500./(XR-XL);S CY(2040)=500./(YT-YB)
31.41 Q
31.50 D 31.52;D 31.53;D 31.54;S H(2040)=FDIS(XP(2040),YP(2040))
31.51 Q
31.52 S XP(2040)=11.+(X-XL)*CX(2040)
31.53 S YP(2040)=-245.+(Y-YB)*CY(2040)
31.54 I (11.-XP(2040))31.55;S XP(2040)=11.;G 31.56
31.55 I (XP(2040)-511.)31.56;S XP(2040)=511.;G 31.56
31.56 I (-245.-YP(2040))31.57;S YP(2040)=-245.;R
31.57 I (YP(2040)-255.)31.58;S YP(2040)=255.;R
31.58 R
31.80 D 31.84;S H(2040)=FPDN(0);D 31.86
31.81 Q
31.84 F ID(2040)=1,270*.5;C
31.86 F ID(2040)=1,270*.5;C
31.90 D 31.94;S H(2040)=FPDN(1);D 31.96
31.91 Q
31.94 F ID(2040)=1,270*.5;C
31.96 F ID(2040)=1,270*.5;C
*@@@@@@@@@@@@@@