(Program
STEIN3,
continued from previous page)
1470 REM
-
sigma
=
mean square deviation
-
1480 REM (Wington eqn
8-29)
1490 SIGMA=SQR(SIGMA/(N-3))
1500 REM
'"
print coefficients and estimated errors
-
1510 REM
(eqns 8-26 and 8-30 of Bevington)
1520
F$="W####AA"A
+I-
##.lWfAAAA
)'
1530 PRlNT
1540 PRINT '1TT
=
C;
1550 PRlNT USlNG FS; C(1), SIGMA'SQR(E(1
,I))
1560 PRINT
-
+
C;
1570 PRlNT USlNG
Ff;
C(2). SIGMA'SQR(E(Z.2));
1580 PRINT
'
'
In (R)'
1590 PRINT
'
+
C;
1600 PRlNT USlNG FS; C(3), SIGMA'SQR(E(3,3));
1610 PRINT'
'
In (R)
"
3'
1620 PRlNT
1630 G$="Key in: C1 C2 C3'
1640 P$=-
#.W #.W 1.W'
1650 PRlNT GS
1660 PRlNT USlNG
F5:
C[lj'lOW!, C[Zj'10000!, C[3]'1E+07
1700
'
1710 PRlNT
1712
PRINT
-
T
T
r
1714 PRINT- R ACTUAL CALC ERROR'
1716PRINT'
=====I==
=I=======
==========
==tn====='
1718 P$=
'
#WUW##
####.## ####.##
####.W
1720 FOR L=l TO N
1730 X=LOG(R(L))
1740
TCALC=lI(C(1)+C(2).X+C(3)'XXXXX)-273.15
1760 PRlNT USlNG
P$;R(L).T(L).TCALC,T(L)-TCALC
1780 NEXT
L
REM
"..""""-".
program end
"""""~
1892 END
2OO0 REM
......
*."
Begin Subroutines
"'"""-
2010 REM
""
Invert 3 x 3 matrix "A' by cofactors
*"
2020 GOSUB 2160
:
GOSUB 2120
2030 DET=SUM
2040 FOR K=l TO
3
:
FOR L-1 TO 3
2050 GOSUB 2160
2060 FOR
J-1 TO 3
:
B(J.L)=O
,
B(K.J)=O
:
NEXT
J
:
B(K,L)=l
2070 GOSUB 2120
2080 REM
""
'E'
=
Inverse
=
transpose of colactor
-
2090 E(L,K)=SUM/DET
2100 NEXT L
:
NEXT
K
2110 RETURN
(Continued on next
page)
Artisan Technology Group - Quality Instrumentation ... Guaranteed | (888) 88-SOURCE | www.artisantg.com
To Next Page
Loading...
Need help?
General
