C-FOCAL , 8/68

01.01 C PROGRAM ONE--LEAST SQUARES FIT TO AN EXPOTENTIAL
01.02 E
01.10 S IN=10;A ?N ?,!;F J=1,N;D 3.1
01.20 A ?NOP ?,!;S NO=1/(10^NOP)
01.3p S T=0;D 10;F J=1,N;D 3.2
01.40 I (T1-T)6.2,6.7,6.7

03.10 T %2.0,"Y",J,":";A Y(J);T "X",J,":",;A X(J),!
03.20 D 5;S T=T+((Y(J)-Y1)^2)/Y(J)

05.10 S Y1=FEXP(AL*X(J))
05.20 S Y1=A*Y1

06.20 S IN=-IN/10
06.70 I (FABS(IN)-NO)9.1,9.1,6.8
06.80 S AL=AL+IN;S T1=T;G 1.3

09.10 T %,"ALPHA =",AL,!;D 10;T "A =",A,!;S T=0;F J=1,1,N;D 3.2
09.20 T "F TEST: ",T/(N-2),!!
00.30 QUIT

10.05 S A=0;F J=1,1,N;D 11.1
10.07 S A2=0;S AL=AL*2;F J=1,N;D 11.2
10.10 S AL=AL/2;S A=A/A2

11.10 D 5.1;S A=A+Y1*Y(J)
11.20 D 5.1;S A2=A2+Y1
*@@@@@@@@@O
C-FOCAL , 8/68

01.01 C PROGRAM TWO--LEAST SQUARES FIT TO AN EXPOTENTIAL
01.02 E
01.10 S IN=10;A ?N ?,!;F J=1,N;D 3.1
01.20 A ?NOP ?,!;S NO=1/(10^NOP);A ?XSTART ?,!?DELTAX ?,!
01.30 S T=0;D 10;F J=1,N;D 3.2
01.40 I (T1-T)6.2,6.7,6.7

03.10 T %2.0,"Y",J,":";A Y(J),!
03.20 D 5;S T=T+((Y(J)-Y1)^2)/Y(J)

05.10 S Y1=FEXP(AL*(XS+J*DE))
05.20 S Y1=A*Y1

06.20 S IN=-IN/10
06.70 I (FABS(IN)-NO)9.1,9.1,6.8
06.80 S AL=AL+IN;S T1=T;G 1.3

09.10 T %,"ALPHA =",AL,!;D 10;T "A=",A,!;S T=0;F J=1,1,N;D 3.2
09.20 T "F TEST: ",T/(N-2),!!
09.30 QUIT

10.05 S A=0;F J=1,1,N;D 11.1
10.07 S A2=0;S AL=AL*2;F J=1,N;D 11.2
10.10 S AL=AL/2;S A=A/A2

11.10 D 5.1;S A=A+Y1*Y(J)
11.20 D 5.1;S A2=A2+Y1
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-FOCAL , 8/68

01.01 C PROGRAM THREE--LEAST SQUARES FIT TO AN EXPOTENTIAL
01.02 E
01.10 S IN=10;A ?N ?,!;F J=1,N;D 2
01.r0 A ?NOP ?,!;S NO=1/(10^NOP)
01.30 S T=0;D 10; F J=1,N;D 3.2
01.40 I (T1-T)6.2,6.7,6.7

02.10 T %2.0,"Y",J,":";A M;T "X",J,":";A L,!
02.20 S L=FNEW(J,L);S M=FNEW(N+J,M)

03.20 D 5;S T=T+((FNEW(N+J)-Y1)^2)/FNEW(N+J)

05.10 S Y1=FEXP(AL*FNEW(J))
05.20 S Y1=A*Y1

06.20 S IN=-IN/10
06.70 I (FABS(IN)-NO)9.1,9.1,6.8
06.80 S AL=AL+IN;S T1=T;G 1.3

09.10 T %,"ALPHA =",AL,!;D 10;T "A =",A,!;S T=0;F J=1,1,N;D 3.2
09.20 T "F TEST: ",T/(N-2),!!
09.30 QUIT

10.05 S A=0;F J=1,1,N;D 11.1
10.07 S A2=0;S AL=AL*2;F J=1,N;D 11.2
10.10 S AL=AL/2;S A=A/A2

11.10 D 5.1;S A=A+Y1*FNEW(N+J)
11.20 D 5.1;S A2=A2+Y1
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
C-FOCAL , 8/68

01.01 C PROGRAM FOUR--LEAST SQUARES FIT TO AN EXPOTENTIAL
01.02 E
01.10 S IN=10;A ?N ?,!;F J=1,N;D 3.1
01.20 A ?NOP ?,!;S NO=1/(10^NOP);A ?XSTART ?,!?DELTAX ?,!
01.30 S T=0;D 10;F J=1,N;D 3.2
01.40 I (T1-T)6.2,6.7,6.7

03.10 T %2.0,"Y",J,":";A M,!;S M=FNEW(J,M)
03.20 D 5;S T=T+((FNEW(J)-Y1)^2)/FNEW(J)

05.10 S Y1=FEXP(AL*(XS+J*DE))
05.20 S Y1=A*Y1

06.20 S IN=-IN/10
06.70 I (FABS(IN)-NO)9.1,9.1,6.8
06.80 S AL=AL+IN;S T1=T;G 1.3

09.10 T %,"ALPHA =",AL,!;D 10;T "A =",A,!;S T=0;F J=1,1,N;D 3.2
09.20 T "F TEST: ",T/(N-2),!!
09.30 QUIT

10.05 S A=0;F J=1,1,N;D 11.1
10.07 S A2=0;S AL=AL*2;F J=1,N;D 11.2
10.10 S AL=AL/2;S A=A/A2

11.10 D 5.1;S A=A+Y1*FNEW(N+J)
11.20 D 5.1;S A2=A2+Y1
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@