IMD 1.17: 21/02/2009 17:00:53 AS-H733A-BC MINC APP PKG V1.0 BIN RX01  @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ߋtv ?BOOT-F-No boot on volume xxx˵)ӵU$ aC}  `  7>wv7& &U@wg @@x 6, N 5z 7  b5@ 7RB  :% eC % C ``` e̋̋ @w,U? <@4e,}`8 p$ 7 55X 7  #̋ 7̋4 (xx~    MINC APP PKGAS-H733A-BC DECRT11A ve< m  '5xg  'd  '芻  '|QT 'df '+y 'df  '; '&<f& { ' { 'Y{  ')&v  'af{  'T  'zTq  GOq<2Oq  G<p  OqԀ- &Ԁz$H   5x &. " &.$TW &.f &.Z F.y{ f._Xf F,Q<g0*2  . e!  . R  . ( N . Z  . Z<.  'f ',Q 'T '  'T 'lx 'x, '(: 'Oq  y}<. R<. y . '   'y::  'Y$;  'fW  ' NAd  '_ 'q '~ 'Uqc 'OqAd '(M 't: 'R '   's  'y}  '&  '3:  &.5x &. " &.$TW &.f &.Z F.y{ f._Xf F,Q<g0*2  . e!  . R  . ( N . Z  . Z<.  y}<. R<. y . '  mOqAd '(M 't: 'R '   's  'y}  '&  '3:  &.5x &. " &.$TW &.f &.Z F.y{ f._Xf F,Q<g0*2  .e!  .R  .( N .Z  .Z<.y}<.R<.y .vL@8 /|QE /}<F3TOq GOqȻ  G|Q:  1dfOq  GyqK  Gzy  GOq   9 12 113.3 10 68 8 12 109.4 10 C=0 \ D=0 1220 REM - LOOP FOR EACH SUBINTERVAL 1230 FOR I=1 TO N-1 1240 IF S=0 THEN A(0,I)=FNC(A+I*H) 1250 IF I/2=INT(I/21000 A=SYS(7,1) \ REM ROOTS OF QUADRATIC EQUATIONS 1010 REM - ENTER COEFFICIENTS A,B,C OF A*X^2+B*X+C 1020 PRINT "COEFFICIENTS) THEN D=D+A(0,I) \ GO TO 1270 \ REM EVEN 1260 C=C+A(0,I) \ REM ODD 1270 NEXT I 1280 PRINT "INTEGRAL =";H/3*(A(0,0)+4*C+2*D+A A,B,C "; \ INPUT A,B,C 1030 D=B^2-4*A*C 1040 E=SQR(ABS(D)) 1050 IF D<0 THEN GO TO 1080 1060 PRINT "ROOTS (REAL): ";(-B-E)/(0,N)) 1290 RETURN (2*A);", ";(-B+E)/(2*A) 1070 GO TO 1090 1080 PRINT "ROOTS (COMPLEX): ";-B/(2*A);"+ OR -";E/(2*A);" I" 1090 PRINT 1100 RETURN 1000 A=SYS(7,1) \ REM INTEGRATION: TRAPEZOIDAL RULE 1010 PRINT 1020 PRINT "FORMULA MUST BE ENTERED AT LINE 20 IN NUM.BAS, PRO1000 A=SYS(7,1) \ REM INTEGRATION: SIMPSON'S RULE 1010 PRINT 1020 PRINT "IS FORMULA KNOWN (AND ENTERED AT LINE 20)"; \ INPUTCEED "; \ INPUT Y$ 1030 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1210 1040 PRINT "INTEGRATION LIMITS (LOWER, UPPER)"; 1050 INPUT A,B  Y$ 1030 PRINT "INTEGRATION LIMITS (LOWER, UPPER)"; \ INPUT A,B 1040 PRINT "NUMBER OF INTERVALS"; \ INPUT N 1050 IF N>A8 THEN PRINT "TOO MANY INTERVALS" \ GO TO 1290 1060 IF N/2<>INT(N/2) THEN PRINT "NUMBER MUST BE EVEN" \ GO TO 1040 1070 H=(B-A)/N 1080 IF SEG$(Y$,1,1)="Y" THEN S=0 \ A(0,0)=FNC(A) \ A(0,N)=FNC(B) \ GO TO 1210 1090 FOR I=0 TO N 1100 PRINT "VALUE OF F(X) AT INTERVAL ";I;" (X=";A+I*H;")"; \ INPUT A(0,I) 1110 NEXT I \ S=1 1120 PRINT "DO YOU WISH TO CHECK/CHANGE THE DATA"; \ INPUT Y$  13 , 5 7 26 6 60 78.5 1 29 15 52 74.3 11 56 8 20 104.3 11 31 8 41130 IF SEG$(Y$,1,1)<>"Y" GO TO 1210 1140 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1150 PRINT "NUMBER OF INTE7 87.6 7 52 6 33 95.9 11 55 9 22 109.2 3 71 17 6 102.7 1 31 22 RVALS = ";N; \ INPUT Y$ 1160 IF SEG$(Y$,1,1)="Y" THEN PRINT "NUMBER OF INTERVALS = "; \ INPUT N\H=(B-A)/N 1170 FOR I=0 TO N 1 44 72.5 2 54 18 22 93.1 21 47 4 26 115.9 1 40 23 34 83.8 11 66 180 PRINT "F(X";I;") = ";A(0,I); \ INPUT Y$ 1190 IF SEG$(Y$,1,1)="Y" THEN PRINT "F(X";I;") = "; \ INPUT A(0,I) 1200 NEXT I 12  TO L 1340 PRINT "INTERPOLATE X =";A(2,I); \ INPUT Y$ 1350 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW X = "; \ INPUT A(2,I) 1360 NEXT I 1370 FOR K=1 TO L 1380 Y=0 1390 FOR J=1 TO P 1400 X=1 1410 FOR I=1 TO P 1420 IF I=J GO TO 1440 1430 X=X*(A(2,K)-A(0,I1000 A=SYS(7,1) \ REM CURVILINEAR INTERPOLATION 1010 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN))/(A(0,J)-A(0,I)) 1440 NEXT I 1450 Y=Y+X*A(1,J) 1460 NEXT J 1470 A(3,K)=Y 1480 NEXT K 1490 PRINT "X-VALUE";TAB(10);"Y-VAL I9=0 \ GO TO 1050 1030 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "NUMBER OUE" 1500 PRINT 1510 FOR I=1 TO L 1520 PRINT A(2,I);TAB(10);A(3,I) 1530 NEXT I 1540 PRINT 1550 RETURN F KNOWN POINTS? "; \ INPUT #I9,P \ IF I9<>0 THEN PRINT P 1060 IF P<=A8 GO TO 1090 1070 PRINT "TOO MANY POINTS" \ IF I9<>0 THEN CLOSE #I9 1080 GO TO 1550 1090 FOR I=1 TO P 1100 REM - ENTER COORDINATES OF KNOWN POINTS ON CURVE 1110 PRINT "X,Y OF POINT"1000 A=SYS(7,1) \ REM MATRIX MULTIPLICATION 1010 FOR L=1 TO 2 1020 PRINT "IS MATRIX";L;"DATA IN A FILE"; \ INPUT Y$ 1030 IF S;I;"? "; 1120 INPUT #I9,A(0,I),A(1,I) \ IF I9<>0 THEN PRINT A(0,I);A(1,I) 1130 NEXT I 1140 PRINT "ENTER X COORDINATES TO BE IEG$(Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1060 1040 PRINT "FILE SPECIFICATION"; \ INPUT Y$ 1050 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 NTERPOLATED (99999 TO FINISH)" 1150 FOR I=1 TO A8 1160 PRINT "INTERPOLATION POINT";I;" X = "; 1170 INPUT #I9,A(2,I) \ IF I9<1060 PRINT "MATRIX";L;"BOUNDS? "; \ INPUT #I9,C,D 1070 IF I9<>0 THEN PRINT C,D 1080 IF L=1 GO TO 1130 1090 IF B=C GO TO 1130 >0 THEN PRINT A(2,I) 1180 IF A(2,I)=99999 THEN L=I-1 \ GO TO 1220 1190 NEXT I 1200 L=A8 1210 PRINT "ARRAY FULL - WILL INTERP 1100 PRINT "NUMBER OF ROWS MUST EQUAL NUMBER OF COLUMNS OF FIRST MATRIX" 1110 IF I9=0 THEN GO TO 1060 1120 CLOSE #I9 \ GO TO 1060 PRINT "NUMBER OF INTERVALS"; 1070 INPUT N 1080 REM - D IS THE SIZE OF EACH INTERVAL 1090 D=(B-A)/N 1100 REM - ADD UP TOLATE THE POINTS ENTERED" 1220 IF I9<>0 THEN CLOSE #I9 \ GO TO 1370 1230 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUTHE AREA OF EACH TRAPEZOID 1110 J=A+D 1120 I=FNC(A)*.5 1130 FOR K=1 TO N-1 1140 I=I+FNC(J) 1150 J=J+D 1160 NEXT K 1170 REM Y$ 1240 IF SEG$(Y$,1,1)<>"Y" GO TO 1370 1250 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1260 PRINT "NUMBER OF - COMPUTE INTEGRAL, PRINT 1180 I=I+FNC(J)*.5 1190 PRINT "INTEGRAL =";I*D 1200 PRINT 1210 RETURN  KNOWN POINTS";P; \ INPUT Y$ 1270 IF SEG$(Y$,1,1)="Y" GO TO 1050 1280 FOR I=1 TO P 1290 PRINT "POINT";I;"= ";A(0,I),A(1,I); \ INPUT Y$ 1300 IF SEG$(Y$,1,1)<>"Y" GO TO 1320 1310 PRINT "NEW POINT";I;"= "; \ INPUT A(0,I),A(1,I) 1320 NEXT I 1330 FOR I=1 T #I9,A(I,J) \ IF I9<>0 THEN PRINT A(I,J) 1190 GO TO 1210 1200 INPUT #I9,B(I,J) \ IF I9<>0 THEN PRINT B(I,J) 1210 NEXT J 1220 NEXT I 1230 IF I9<>0 THEN CLOSE #I9 1240 IF L=1 THEN A=C \ B=D 1250 NEXT L 1260 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1270 IF SEG$(Y$,1,1)<>"Y" GO TO 1490 1280 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1290 R=A \1000 A=SYS(7,1) \ REM TRIGONOMETRIC POLYNOMIAL 1010 PRINT 1020 PRINT "NUMBER OF PAIRS OF TERMS IN POLYNOMIAL "; \ INPUT N 10 K=B 1300 FOR L=1 TO 2 1310 PRINT 1320 PRINT "MATRIX";L 1330 PRINT \ PRINT "BOUNDS";R;",";K; \ INPUT Y$ 1340 IF SEG$(Y$,130 FOR I=1 TO N 1040 PRINT "TERMS";I;"= "; 1050 INPUT A(0,I),A(1,I) 1060 NEXT I 1070 PRINT "ANGLE"; \ INPUT X 1080 Y=0 109,1)="Y" GO TO 1060 1350 FOR I=1 TO R 1360 PRINT "ROW";I 1370 FOR J=1 TO K 1380 PRINT "COLUMN";J;"="; 1390 IF L=1 THEN PRINT0 FOR I=1 TO N 1100 Y=Y+A(0,I)*SIN(I*X)+A(1,I)*COS(I*X) 1110 NEXT I 1120 PRINT "F(";X;")=";Y 1130 PRINT "ANOTHER ANGLE"; \ I A(I,J); \ GO TO 1410 1400 PRINT B(I,J); 1410 INPUT Y$ \ IF SEG$(Y$,1,1)<>"Y" GO TO 1450 1420 PRINT "NEW VALUE"; 1430 IF L=1NPUT Y$ 1140 IF SEG$(Y$,1,1)="Y" GO TO 1070 1150 RETURN  THEN INPUT A(I,J) \ GO TO 1450 1440 INPUT B(I,J) 1450 NEXT J 1460 NEXT I 1470 R=C \ K=D 1480 NEXT L 1490 FOR I=1 TO A 151000 A=SYS(7,1) \ REM INTEGRATION: GAUSIAN QUADRATURE 1010 PRINT 1020 PRINT "FORMULA MUST BE ENTERED AT LINE 20 IN NUM.BAS, P00 FOR J=1 TO D 1510 S=0 1520 FOR K=1 TO B 1530 S=S+A(I,K)*B(K,J) 1540 NEXT K 1550 A(I-1,J)=S 1560 NEXT J 1570 NEXT I 15ROCEED "; \ INPUT Y$ 1030 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1250 1040 REM - ABSCISSAS AND WEIGHT FACTORS FOR 20-POINT GAUSSIAN 80 PRINT "WHERE ARE THE RESULTS TO BE PRINTED/SAVED - TT:, LP:, OR FILESPEC"; 1590 INPUT Y$ \ L=0 \ IF SEG$(Y$,1,3)="TT:" THEN INTEGRATION 1050 DATA .076526521,.15275339,.22778585,.14917299,.37370609 1060 DATA .14209611,.510867,.13168864,.63605368,.1181I9=0 \ GO TO 1630 1600 IF LEN(Y$)=0 THEN I9=0 \ GO TO 1630 1610 I9=3 \ OPEN Y$ FOR OUTPUT AS FILE #3 1620 IF SEG$(Y$,1,3)<>"L9453 1070 DATA .74633191,.10193012,.83911697,.083276742,.91223443 1080 DATA .062672048,.96397193,.04060143,.9931286,.017614007P:" THEN L=1 \ PRINT #I9,A;",";D 1630 FOR I=0 TO A-1 1640 FOR J=1 TO D 1650 PRINT #I9,A(I,J); \ IF L=1 THEN PRINT #I9 1660 N 1090 PRINT "INTEGRATION LIMITS (LOWER,UPPER)"; \ INPUT X,Y 1100 PRINT "NUMBER OF INTERVALS"; \ INPUT N 1110 S=(Y-X)/N/2 \ T=EXT J 1670 PRINT #I9 1680 NEXT I 1690 IF I9<>0 THEN CLOSE #I9 1700 PRINT "ANYWHERE ELSE "; \ INPUT Y$ 1710 IF SEG$(Y$,1,1)=X+S \ R=0 1120 REM - COMPUTE INTERGRAL FOR EACH SUBINTERVAL 1130 FOR I=1 TO N 1140 P=0 1150 REM - COMPUTE SUMMATION FACTOR F 1720 1130 FOR I=1 TO C 1140 PRINT "ROW";I 1150 FOR J=1 TO D 1160 PRINT "COLUMN";J;"? "; 1170 IF L<>1 GO TO 1200 1180 INPU"Y" THEN GO TO 1580 1720 RETURN  30 NEXT I 1240 A(0,0)=A(0,N) \ A(1,0)=A(1,N) 1250 X=0 1260 REM - CALCULATE AREA, PRINT 1270 FOR I=1 TO N 1280 X=X+(A(0,I-1)ST RATE ON INVESTMENT,NOMIRT.BAS EFFECTIVE INTEREST RATE ON INVESTMENT,EINTRT.BAS EARNED INTEREST TABLE,EINTTB.BAS DEPRECIATI+A(0,I))*(A(1,I-1)-A(1,I)) 1290 NEXT I 1300 PRINT "AREA =";X/2 1310 PRINT 1320 RETURN ON RATE,DEPRT.BAS DEPRECIATION AMOUNT,DEPRCN.BAS SALVAGE VALUE,SLVGE.BAS PRINCIPAL ON A LOAN,PRINC.BAS REGULAR PAYMENTS ON AOR EACH SUBINTERVAL 1160 FOR J=1 TO 10 1170 READ A,B \ C=S*A+T \ D=T-S*A 1180 P=P+B*(FNC(C)+FNC(D)) 1190 NEXT J 1200 RESET  1210 R=R+P*S \ T=T+2*S 1220 NEXT I 1230 PRINT "INTEGRAL =";R 1240 PRINT 1250 RETURN 1000 A=SYS(7,1) \ REM DERIVATIVE 1010 PRINT "FORMULA MUST BE ENTERED AT LINE 20 IN NUM.BAS, PROCEED "; \ INPUT Y$ 1020 IF SEG$1000 A=SYS(7,1) \ REM "AREA OF A POLYGON" 1010 PRINT 1020 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1030 IF SEG$(Y$,1,1)<>"(Y$,1,1)<>"Y" THEN GO TO 1280 1030 PRINT "ENTER X COORDINATES (99999 TO FINISH)" 1040 FOR I=1 TO A8 1050 PRINT "POINT X(";IY" THEN I9=0 \ GO TO 1060 1040 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1050 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1060 PRINT "N;") = "; \ INPUT A(2,I) 1060 IF A(2,I)=99999 THEN L=I-1 \ GO TO 1100 1070 NEXT I 1080 L=A8 1090 PRINT "VECTOR FULL - WILL COUMBER OF VERTICES ? "; \ INPUT #I9,N \ IF I9<>0 THEN PRINT N 1070 IF N>A8 THEN PRINT "TOO MANY VERTICES" \ GO TO 1320 1080 PRIMPUTE DERIVATIVES FOR X-VALUES ENTERED" 1100 FOR I=0 TO L 1110 D=0 \ B=A(2,I) 1120 REM - CALCULATE DIFFERENCE QUOTIENTS FOR PNT "COORDINATES OF VERTEX 1 ? "; \ INPUT #I9,A(0,1),A(1,1) 1090 IF I9<>0 THEN PRINT A(0,1);",";A(1,1) 1100 FOR I=2 TO N 1110 OINTS APPROACHING X 1130 FOR N=1 TO 10 1140 C=D \ E=.5^N 1150 A=B+E 1160 D=(FNC(A)-FNC(B))/E 1170 IF ABS(D-C)0 THEN PRINT A(0,I);",";A(1,I) 1130 NEXT I \ IF I9<>0 THE-06) GO TO 1200 1180 NEXT N 1190 REM - APPROXIMATE DERIVATIVE OF FUNCTION AT X 1200 A(3,I)=2*D-C 1210 NEXT I 1220 PRINT "XEN CLOSE #I9 1140 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1150 IF SEG$(Y$,1,1)<>"Y" GO TO 1240 1160 PRINT "T-VALUE";TAB(10);"Y-VALUE" 1230 PRINT 1240 FOR I=1 TO L 1250 PRINT A(2,I);TAB(10);A(3,I) 1260 NEXT I 1270 PRINT 1280 RETUYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1170 PRINT "NUMBER OF VERTICES = ";N; \ INPUT Y$ 1180 IF SEG$(Y$,1,1)="Y" THRN EN PRINT "NUMBER OF VERTICES = "; \ INPUT N 1190 PRINT "VERTEX";TAB(10);"X-VALUE";TAB(20);"Y-VALUE" 1200 FOR I=1 TO N 1210 PR19 FUTURE VALUE OF INVESTMENT,FVLINV.BAS FUTURE VALUE OF REGULAR DEPOSITS,FVALRD.BAS REGULAR DEPOSITS,REGDEP.BAS REGULAR WITINT I;TAB(10);A(0,I);TAB(20);A(1,I); \ INPUT Y$ 1220 IF SEG$(Y$,1,1)="Y" THEN PRINT "VERTEX";I;" = "; \ INPUT A(0,I),A(1,I) 12HDRAWALS FROM INVESTMENT,RWFINV.BAS INITIAL INVESTMENT,INIINV.BAS MINIMUM INVESTMENT FOR WITHDRAWALS,MINVW.BAS NOMINAL INTERE  IF A(I)=0 THEN 1190 1140 A=0 1150 FOR J=1 TO A(I) 1160 A=A+LOG(J) 1170 NEXT J 1180 A(I)=A 1190 NEXT I 1200 REM - USING TTHEN I9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "(0-STANDARD, 1HE COMPUTED FACTORIALS, COMPUTE PROBABILITY 1210 R=EXP(A(1)-A(2)-A(3)+X*LOG(P)+(N-X)*LOG(1-P)) 1220 REM BINOMIAL DISTRIBUTION =NON-STANDARD)" 1060 PRINT "WHICH TYPE OF VARIABLE ? "; 1070 INPUT #I9,S \ IF I9<>0 THEN PRINT S 1080 IF S>=0 THEN IF S<=1 GO- PRINT 1230 PRINT 1240 PRINT "PROBABILITY OF";X;"SUCCESSES IN";N;"TRIALS =";R 1250 PRINT 1260 RETURN  TO 1100 1090 IF I9=0 GO TO 1050 \ PRINT "ILLEGAL INPUT" \ CLOSE #I9 \ GO TO 1590 1100 IF S=0 THEN 1170 1110 REM - REQUEST 'NON-STANDARD' VARIABLE DATA 1120 PRINT "MEAN ? "; 1130 INPUT #I9,M \ IF I9<>0 THEN PRINT M 1140 PRINT "STANDARD DEVIATION ? "1000 A=SYS(7,1) \ REM CHI-SQUARE DISTRIBUTION- 1010 PRINT "DEGREES OF FREEDOM"; \ INPUT V 1020 PRINT "CHI-SQUARE"; \ INPUT W ; 1150 INPUT #I9,S \ IF I9<>0 THEN PRINT S 1160 GO TO 1180 1170 S=1 1180 PRINT 1190 PRINT "NUMBER OF ELEMENTS (MAXIMUM= "; LOAN,REGPAY.BAS LAST PAYMENT ON A LOAN,LSTPAY.BAS REMAINING BALANCE ON A LOAN,REMBLC.BAS TERM OF A LOAN,TRMOLN.BAS ANNUAL I1030 REM - R=THE DENOMINATOR PRODUCT 1040 R=1 1050 FOR I=V TO 2 STEP -2 1060 R=R*I 1070 NEXT I 1080 REM - K=THE NUMERATOR PNTEREST RATE ON A LOAN,AINTLN.BAS MORTGAGE AMMORTISATION TABLE,MRTAMO.BAS END RODUCT 1090 K=W^(INT((V+1)/2))*EXP(-W/2)/R 1100 REM - THE PI FACTOR IS USED ONLY WHEN DEG. FREEDOM ARE ODD 1110 IF INT(V/2)<>V/2 THEN J=SQR(2/W/3.14159) \ GO TO 1140 1120 REM - CALCULATE L (SUMMATION FACTOR) 1130 J=1 1140 L=1 1150 M=1 1160 V=V+2  1170 M=M*W/V 1180 REM - CHECK FOR END OF SUMMATION 1190 IF M<1.00000E-07 THEN 1220 1200 L=L+M 1210 GO TO 1160 1220 REM 121000 A=SYS(7,1) \ REM-BINOMIAL DISTRIBUTION - 1010 PRINT "NUMBER OF TRIALS"; \ INPUT N 1020 IF N<=0 THEN PRINT "NUMBER MUST B30 REM CHI-SQUARE DISTRIBUTION-PRINT 1240 PRINT 1250 PRINT "TAIL END VALUE =";1-J*K*L 1260 PRINT 1270 RETURN E >=0" \ GO TO 1010 1030 PRINT "EXACT NUMBER OF SUCESSES"; \ INPUT X 1040 IF X>N THEN PRINT "SUCESSES MUST BE <= TRIALS" \ GO TO 1030 1050 PRINT "PROBABILITY OF SUCCESS"; \ INPUT P 1060 IF P>0 THEN IF P<1 GO TO 1080 1070 PRINT "PROBABILITY MUST LIE BETWEEN 0 AND 1" \ GO TO 1050 1080 REM - COMPUTE THE FACTORIALS 1090 A(1)=N 1100 A(2)=X 1110 A(3)=N-X 1120 FOR I=1 TO 3 11301000 A=SYS(7,1) \ M=0 \ REM - "NORMAL DISTRIBUTION" 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y"  S = "; \ INPUT M,S 1340 PRINT "NUMBER OF ELEMENTS = ";N; \ INPUT Y$ 1350 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUHEN 1250 1220 X=.5/(1+L*(.196854+L*(.115194+L*(3.44000E-04+L*.019527))))^4 1230 X=INT(X*10000+.5)/10000 1240 GO TO 1280 1250T N 1360 IF N>0 THEN IF N<=R8 GO TO 1380 1370 GO TO 1340 1380 FOR I=1 TO N 1390 PRINT "ELEMENT";I;" = ";R(4,I); \ INPUT Y$  L=L*(1+.08*L^4/R^3) 1260 GO TO 1220 1270 REM - ADJUST IF INVERSE WAS COMPUTED 1280 IF T>=1 THEN 1300 1290 X=1-X 1300 PRINT1400 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT R(4,I) 1410 NEXT I 1420 REM-ADJUST FOR NON STANDARD VALUES 1430 F "PERCENTILE =";1-X 1310 PRINT 1320 RETURN OR I=1 TO N 1440 R(1,I)=(R(4,I)-M)/S 1450 REM - COMPUTE FREQUENCY (Y COORDINATE 1460 R(2,I)=EXP(-R(1,I)^2/2)/2.50663 1470 Z=R(1,I) 1480 REM - APPROXIMATE PROBABILITY (AREA UNDER CURVE) 1490 R(3,I)=1/(1+.33267*ABS(R(1,I))) 1500 R(3,I)=.5-R(2,I)*(.4361000 A=SYS(7,1) \ REM - "POISSON DISTRIBUTION" 1010 PRINT "CALCULATED FREQUENCY"; 1020 INPUT A(1) 1030 IF A(1)<=0 THEN PRINT184*R(3,I)-.120168*R(3,I)^2+.937298*R(3,I)^3) 1510 IF R(1,I)>=0 THEN R(3,I)=R(3,I)+.5 \ GO TO 1530 1520 R(3,I)=.5-R(3,I) 1530 "FREQUENCY MUST BE >0" \ GO TO 1010 1040 PRINT "TEST FREQUENCY"; 1050 INPUT A(2) 1060 IF A(2)<0 THEN PRINT "TEST FREQUENCY M NEXT I 1540 PRINT "X-VALUE";TAB(10);"FREQUENCY(Y-VALUE)";TAB(32);"PROBABILITY" 1550 FOR I=1 TO N 1560 PRINT R(4,I);TAB(10);RUST BE >= 0" \ GO TO 1040 1070 REM - COMPUTE FACTORIAL 1080 A=0 1090 FOR I=1 TO A(2) 1100 A=A+LOG(I) 1110 NEXT I 1120 REM (2,I);TAB(32);R(3,I) 1570 NEXT I 1580 PRINT 1590 RETURN - COMPUTE PROBABILITY 1130 A=EXP(-A(1)+A(2)*LOG(A(1))-A) 1140 PRINT "PROBABILITY OF";A(2);"OCCURRENCES =";A 1150 PRINT 1160R8;") ? "; 1200 INPUT #I9,N \ IF I9<>0 THEN PRINT N 1210 IF N>0 THEN IF N<=R8 GO TO 1230 1220 IF I9=0 GO TO 1190 \ PRINT "ILLEGAL INPUT" \ CLOSE #I9 \ GO TO 1590 1230 FOR I=1 TO N 1240 PRINT "VALUE ";I;" = ? "; 1250 INPUT #I9,R(4,I) \ IF I9<>0 THEN P1000 A=SYS(7,1) \ REM - "F-DISTRIBUTION" 1010 PRINT "F-VALUE"; 1020 INPUT T 1030 PRINT "DEGREES OF FREEDOM IN NUMERATOR"; 10RINT R(4,I) 1260 NEXT I \ IF I9<>0 THEN CLOSE #I9 1270 REM-CHECK THE DATA 1280 PRINT "DO YOU WISH TO CHECK/CHANGE THE DATA"; 40 INPUT D1 1050 PRINT "DEGREES OF FREEDOM IN DENOMINATOR"; 1060 INPUT D2 1070 X=1 1080 REM - COMPUTE USING INVERSE FOR SMAL\ INPUT Y$ 1290 IF SEG$(Y$,1,1)<>"Y" GO TO 1420 1300 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK " 1310 IF M=0 TL T-VALUES 1090 IF T<1 THEN 1140 1100 S=D1 1110 R=D2 1120 Z=T 1130 GO TO 1170 1140 S=D2 1150 R=D1 1160 Z=1/T 1170 J=2/9HEN GO TO 1340 1320 PRINT "MEAN AND STANDARD DEVIATION = ";M;",";S; \ INPUT Y$ 1330 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE/S 1180 K=2/9/R 1190 REM - COMPUTE USING APPROXIMATION FORMULAS 1200 L=ABS((1-K)*Z^(1/3)-1+J)/SQR(K*Z^(2/3)+J) 1210 IF R<4 T  NS=";N; \ INPUT Y$ 1330 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT N 1340 IF K=1 THEN GO TO 1410 1350 PRINT "OBSERVATION","ITEM","FREQUENCY" 1360 FOR I=1 TO N 1370 PRINT I,R(1,I),R(2,I); \ INPUT Y$ 1380 IF SEG$(Y$,1,1)="Y" THEN PRINT "NE1000 A=SYS(7,1) \ REM - "MEAN, VARIANCE, STANDARD DEVIATION" 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1W VALUES = "; \ INPUT R(1,I),R(2,I) 1390 NEXT I 1400 GO TO 1450 1410 FOR I=1 TO N 1420 PRINT "ITEM ";I;"= ";R(1,I); \ INPUT ,1)<>"Y" THEN I9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "WHICHY$ 1430 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT R(1,I) 1440 NEXT I 1450 REM 1460 R=0 1470 M=0 1480 P=0 149 METHOD (0=POPULATION,1=SAMPLE) ? "; 1060 INPUT #I9,S \ IF I9<>0 THEN PRINT S 1070 IF S>=0 THEN IF S<=1 GO TO 1090 1080 IF I90 IF K=1 THEN 1640 1500 REM - FOR GROUPED DATA 1510 FOR I=1 TO N 1520 REM - ACCUMULATE ENTERED VALUES 1530 R=R+R(2,I)*R(1,I)=0 GO TO 1050 \ CLOSE #I9 \ PRINT "ILLEGAL INPUT" \ GO TO 1790 1090 PRINT "KIND OF DATA (0=GROUPED,1=UNGROUPED) ? "; 1100 INPU 1540 REM - CALCULATE INTERMEDIATE VALUES FOR VARIANCE 1550 P=P+R(2,I) 1560 M=M+R(2,I)*R(1,I)^2 1570 NEXT I 1580 REM - CALCT #I9,K \ IF I9<>0 THEN PRINT K 1110 IF K>=0 THEN IF K<=1 GO TO 1130 1120 IF I9=0 GO TO 1090 \ CLOSE #I9 \ PRINT "ILLEGAL INPUULATE MEAN AND VARIANCE 1590 R=R/P 1600 V=(M-P*R^2)/(P-S) 1610 REM - PRINT RESULTS 1620 GO TO 1730 1630 REM - FOR UNGROUPEDT" \ GO TO 1790 1130 PRINT "NUMBER OF OBSERVATIONS (MAXIMUM=";R8;") ? "; 1140 INPUT #I9,N \ IF I9<>0 THEN PRINT N 1150 IF N>0 DATA 1640 FOR I=1 TO N 1650 REM - ACCUMULATE ENTERED VALUES 1660 P=P+R(1,I) 1670 REM - ACCUMULATE INTERMEDIATE VALUES FOR V THEN IF N<=R8 GO TO 1170 1160 IF I9=0 GO TO 1130 \ PRINT "ILLEGAL INPUT" \ CLOSE #I9 \ GO TO 1790 1170 IF K=1 THEN GO TO 123ARIANCE 1680 M=M+R(1,I)^2 1690 NEXT I 1700 REM - CALCULATE MEAN AND VARIANCE, PRINT 1710 R=P/N 1720 V=(M-N*R^2)/(N-S) 17300 1180 FOR I=1 TO N 1190 PRINT "ITEM, FREQUENCY ";I;"? "; 1200 INPUT #I9,R(1,I),R(2,I) \ IF I9<>0 THEN PRINT R(1,I);" ";R(2, REM 1740 PRINT 1750 REM - PRINT RESULTS 1760 PRINT "MEAN","VARIANCE","STANDARD DEVIATION" 1770 PRINT R,V,SQR(V) 1780 PRINI) 1210 NEXT I 1220 GO TO 1270 1230 FOR I=1 TO N 1240 PRINT "ITEM ";I;"? "; 1250 INPUT #I9,R(1,I) \ IF I9<>0 THEN PRINT R(1T 1790 RETURN  RETURN ,I) 1260 NEXT I 1270 REM 1280 IF I9<>0 THEN CLOSE #I9 1290 PRINT "DO YOU WISH TO CHECK\CHANGE THE DATA "; \ INPUT Y$ 1300 IF SEG$(Y$,1,1)<>"Y" GO TO 1460 1310 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1320 PRINT "NUMBER OF OBSERVATIO  ); 1310 INPUT Y$ 1320 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUES = "; \ INPUT R(1,I),R(2,I) 1330 NEXT I 1340 REM-CALCULATION 1350 Z=0 1360 REM - INCLUDE EACH COMPONENT IN RELIABILITY 1370 FOR I=1 TO N 1380 Z=Z+1/R(1,I)+R(2,I) 1390 NEXT I 1400 PRINT 1410 REM - CALCUALTE RELIABILITY, PRINT 1420 Z=EXP(-Z*T) 1430 PRINT "SYSTEM RELIABILITY =";Z 1440 PRINT 1450 RETURN 1000 A=SYS(7,1) \ REM SYSTEM RELIABILITY 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT 1060 PRINT "OPERATING TIME IN HOURS ? "; 1070 INPUT #I9,T \ IF I9<>0 THEN PRINT T 1080 PRINT "NUMBER OF COMPONENTS (MAXIMUM=";R8;")? "; 1090 INPUT #I9,N \ IF I9<>0 THEN PRINT N 1100 IF N>0 THEN IF N<=R8 GO TO 1120 1110 IF I9=0 GO TO 1080 \ PRINT "ILLEGAL DATA" \ CLOSE #I9 \ G1000 A=SYS(7,1) \ REM - "PERMUTATIONS AND COMBINATIONS" 1010 PRINT "TOTAL NUMBER OF OBJECTS"; 1020 INPUT N 1030 PRINT "SIZE O TO 1450 1120 REM - ENTER DATA FOR EACH COMPONENT 1130 FOR I=1 TO N 1140 PRINT "COMPONENT";I 1150 PRINT " AVERAGE WEAROUT TOF SUBGROUP"; 1040 INPUT D 1050 REM - SIZE OF SUBGROUP CANNOT BE LARGER THAN SIZE OF GROUP 1060 IF D<=N THEN GO TO 1100 107IME ? "; 1160 INPUT #I9,R(1,I) \ IF I9<>0 THEN PRINT R(1,I) 1170 PRINT " AVERAGE FAILURE RATE ?"; 1180 INPUT #I9,R(2,I) \ IF 0 PRINT "SUBGROUP TOO LARGE" 1080 PRINT 1090 GO TO 1010 1100 REM 1110 REM - COMPUTE PERMUTATIONS 1120 P=1 1130 C=1 1140I9<>0 THEN PRINT R(2,I) 1190 NEXT I 1200 IF I9<>0 THEN CLOSE #I9 1210 REM CHECK 1220 PRINT "DO YOU WISH TO CHECK/CHANGE THE  FOR I=N-D+1 TO N 1150 REM - DON'T ALLOW NUMBER SIZE TO OVERFLOW MACHINE CAPACITY 1160 REM IF P>9.9*10^38/I THEN PRINT "MORE TDATA"; \ INPUT Y$ 1230 IF SEG$(Y$,1,1)<>"Y" GO TO 1340 1240 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK " 1250 PHAN 9.9E38 PERMUTATIONS" \ REM RETURN 1170 P=P*I 1180 NEXT I 1190 REM - COMPUTE INTERMEDIATE FACTORIAL FOR COMBINATIONS 1RINT "NUMBER OF COMPONENTS =";N; \ INPUT Y$ 1260 IF SEG$(Y$,1,1)<>"Y" GO TO 1280 1270 PRINT "NEW NUMBER OF COMPONENTS = "; \ I200 FOR J=2 TO D 1210 C=C*J 1220 NEXT J 1230 REM 1240 PRINT P;"PERMUTATIONS" 1250 PRINT P/C;"COMBINATIONS" 1260 PRINT 12NPUT N 1280 FOR I=1 TO N 1290 PRINT "COMPONENT ";I 1300 PRINT "AVERAGE WEAROUT TIME =";R(1,I);" AVERAGE FAILURE RATE =";R(2,I70 RETURN  F SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT D(I) 1290 NEXT I 1300 REM - ITERATIVELY COMPUTE MEAN 1310 FOR I=1 TO N 300 1200 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1210 FOR I=1 TO 2 1220 PRINT "SAMPLE ";I;", SIZE = ";T(I) 1320 M=M*D(I)^P 1330 REM - ACCUMULATE INTERMEDIATE TERM FOR DEVIATION 1340 Q=Q+LOG(D(I))^2 1350 NEXT I 1360 REM - COMPUTE D; \ INPUT Y$ 1230 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT T(I) 1240 FOR J=1 TO T(I) 1250 PRINT "ITEM ";J;" = "EVIATION 1370 IF N=1 THEN R=1 \ GO TO 1390 1380 R=EXP(SQR(Q/(N-1)-(N/(N-1)*(LOG(M))^2))) 1390 PRINT "GEOMETRIC MEAN =";M 140;R(I,J); \ INPUT Y$ 1260 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT R(I,J) 1270 NEXT J 1280 NEXT I 1290 REM 1300 PRINT "GEOMETRIC DEVIATION =";R 1410 PRINT 1420 RETURN 0 REM -SORT EACH SAMPLE 1310 X=0 \ Y=0 1320 FOR I=1 TO 2 1330 FOR J=1 TO T(I) 1340 FOR K=1 TO T(I)-J 1350 C=R(I,K) 1360 D=1000 A=SYS(7,1) \ REM - GEOMETRIC MEAN AND DEVIATION 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "NUMBER OF OBSERVATIONS (MAXIMUM=";D8;") ? "; 1060 INPUT #I9,N \ IF I9<>0 THEN PRINT N 1070 IF N>0 THEN IF N<=D8 GO TO 1090 1080 IF I9=0 GO TO 1050 \ CLOSE #I9 \ PRINT "ILLEGAL INPUT" \ GO TO 1420 1090 REM - COMPUTE WHICH ROOT TO USE 1100 P=1/N 1110 M=1 \ Q=0 11201000 A=SYS(7,1) \ REM - "MANN-WHITNEY U-TEST" 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I FOR I=1 TO N 1130 PRINT "ITEM";I;"? "; 1140 INPUT #I9,D(I) \ IF I9<>0 THEN PRINT D(I) 1150 IF D(I)<>0 GO TO 1180 1160 PRINT9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 REM - INPUT THE TWO SAMPLES "ZERO INPUT GIVES ZERO MEAN AND INFINITE DEVIATION" 1170 IF I9=0 GO TO 1130 \ CLOSE #I9 \ GO TO 1420 1180 NEXT I 1190 IF I9< 1060 FOR I=1 TO 2 1070 PRINT "SAMPLE";I;":" 1080 PRINT " SIZE (MAXIMUM=";R8;") ? "; 1090 INPUT #I9,T(I) \ IF I9<>0 THEN PRI>0 THEN CLOSE #I9 1200 PRINT "DO YOU WISH TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1210 IF SEG$(Y$,1,1)<>"Y" GO TO 1300 1220 PRINT T(I) 1100 IF T(I)>0 THEN IF T(I)<=R8 GO TO 1120 1110 IF I9=0 GO TO 1080 \ PRINT "ILLEGAL INPUT" \ CLOSE #I9 \ GO TO 1820 1NT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK " 1230 PRINT "NUMBER OF OBSERVATIONS = ";N; \ INPUT Y$ 1240 IF SEG$(Y$,1120 FOR J=1 TO T(I) 1130 PRINT " DATA ";J;" ? "; 1140 INPUT #I9,R(I,J) \ IF I9<>0 THEN PRINT R(I,J) 1150 NEXT J 1160 NEXT I,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT N 1250 P=1/N 1260 FOR I=1 TO N 1270 PRINT "ITEM ";I;" = ";D(I); \ INPUT Y$ 1280 I 1170 IF I9<>0 THEN CLOSE #I9 1180 PRINT "DO YOU WISH TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1190 IF SEG$(Y$,1,1)<>"Y" GO TO 1  0 1700 X=X+R 1710 I=I+1 1720 R=R+1 1730 GO TO 1470 1740 REM -U1=NUMBER OF TIMES SAMPLE 1 SCORES PRECEED SAMPLE 2 SCORES 17 "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK " 1310 PRINT "NUMBER OF ROWS , NUMBER OF COLUMNS = ";R;",";C; \ INPUT Y$ 150 U1=T(1)*T(2)+T(1)*(T(1)+1)/2-X 1760 REM - U=NUMBER OF TIMES SAMPLE 2 SCORES PRECEED SAMPLE 1 SCORES 1770 U2=T(1)*T(2)+T(2)320 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT R,C 1330 FOR I=1 TO R 1340 PRINT "ROW ";I 1350 FOR J=1 TO C 1360 *(T(2)+1)/2-Y 1780 REM 1790 PRINT 1800 PRINT "FIRST SAMPLE PRECEEDING, U =";U1 1810 PRINT "SECOND SAMPLE PRECEEDING, U =";UPRINT "ELEMENT ";J;" = ";X(I,J); \ INPUT Y$ 1370 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE = "; \ INPUT X(I,J) 1380 NEXT J \ P2 1820 RETURN RINT 1390 NEXT I 1400 REM 1410 REM - ADD UP MARGINAL FREQUENCIES FOR EACH ROW 1420 L=0 \ M=1 1430 FOR I=1 TO R 1440 FOR J1000 A=SYS(7,1) \ REM CHI-SQUARE TEST-INPUT 1010 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9==1 TO C 1450 R(2,I)=R(2,I)+X(I,J) 1460 NEXT J 1470 L=L+R(2,I) 1480 NEXT I 1490 REM - ADD UP MARGINAL FREQUENCIES FOR EACH C0 \ GO TO 1050 1030 PRINT "FILE SPECIFICATION"; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 M=X7 \ IF R80 THEN PRINT R 1080 IF R>0 THEN IF R<=M GO TOTEST-PRINT 1570 PRINT "OBSERVED VALUE EXPECTED VALUE CHI^2 CONTRIBUTION" 1580 FOR I=1 TO C 1590 PRINT " COLUMN";I 1600 FOR JR(I,K+1) 1370 IF R(I,K)T(1) THEN 1660 1480 IF J>T(2) THEN 1700 1ER OF COLUMNS (MAXIMUM=";L;") ? "; 1120 INPUT #I9,C \ IF I9<>0 THEN PRINT C 1130 IF C>0 THEN IF C<=M GO TO 1150 1140 IF I9=0 490 IF R(1,I)T(1) THEN 1580 1560 IF R(1,I)<>R(1,I-1) THEN 1580 1570 I=I+1 \ T "ELEMENT";J;"? "; \ INPUT #I9,X(I,J) \ IF I9<>0 THEN PRINT X(I,J) 1190 NEXT J 1200 NEXT I \ IF I9<>0 THEN CLOSE #I9 1210 REGO TO 1610 1580 IF J>T(2) THEN 1630 1590 IF R(2,J)<>R(2,J-1) THEN 1630 1600 J=J+1 1610 R1=R1+R \ R=R+1 \ K=K+1 1620 GO TO 1M CLEAR THE MARGINAL FREQUENCY ARRAY 1220 FOR I=1 TO 2 1230 FOR J=1 TO R8 1240 R(I,J)=0 1250 NEXT J 1260 NEXT I 1270 REM C550 1630 X=X+(I-M)*R1/K 1640 Y=Y+(J-L)*R1/K 1650 GO TO 1470 1660 IF J>T(2) THEN 1750 1670 Y=Y+R 1680 J=J+1 1690 GO TO 172HECK DATA 1280 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1290 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1410 1300 PRINT  R OF DAYS FROM 1ST OF YEAR TO 1ST OF EACH MONTH 1410 REM - GET NUMBER OF DAYS FROM JAN 1 TO 1ST OF MONTH FROM DATA TABLE 1420 1000 A=SYS(7,1) \ REM - "DAYS BETWEEN TWO DATES" \ PRINT 1010 DATA 0,31,59,90,120,151,181,212,243,273,304,334 1020 FOR I=1A=X(M) 1430 REM - COMPUTE NUMBER OF DAYS FROM 0,0,0 TO M,D,Y 1440 A=A+Y*365+INT(Y/4)+D+1-INT(Y/100)+INT(Y/400) 1450 REM - POS TO 12 1030 READ X(I) 1040 NEXT I 1050 RESTORE 1060 PRINT "FIRST DATE"; \ INPUT M1,D1,Y1 \ IF Y1<1582 GO TO 1560 1070 PRINSIBLY A LEAP YEAR? 1460 IF INT(Y/4)=Y/4 THEN 1490 1470 GO TO 1570 1480 REM - CONTINUE TEST FOR LEAP YEAR 1490 IF Y/400=INT(YT "SECOND DATE"; \ INPUT M2,D2,Y2 \ IF Y2<1582 GO TO 1560 1080 REM - SET VARIABLES TO BE USED IN SUBROUTINE 1090 M=M1 \ D=D1 \/400) THEN 1530 1500 IF Y/100=INT(Y/100) THEN 1570 1510 REM - YEAR IS A LEAP YEAR; 1520 REM - IF MONTH IS JAN OR FEB ADJUST C Y=Y1 1100 GOSUB 1220 1110 REM - SAVE COMPUTED NUMBER OF DAYS IN N 1120 REM - SET VARIABLES TO BE USED IN SUBROUTINE 1130 N=ALCULATED NUMBER OF DAYS 1530 IF M>2 THEN 1570 1540 A=A-1 \ GO TO 1570 1550 REM - END OF SUBROUTINE, RETURN TO MAIN PROGRAM A \ M=M2 \ D=D2 \ Y=Y2 1140 GOSUB 1220 1150 REM - CALCULATE DIFFERENCE AND PRINT 1160 N=A-N 1170 PRINT "DIFFERENCE =";N;"DAY1560 PRINT "YEAR IS PRIOR TO START OF GREGORIAN CALENDAR - 1582" 1570 RETURN S" \ PRINT 1180 GO TO 1570 1190 REM - SUBROUTINE TO COMPUTE NUMBER OF DAYS SINCE 0,0,0 TO M,D,Y 1200 REM - START WITH TEST FOR UNREAL DATE 1210 REM - GO TO CORRECT TEST DEPENDING ON NUMBER OF DAYS IN MONTH 1220 ON M GO TO 1280,1310,1280,1390,1280,1390,1280,1280,1390,1280,1390,1280 1230 REM - IF THIS MESSAGE IS PRINTED THE ANSWER IS PROBABLY INCORRECT 1240 PRINT "UNREAL DATE1000 A=SYS(7,1) \ REM - DAY OF THE WEEK 1010 DATA 0,3,3,6,1,4,6,2,5,0,3,5 1020 FOR I=1 TO 12 1030 READ X(I) 1040 NEXT I \RES" 1250 REM - STOP CALCULATIONS, RETURN TO MAIN PROGRAM 1260 RETURN 1270 REM - MONTH HAS 31 DAYS 1280 IF D>31 THEN 1240 1290TORE 1050 PRINT \ PRINT "ENTER DATE AS MONTH,DAY,YEAR"; \ INPUT M,D,Y 1060 I1=INT((Y-1500)/100) 1070 IF Y-1582<0 THEN 1330 =1 TO R 1610 REM - P=EXPECTED CELL VALUE 1620 P=R(2,J)*R(1,I)/L 1630 REM - Y=CHI-SQUARE CONTRIBUTION FROM THIS CELL 1640 Y=( GO TO 1410 1300 REM - MONTH IS FEBRUARY; A LEAP YEAR? 1310 IF Y/4=INT(Y/4) THEN IF Y/400=INT(Y/400) THEN 1360 1320 IF Y/100=X(J,I)-P)^2/P 1650 REM - Z=TOTAL CHI-SQUARE VALUE 1660 Z=Z+Y 1670 PRINT " ";X(J,I),P,Y 1680 NEXT J 1690 NEXT I 1700 PRINTINT(Y/100) THEN 1340 1330 REM - NOT A LEAP YEAR; MONTH HAS 28 DAYS 1340 IF D>28 THEN 1240 1350 REM - A LEAP YEAR; MONTH HAS 2 1710 PRINT "CHI-SQUARE =";Z 1720 PRINT "DEGREES OF FREEDOM =";(C-1)*(R-1) 1730 RETURN 9 DAYS 1360 IF D>29 THEN 1240 1370 GO TO 1410 1380 REM - MONTH HAS 30 DAYS 1390 IF D>30 THEN 1240 1400 REM - TABLE OF NUMBE 1080 A=I1*5+(I1+3)/4 1090 I2=INT(A-INT(A/7)*7) 1100 Y2=INT(Y/100) 1110 Y3=INT(Y-Y2*100) 1120 A=Y3/4+Y3+D+X(M)+I2 1130 B=INTRINT 1120 REM - DEPRECIATION RATE CONVERTED TO DECIMAL FOR USE IN CALCULATIONS 1130 PRINT "VALUE = $";INT(100*P*(1-I/100)^Y+.5(A-INT(A/7)*7)+1 1140 IF M>2 THEN 1220 1150 IF Y3=0 THEN 1200 1160 T1=INT(Y-INT(Y/4)*4) 1170 IF T1<>0 THEN 1220 1180 IF B=0)/100 1140 PRINT 1150 REM - RETURN FOR NEXT YEAR NUMBER 1160 GO TO 1070 1170 RETURN  THEN B=6 1190 B=B-1 \ GO TO 1220 1200 A=I1-1 \ T1=INT(A-INT(A/4)*4) 1210 IF T1=0 THEN 1180 1220 IF B=0 THEN B=7 1230 PRINT M;",";D;",";Y;" A "; 1240 ON B GO TO 1250,1260,1270,1280,1290,1300,1310 1250 PRINT "SUNDAY" \ GO TO 1340 1260 PRINT "MONDAY" \ GO TO 1340 1270 PRINT "TUESDAY" \ GO TO 1340 1280 PRINT "WEDNESDAY" \ GO TO 1340 1290 PRINT "THURSDAY" \ GO TO 1340 1300 PRINT "FRIDAY" \ GO TO 1340 1310 PRINT "SATURDAY" \ GO TO 1340 1320 PRINT "ILLEGAL DATE " \ GO TO 1340 1330 PRINT "DATE ENTER1000 A=SYS(7,1) \ REM - "FUTURE VALUE OF AN INVESTMENT" 1010 PRINT 1020 PRINT "INITIAL INVESTMENT"; 1030 INPUT P 1040 PRINED IS PRIOR TO START OF GREGORIAN CALENDAR - 1582 " 1340 RETURN T "NOMINAL INTEREST RATE"; 1050 INPUT I 1060 PRINT "NUMBER OF COMPOUNDING PERIODS PER YEAR"; 1070 INPUT N 1080 PRINT "NUMBER OF YEARS"; 1090 INPUT Y 1100 REM - CALCULATE INTEREST RATE PER PERIOD, 1110 REM - CONVERT FROM PERCENT TO DECIMAL 1120 I=I/N/100 1130 REM _ CALCULATE FUTURE VALUE BY FORMULA 1140 T=P*(1+I)^(N*Y) 1150 REM - ROUND OFF TO NEAREST CENT PRINT 1160 PRINT "FUTURE VALUE=$",INT(T*100+.5)/100 1170 REM - PRINT BLANK LINE TO SEPARATE DATA FROM QUESTION 1180 PRINT 1190 RETURN 1000 A=SYS(7,1) \ REM - "SALVAGE VALUE" 1010 PRINT 1020 PRINT "ORIGINAL PRICE"; 1030 INPUT P 1040 PRINT "DEPRECIATION RATE"; 1050 INPUT I 1060 PRINT "--(ENTER YEARS=0 WHEN NO MORE VALUES DESIRED FOR THIS ITEM)--" 1070 PRINT "YEARS"; 1080 INPUT Y  1090 REM - CALCULATE ANOTHER SALVAGE VALUE? 1100 IF Y=0 THEN 1170 1110 REM - CALCULATE SLAVAGE VALUE BY FORMULA, ROUND OFF, P  RATE"; 1090 INPUT I 1100 REM - CALCULATE INTEREST RATE PER PERIOD, 1110 REM - CONVERT FROM % TO DECIMAL 1120 I=I/N/100 1130 P=R*N*(1-1/((I/100)/N+1)^(N*Y))/(I/100) 1130 REM - ROUND OFF TO NEAREST CENT, PRINT 1140 PRINT "PRINCIPAL = $";INT(P*100+.5)0 REM - CALCULATE INITIAL INVESTMENT BY FORMULA 1140 P=T/(1+I)^(N*Y) 1150 REM - ROUND OFF TO NEAREST CENT PRINT 1160 PRINT "/100 1150 REM - PRINT BLANK LINE TO SEPERATE DATA FROM QUESTION 1160 PRINT 1170 RETURN INITIAL INVESTMENT = $";INT(P*100+.5)/100 1170 REM - PRINT BLANK LINE TO SEPARATE DATA FROM QUESTION 1180 PRINT 1190 RETURN 1000 A=SYS(7,1) \ REM - "REGULAR DEPOSITS" 1010 PRINT 1020 PRINT "TOTAL VALUE AFTER Y YEARS"; 1030 INPUT T 1040 PRINT "NOM INAL INTEREST RATE"; 1050 INPUT I 1060 PRINT "NUMBER OF DEPOSITS PER YEAR"; 1070 INPUT N 1080 PRINT "NUMBER OF YEARS"; 1090 INPUT Y 1100 REM - CALCULATE INTEREST RATE PER DEPOSIT, 1110 REM - CONVERT FROM PERCENT TO DECIMAL 1120 I=I/N/100 1130 REM _ CALCULATE AMOUNT OF REGULAR DEPOSIT BY FORMULA 1140 R=T*I/((I+1)^(N*Y)-1) 1150 REM - ROUND OFF TO NEAREST CENT, PRINT 1160 1000 A=SYS(7,1) \ REM - "NOMINAL INTEREST RATE ON INVESTMENTS" 1010 PRINT 1020 PRINT "PRINCIPAL"; 1030 INPUT P 1040 PRINT PRINT "REGULAR DEPOSITS = $";INT(R*100+.5)/100 1170 REM - PRINT BLANK LINE TO SEPARATE DATA FROM QUESTION 1180 PRINT 1190 RE"TOTAL VALUE"; 1050 INPUT T 1060 PRINT "NUMBER OF YEARS"; 1070 INPUT Y 1080 PRINT "NUMBER OF COMPOUNDING PERIODS PER YEAR"; TURN  1090 INPUT N 1100 REM - CALCULATE NOMINAL INTEREST RATE BY FORMULA, PRINT 1110 I2=N*((T/P)^(1/(N*Y))-1)*100 1120 PRINT "NOMINAL INTEREST RATE=";I2;"%" 1130 REM - PRINT BLANK LINE TO SEPARATE DATA FORM QUESTION 1140 PRINT 1150 RETURN 1000 A=SYS(7,1) \ REM - "PRINCIPAL ON A LOAN" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; 1030 INPUT R 1040 PRINT "TERM IN YE1000 A=SYS(7,1) \ REM - "INITIAL INVESTMENT" 1010 PRINT 1020 PRINT "TOTAL VALUE AFTER Y YEARS"; 1030 INPUT T 1040 PRINT "NARS"; 1050 INPUT Y 1060 PRINT "ANNUAL INTEREST RATE"; 1070 INPUT I 1080 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1090 INPUT N UMBER OF COMPOUNDING PERIODS PER YEAR"; 1050 INPUT N 1060 PRINT "NUMBER OF YEARS"; 1070 INPUT Y 1080 PRINT "NOMINAL INTEREST1100 REM - CALCULATE AMOUNT OF PRINCIPAL BY FORMULA 1110 REM - INTEREST CONVERTED FROM PERCENT TO DECIMAL FOR CALCULATIONS 112 ANCE REMAINING DECREASES WITH EACH PAYMENT 1220 B0=B0-A 1230 NEXT J1 1240 REM - CALCULATE LAST PAYMENT,ROUND OFF, PRINT 1250TS DESIRED FOR THIS ITEM)--" 1090 PRINT "YEAR"; 1100 INPUT Y 1110 REM - THROUGH CALCULATING FOR THIS ITEM? 1120 IF Y=0 THEN  PRINT "LAST PAYMENT=$";INT((R+B0)*100+.5)/100 1260 PRINT 1270 RETURN 1200 1130 REM - CALCULATE DEPRECIATION AMOUNT BY FORMULA 1140 D=P*I*(1-I)^(Y-1) 1150 REM - ROUND OFF TO NEAREST CENT, PRINT 1160 PRINT "DEPRECIATION = $";INT(D*100+.5)/100 1170 PRINT 1180 REM - RETURN FOR NEXT YEAR NUMBER 1190 GO TO 1090 1200 RETU1000 A=SYS(7,1) \ REM - "TERM OF A LOAN" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; 1030 INPUT R 1040 PRINT "PRINCIPAL"; 10RN 50 INPUT P 1060 PRINT "ANNUAL INTEREST RATE"; 1070 INPUT I 1080 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1090 INPUT N 1100 REM - CALCULATE TERM IN YEARS BY FORMULA, 1110 REM - INTEREST CONVERTED FROM PERCENT TO DECIMAL FOR CALCULATION 1120 Y=-(LOG(1-(P*(I/100))/(N*R))/(LOG(1+I/100/N)*N)) 1130 REM - ROUND OFF TO NEAREST TENTH, PRINT 1140 PRINT "TERM =";INT(Y*10+.5)/10;"YEARS" 1150 PRINT 1160 RETURN 1000 A=SYS(7,1) \ REM - "LAST PAYMENT ON A LOAN" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; 1030 INPUT R 1040 PRINT "PRINCIPAL"; 1050 INPUT P 1060 PRINT "TERM IN YEARS"; 1070 INPUT Y 1080 PRINT "ANNUAL INTEREST RATE"; 1090 INPUT I 1100 REM - CONVERT INTEREST FROM PERCENT TO DECIMAL 1110 I=I/100 1120 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1130 INPUT N 1140 B0=P 1150 REM - COMPUTE ALL PAYMENTS, BALANCES THROUGH LAST PAYMENT USING R 1160 FOR J1=1 TO N*Y 1170 REM - ROUND OFF INTEREST PAID TO NEAR1000 A=SYS(7,1) \ REM - "DEPRECIATION AMOUNT" 1010 PRINT 1020 PRINT "ORIGINAL PRICE"; 1030 INPUT P 1040 PRINT "DEPRECIATIOEST CENT 1180 I1=INT((B0*I/N)*100+.5)/100 1190 REM - CALCULATE AMOUNT AMORTIZED WITH EACH PAYMENT 1200 A=R-I1 1210 REM - BALN RATE"; 1050 INPUT I 1060 REM - CONVERT FROM PERCENT TO DECIMAL 1070 I=I/100 1080 PRINT "--(ENTER YEAR=0 WHEN NO MORE AMOUN AREST CENT 1170 R1=INT(R1*100+.5)/100 1180 REM - I3=NUMBER USED TO CLOSE IN ON INTEREST RATE 1190 I3=ABS(I-I2)/2 1200 REM - 00 I1=INT((B0*I/N)*100+.5)/100 1210 REM - CALCULATE AMOUNT AMORTIZED WITH EACH PAYMENT 1220 A=R-I1 1230 REM - CALCULATE REMAISAVE THIS GUESS 1210 I2=I 1220 REM - COMPARE COMPUTED PAYMENT (R1) TO INPUT PAYMENT (R) 1230 REM - IF THEY'RE EQUAL, LAST RATNING BALANCE ON PRINCIPAL 1240 B0=B0-A 1250 NEXT J1 1260 REM - ROUND OFF, PRINT 1270 PRINT "REMAINING BALANCE = $";INT(B0*10E GUESSED=APPROXIMATE REG. PAYMENT 1240 IF R1=R THEN 1350 1250 IF R1>R THEN 1310 1260 REM - R1R, RATE MUST BE LOWER THAN LAST GUESS 11000 A=SYS(7,1) \ REM "REGULAR PAYMENT ON A LOAN" 1010 PRINT 1020 PRINT "TERM OF YEARS"; 1030 INPUT Y 1040 PRINT "PRINCIPAL310 I=I-I3 1320 REM - RETEST WITH NEW GUESS 1330 GO TO 1150 1340 REM - COMPUTE INTEREST TO PROPER PROPORTIONS, ROUND OFF, PRI"; 1050 INPUT P 1060 PRINT "ANNUAL INTEREST RATE"; 1070 INPUT I 1080 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1090 INPUT N 110NT 1350 I=((INT((I*1000)*100+.5))/100)/1000 1360 PRINT "ANNUAL INTEREST RATE =";I*100;"%" 1370 PRINT 1380 RETURN 0 REM - CALCULATE AMOUNT OF REGULAR PAYMENT BY FORMULA, 1110 REM - INTEREST CONVERTED FROM PERCENT TO DECIMAL FOR CALCULATIONS  1120 R=((I/100)*P/N)/(1-1/((I/100)/N+1)^(N*Y)) 1130 REM - ROUND OFF TO NEAREST CENT, PRINT 1140 PRINT "REGULAR PAYMENT = $";INT(R*100+.5)/100 1150 REM - PRINT BLANK LINE TO SEPARATE DATA FORM QUESTION 1160 PRINT 1170 RETURN 1000 A=SYS(7,1) \ REM - "REMAINING BALANCE ON A LOAN" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; 1030 INPUT R 1040 PRINT "PR1000 A=SYS(7,1) \ REM - "ANNUAL RATE ON A LOAN" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; 1030 INPUT R 1040 PRINT "TERM IN INCIPAL"; 1050 INPUT P 1060 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1070 INPUT N 1080 PRINT "ANNUAL INTEREST RATE"; 1090 INPUTYEARS"; 1050 INPUT Y 1060 PRINT "PRINCIPAL"; 1070 INPUT P 1080 PRINT "NUMBER OF PAYMENTS PER YEAR"; 1090 INPUT N 1100 REM  I 1100 REM - CONVERT FROM PERCENT TO DECIMAL 1110 I=I/100 1120 REM - ENTER THE PAYMENT NUMBER WITHIN THE YEAR, I.E. N1<=N - GUESS AN INTEREST RATE (10%) TO INITIATE TESTING 1110 I=10 1120 REM - I2=LAST GUESS OR ESTIMATE (START WITH 0) 1130 I2=0 11130 PRINT "LAST PAYMENT MADE (PAYMENT NO, YEAR)"; 1140 INPUT N1,Y 1150 REM - INITIALIZE REMAINING BALANCE 1160 B0=P 1170 RE140 REM - COMPUTE REGULAR PAYMENT USING GUESSED INTEREST RATE 1150 R1=(I*P/N)/(1-1/((I/N+1)^(N*Y))) 1160 REM - ROUND OFF TO NEM - LOOP TO ACCUMULATE AMOUNT PAID SO FAR 1180 FOR J1=1 TO N*(Y-1)+N1 1190 REM - CALCULATE INTEREST PAID WITH EACH PAYMENT 12 CENT, PRINT 1160 PRINT "AMOUNT OF EACH WITHDRAWAL = $";INT(R*100+.5)/100 1170 REM - PRINT BLANK LINE TO SEPARATE QUESTION FROM1000 A=SYS(7,1) \ REM - "FUTURE VALUE OF REGULAR DEPOSITS (ANNUITY)" 1010 PRINT 1020 PRINT "AMOUNT OF REGULAR DEPOSITS"; 10 DATA 1180 PRINT 1190 RETURN 30 INPUT R 1040 PRINT "NOMINAL INTEREST RATE"; 1050 INPUT I 1060 PRINT "NUMBER OF DEPOSITS PER YEAR"; 1070 INPUT N 1080 PRINT "NUMBER OF YEARS"; 1090 INPUT Y 1100 REM - CALCULATE INTEREST RATE PER DEPOSIT, 1110 REM - CONVERT FROM PERCENT TO DECIMAL 1120 I=I/N/100 1130 REM - CLACULATE FUTURE VALUE BY FORMULA 1140 T=R*((1+I)^(N*Y)-1)/I 1150 REM - ROUND OFF TO NEAREST CENT1000 A=SYS(7,1) \ REM - "MINIMUM INVESTMENT FOR WITHDRAWALS" 1010 PRINT 1020 PRINT "AMOUNT OF WITHDRAWALS"; 1030 INPUT R 1, PRINT 1160 PRINT "FUTURE VALUE = $";INT(T*100+.5)/100 1170 REM - PRINT BLANK LINE TO SEPARATE DATA FROM QUESTION 1180 PRINT040 PRINT "NOMINAL INTEREST RATE"; 1050 INPUT I 1060 PRINT "NUMBER OF WITHDRAWALS PER YEAR"; 1070 INPUT N 1080 PRINT "NUMBER 1190 RETURN  OF YEARS"; 1090 INPUT Y 1100 REM - CONVERT FROM PERCENT TO DECIMAL 1110 I=I/100 1120 REM - CALCULATE MINIMUM INVESTMENT BY FORMULA 1130 P=R*N/I*(1-1/((1+I/N)^(N*Y))) 1140 REM - ROUND OFF TO NEAREST CENT, PRINT 1150 PRINT "MINIMUM INVESTMENT = $";INT(100*P+.5)/100 1160 REM - PRINT BLANK LINE TO SEPARATE DATA FROM QUESTION 1170 PRINT 1180 RETURN 1000 A=SYS(7,1) \ REM - "REGULAR WITHDRAWALS FROM AN INVESTMENT" 1010 PRINT 1020 PRINT "INITIAL INVESTMENT"; 1030 INPUT P 1040 PRINT "NOMINAL INTEREST RATE"; 1050 INPUT I 1060 PRINT "NUMBER OF WITHDRAWALS PER YEAR"; 1070 INPUT N 1080 PRINT "NUMBER OF YEARS"; 1090 INPUT Y 1100 REM - CALCULATE INTEREST RATE PER WITHDRAWAL, 1110 REM - CONVERT FROM PERCENT TO DECIMAL 1120 I=I/N/100 1130 REM - CALCULATE REGULAR WITHDRAWAL BY FORMULA 1140 R=P*(I/((1+I)^(N*Y)-1)+I) 1150 REM - ROUND OFF TO NEAREST  AS FILE #1 \ INPUT #1,P9 80 PRINT "FUNCTION CODES" 90 FOR I=1 TO P9 100 INPUT #1,F$(0) 110 PRINT I;" = ";F$(0) 120 NEXT I \ 9000 1000 REM 9000 END  RESET #1 \ INPUT #1,P9 130 PRINT "FUNCTION CODE "; \ INPUT F9 140 IF F9>0 THEN IF F9<=P9 THEN P8=1 \ GO TO 160 150 PRINT "IL1000 A=SYS(7,1) \ REM - "EFFECTIVE INTEREST TATE ON INVESTMENTS" 1010 PRINT 1020 PRINT "INITIAL INVESTMENT"; 1030 INPUT P LEGAL FUNCTION CODE" \ GO TO 80 160 FOR I=1 TO F9 170 INPUT #1,F$(0),F$(1) 180 NEXT I \ CLOSE #1 190 PRINT "LOADING ";F$(0);1040 PRINT "TOTAL VALUE AFTER Y YEARS"; 1050 INPUT T 1060 PRINT "NUMBER OF YEARS"; 1070 INPUT Y 1080 REM - CALCULATE EFFECTI"...." 200 OVERLAY D9$+F$(1) 210 GOSUB 1000 220 PRINT "MORE DATA "; \ INPUT Y$ 230 IF SEG$(Y$,1,1)="Y" THEN GO TO 210 240 VE INTEREST RATE, PRINT 1090 PRINT "ANNUAL INTEREST RATE =";((T/P)^(1/Y)-1)*100;"%" 1100 REM - PRINT BLANK LINE TO SEPARATE DAPRINT "DO YOU WANT ANOTHER FUNCTION "; \ INPUT Y$ 250 IF SEG$(Y$,1,1)='Y' THEN GO TO 70 260 GO TO 9000 1000 REM 9000 END TA FROM QUESTION 1110 PRINT 1120 RETURN 10 DIM F$(2),A$(100),X(100),Y(100) 20 DEF FNC(X)=X 30 A=SYS(7,1) \ REM -ENSURE ALL REPLIES INTERPRETED AS UPPER CASE 40 DEF F1000 A=SYS(7,1) \ REM - "DEPRECIATION RATE" 1010 PRINT 1020 PRINT "ORIGINAL PRICE"; 1030 INPUT P 1040 PRINT "RESALE PRICE"NR(X)=INT(X*100+.5)/100 50 PRINT "DEVICE ON WHICH PACKAGE RESIDES "; \ INPUT D9$ 60 Y$=SEG$(D9$,1,4) 70 IF Y$<>" " THEN IF Y$; 1050 INPUT T 1060 PRINT "YEARS"; 1070 INPUT Y 1080 REM - CALCULATE DEPRECIATION RATE BY FORMULA CONVERT TO PERCENT 1090 D<>"SY:" THEN IF Y$<>"SY1:" THEN IF Y$<>"SY0:" GO TO 40 80 REM SET UP THE ARRAY OF FUNCTION NAMES 90 OPEN D9$+"FIN.DAT" FOR INP=100*(1-(T/P)^(1/Y)) 1100 REM - ROUND OFF, PRINT 1110 PRINT "DEPRECIATION RATE =";INT(1000*D+.5)/1000;"%" 1120 PRINT 1130 RUT AS FILE #1 \ INPUT #1,P9 100 PRINT "FUNCTION CODES " 110 FOR I=1 TO P9 120 INPUT #1,F$(0) 130 PRINT I;" = ";F$(0) 140 NEETURN XT I \ RESTORE #1 \ INPUT #1,P9 150 PRINT "FUNCTION CODE "; \ INPUT F9 160 IF F9>0 THEN IF F9<=P9 GO TO 190 170 IF F9=99 THEN10 DIM F$(1),T(2),A(3),S(25,2),R(4,25),D(25),C(2,2) \ R8=25 \ D8=25 \ S8=25 20 DIM X(10,10) \ X7=10 \ X8=10 \ REM X(X7,X8) 30 GO TO 9000 180 PRINT "ILLEGAL FUNCTION CODE" \ GO TO 100 190 FOR I=1 TO F9 200 INPUT #1,F$(0),F$(1) 210 NEXT I \ CLOSE #1  A=SYS(7,1) \ REM-ENSURE ALL REPLIES INTERPRETED AS UPPER CASE 40 PRINT "DEVICE ON WHICH PACKAGE RESIDES "; \ INPUT D9$ 50 Y$= 220 PRINT "LOADING ";F$(0);"...." 230 OVERLAY D9$+F$(1) 240 GOSUB 1000 250 PRINT "MORE DATA "; \ INPUT Y$ 260 IF SEG$(Y$,1,SEG$(D9$,1,4) 60 IF Y$<>" " THEN IF Y$<>"SY:" THEN IF Y$<>"SY0:" THEN IF Y$<>"SY1:" GO TO 40 70 OPEN D9$+"STAT.DAT" FOR INPUT 1)="Y" THEN GO TO 240 270 PRINT "DO YOU WANT ANOTHER FUNCTION "; \ INPUT Y$ 280 IF SEG$(Y$,1,1)='Y' THEN GO TO 90 290 GO TO EN B=A \ A=X \ REM MAKE A POSITIVE FNC VALUE 1170 X=(A+B)/2 \ REM MIDPOINT 1180 Y=SGN(FNC(X)) 1190 IF Y=0 THEN 1290 \ REM ROON OF THE CURRENT ROW 1300 Y=1/A(J,J) 1310 FOR K=0 TO N 1320 A(J,K)=Y*A(J,K) 1330 NEXT K 1340 REM - SUBTRACT THE CURRENT EQT FOUND 1200 IF Y<0 THEN B=X \ GO TO 1220 \ REM NEWLIMIT 1210 A=X 1220 IF ABS(B-A)/(ABS(A)+ABS(B))<5.00000E-06 THEN 1290 \ REUATION FROM 1350 REM - THE OTHER ROWS 1360 FOR I=1 TO N 1370 IF I=J THEN 1420 1380 Y=-A(I,J) 1390 FOR K=0 TO N 1400 A(I,K)M NEAR ROOT? 1230 GO TO 1170 1240 REM - ROOT AT AN INTERVAL LIMIT; FIND WHICH LIMIT, PRINT 1250 IF C=0 THEN 1280 1260 X=D 1=A(I,K)+Y*A(J,K) 1410 NEXT K 1420 NEXT I 1430 REM - THIS PROCESS IS REPEATED FOR ALL EQUATIONS 1440 NEXT J 1450 PRINT 146270 GO TO 1290 1280 X=C 1290 PRINT "ROOT = ";X 1300 RETURN 0 FOR J=N TO 1 STEP -1 1470 B(J,0)=A(J,0) 1480 FOR I=J+1 TO N 1490 B(J,0)=B(J,0)-A(J,I)*B(I,0) 1500 NEXT I 1510 B(J,0)=B(J,0)/A(J,J) 1520 NEXT J 1530 FOR I=1 TO N 1540 PRINT "X";I;"= ";INT(B(I,0)*1000+.5)/1000 1550 NEXT I 1560 RETURN 1000 A=SYS(7,1) \ REM SIMULTANEOUS EQUATIONS 1010 PRINT 1020 PRINT "NUMBER OF EQUATIONS"; \ INPUT N 1030 IF N>A8 THEN PRINT 1000 A=SYS(7,1) \ REM ROOTS OF POLYNOMIAL - HALF INTERVAL 1010 PRINT "FORMULA MUST BE ENTERED AT LINE 20 IN NUM.BAS, PROCEED ";"TOO MANY EQUATIONS" \ GO TO 1560 1040 PRINT "COEFFICIENT MATRIX" 1050 FOR J=1 TO N 1060 PRINT "EQUATION";J 1070 FOR I=1 TO  \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1300 1030 PRINT "INTERVAL (LOWER, UPPER LIMIT)"; \ INPUT A,B 1040 C=SGN(FNCN 1080 PRINT "COEFFICIENT";I; \ INPUT A(J,I) 1090 NEXT I 1100 PRINT "CONSTANT"; \ INPUT A(J,0) 1110 NEXT J 1120 FOR J=1 TO (A)) \ D=SGN(FNC(B)) 1050 IF C*D=0 THEN 1250 \ REM ROOT AT END 1060 IF C*D<0 THEN X=B \ GO TO 1160 \ REM OPPOSITE SIGNS AT ENDN 1130 REM - FIND THE FIRST EQUATION WITH A 1140 REM - NON-ZERO COEFFICIENT FOR THE CURRENT COLUMN 1150 K=-9999 \ L=J 1160 S 1070 REM - LOOP TO SEARCH 1000 NUMBERS FOR OPPOSITE SIGNS IN FUNCTION 1080 FOR I=1 TO 1000 1090 X=A+RND(2)*(B-A) \ Y=SGN(FNFOR I=J TO N 1170 IF A(I,J)=0 THEN 1190 1175 IF A(I,J)<=K THEN 1190 1180 L=I \ K=A(I,J) 1190 NEXT I 1200 IF K>-9999 THEN 12C(X)) 1100 IF Y=0 THEN 1290 \ REM ROOT FOUND 1110 IF C*Y<0 THEN 1150 \ REM OPPOSITE SIGNS AT RANDOM AND LOWER LIMIT 1120 NEXT30 1210 PRINT "NO UNIQUE SOLUTION" \ GO TO 1560 1220 REM - MOVE THAT EQUATION UP TO THE CURRENT ROW 1230 FOR K=0 TO N 1240 X I 1130 PRINT "NO CHANGE OF SIGN FOUND" \ GO TO 1300 1140 REM - CHANGE OF SIGN FOUND; CALCULATE ROOT 1150 B=X 1160 IF C<0 TH=A(J,K) 1250 A(J,K)=A(L,K) 1260 A(L,K)=X 1270 NEXT K 1280 REM - GET A 1 COEFFICIENT IN THE FIRST 1290 REM - NON-ZERO COLUM COMPOUNDING PERIODS PER YEAR"; \ INPUT N 1140 N1=0 \ L2=4 \ REM - PRINT FOUR TIMES PER YEAR 1150 PRINT "START WITH WHAT YEAR"; NEW BALANCE 1550 N2=N2+1 \ REM - COUNT DEPOSITS/WITHDRAWALS MADE PER YEAR 1560 B2=B0*(1+I/N) 1570 I1=B2-B0 \ REM - I1=AMOUNT \ INPUT X 1160 PRINT "END PRINTING WITH WHAT YEAR"; \ INPUT Y 1170 X=INT(X) \ REM - START PRINTING AT THE BEGINNING OF A YEAR INTEREST WITH EACH COMPOUNDING PERIOD 1580 I3=I3+I1 \ REM - I3=AMOUNT INTEREST ACCUMULATED BETWEEN POSTING 1590 I2=I2+I1 \ RE 1180 REM - INITIATE RUNNING TOTALS 1190 B0=P \ I1=0 \ I2=0 \ I3=0 \ K=66 \ P1=4 1200 PRINT "WHERE ARE THE RESULTS TO BE PRINM - I2=TOTAL INTEREST ACCUMULATED TO DATE 1610 IF P2/P1>J1/N THEN 1640 1620 P2=P2+1 1630 REM - YEAR TO START PRINTING? 1640 TED/SAVED - TT:,LP: OR FILESPEC "; 1210 INPUT P1$ \ IF SEG$(P1$,1,3)="TT:" THEN I9=0 \ GO TO 1240 1220 IF LEN(P1$)=0 THEN GOIF J0=Y THEN 1810 1750 NEXT J1 1760 REM - START PRINTING? 1770 IF J0N1 THEN 1560 1520 REM - TIME TO MAKE DEPOSIT/WITHDRAWAL? 1530 IF N2/N1>J1/N THEN 1560 1540 B0=B0+R \ REM - CALCULATE 1000 A=SYS(7,1) \ REM - "MORTGAGE AMORTIZATION TABLE" 1010 PRINT 1020 PRINT "REGULAR PAYMENT"; \ INPUT R1 1030 PRINT "TERM E SO FAR 1410 A=R-I1 \ REM - CALCULATE AMOUNT AMORTIZED THIS PAYMENT 1420 A1=A1+A \ REM - SUM AMOUNT AMORTIZED TO DATE 1430 BIN YEARS"; \ INPUT Y 1040 PRINT "PRINCIPAL"; \ INPUT P 1050 PRINT "ANNUAL INTEREST RATE"; \ INPUT I 1060 I=I/100 \ REM - CONV0=P-A1 \ REM - CALCULATE BALANCE DUE 1440 REM - LAST PAYMENT? IF YES, CALCULATE AMOUNT SO THAT THE 1450 REM - BALANCE DUE EQUAERT FROM PERCENT TO DECIMAL 1070 PRINT "NUMBER OF PAYMENTS PER YEAR"; \ INPUT N 1080 PRINT "START PRINTING WITH WHAT YEAR"; \ LS $00 00 AFTER THIS PAYMENT 1460 IF C1<>N*Y THEN 1480 1470 R=R+B0 \ A=A+B0 \ A1=A1+B0 \ B0=0 1480 I2=I2+I1 \ REM - SUM INTERINPUT X 1090 X=INT(X) \ REM - START PRINTING AT BEGINNING OF A YEAR 1100 REM - INITIALIZE VARIABLES 1110 C1=0 \ I2=0 \ I3=0 \EST PAID TO DATE 1490 I3=I3+I1 \ REM - SUM INTEREST PAID THIS YEAR 1500 A2=A2+A \ REM - SUM AMOUNT AMORTIZED THIS YEAR 1510 R J0=0 \ N1=N \ K=66 \ B0=P \ A1=0 \ A2=0 \ R=R1 1120 IF INT(Y)<1 THEN 1720 \ REM TERM LESS THEN ONE YEAR? 1130 PRINT "WHERE DOEM - STARTED PRINTING? IF YES, PRINT COMPUTED VALUES IN TABLE 1520 IF J0N*Y THEN 1590 1570 PRINT #I9," LAST PAYMENT = $";(INT(R*150 IF LEN(P1$)=0 THEN GO TO 1130 1160 OPEN P1$ FOR OUTPUT AS FILE #3 \ I9=3 1170 REM - LOOP FOR EACH YEAR 1180 FOR J0=1 TO I00+.5))/100 1580 REM - STARTED PRINTING? IF YES, PRINT YEARLY TOTALS 1590 IF J0Y THEN 1750 1650 REM - REINITIALIZE YEARLY VARIABLES 16 REM - SPACE TO TOP OF NEXT PAGE (ASSUME66LINESPERPAGE) 1240 FOR K1=K TO 66 1250 PRINT #I9, 1260 NEXT K1 1270 PRINT #I9, 1260 I3=0 1670 A2=0 1680 NEXT J0 1690 REM - NEED TO PRINT A PARTIAL YEAR? 1700 IF Y=J0 THEN 1750 1710 REM - ADJUST VARIABLES 80 REM - PRINT PAGE HEADINGS 1290 PRINT #I9," MORTGAGE AMORTIZATION TABLE" 1300 PRINT #I9," PRINCIPAL $";P;" AT";I*100;"% FORTO PRINT A PARTIAL YEAR 1720 N1=((Y-INT(Y))*12)/12*N 1730 J0=J0+1 1740 GO TO 1200 1750 IF I9<>0 THEN CLOSE #I9 1760 PRINT "";Y;"YEARS" 1310 PRINT #I9," REGULAR PAYMENT = $";R 1320 PRINT #I9, 1330 PRINT #I9,"NO ";"INTEREST","AMORTIZED","BALANCE",ANYWHERE ELSE "; \ INPUT Y$ 1770 IF SEG$(Y$,1,1)="Y" THEN GO TO 1090 1780 PRINT 1790 RETURN "ACCUM INTEREST" 1340 REM - COUNT LINES PRINTED ON EACH PAGE IN K 1350 K=7 1360 K=K+N+3 1370 FOR J1=1 TO N1 1380 REM - CALCULATE INTEREST PAID THIS PAYMENT, ROUND OFF 1390 I1=INT((B0*I/N)*100+.5)/100 1400 C1=C1+1 \ REM - COUNT NUMBER OF PAYMENTS MAD  1440 X(1)=SQR(1-I^2) 1450 X(1)=ATN(X(1)/I) 1460 GO TO 1210 1470 PRINT 1480 REM - PRINT RESULTS 1490 FOR I=1 TO 3 1500 RERIGIN 1320 PRINT "R =0, A = 0" 1330 GO TO 1040 1340 REM - POINT IS ON X-AXIS 1350 PRINT "R =";ABS(X);","; 1360 REM - IS POI1000 PRINT \ A=SYS(7,1) \ REM - "PARTS OF A TRIANGLE" 1010 REM - ENTER NUMBER OF PROBLEM TYPE ACCORDING TO KNOWN PARTS 1020 M - THE ANGLE OF A TRIANGLE CANNOT BE LESS THAN ZERO 1510 IF X(I)<0 THEN 1560 1520 PRINT "SIDE";I;"=";INT(Y(I)*1000+.5)/1000 REM - OF THE TRIANGLE WHERE A=ANGLE, S=LENGTH OF SIDE 1030 PRINT "PROBLEM TYPES: 1=ASA,2=SAS,3=AAS,4=SSA,5=SSS,6=END PROGRAM" 1530 PRINT "OPPOSITE ANGLE=";INT(X(I)*1000+.5)/1000;"RADIANS" 1540 NEXT I 1550 PRINT \ GO TO 1580 1560 PRINT 1570 PRINT \1040 PRINT "ENTER PROBLEM TYPE"; \ INPUT N 1050 REM - DIRECT PROGRAM TO PROPER CALCULATIONS 1060 IF N=6 THEN 1580 1070 IF N=5 PRINT "NO SOLUTION" 1580 RETURN  THEN 1360 1080 IF N=4 THEN 1290 1090 IF N=3 THEN 1250 1100 IF N=2 THEN 1180 1110 PRINT "ENTER ANGLE,SIDE,ANGLE"; 1120 INPU1000 PRINT \ A=SYS(7,1) \ REM - "COORDINATE CONVERSION" 1010 PRINT " ( 1=CARTESIAN TO POLAR)" 1020 PRINT " (-1=POLAR TO CAT X(1),Y(3),X(2) 1130 X(3)=PI-X(1)-X(2) 1140 I=Y(3)/SIN(X(3)) 1150 Y(1)=SIN(X(1))*I 1160 Y(2)=SIN(X(2))*I 1170 GO TO 1470 RTESIAN)" 1030 PRINT " ( 0=END PROGRAM)" 1040 PRINT "WHICH DIRECTION"; 1050 INPUT D 1060 REM - END PROGRAM? 1070 IF D=0 TH1180 PRINT "ENTER SIDE,ANGLE,SIDE"; 1190 INPUT Y(3),X(1),Y(2) 1200 Y(1)=SQR(Y(3)^2+Y(2)^2-2*Y(3)*Y(2)*COS(X(1))) 1210 X(2)=SIEN 1520 1080 REM - DIRECT PROGRAM TO PERFORM PROPER CONVERSION 1090 IF D=-1 THEN 1440 1100 REM - CONVERT FROM CARTESIAN COORDN(X(1))*Y(2)/Y(1) 1214 IF X(2)<1 THEN 1220 1216 X(2)=1.5708 \ GO TO 1230 1220 X(2)=ATN(X(2)/SQR(1-X(2)^2)) 1230 X(3)=PI-X(1)INATES TO POLAR COORDINATES 1110 REM - ENTER CARTESIAN COORDINATES (ABSCISSA,ORDINATE) 1120 PRINT "X,Y"; 1130 INPUT X,Y 1140-X(2) 1240 GO TO 1470 1250 PRINT "ENTER ANGLE,ANGLE,SIDE"; 1260 INPUT X(3),X(2),Y(3) 1270 X(1)=PI-X(2)-X(3) 1280 GO TO 1140 REM - POINT ON Y-AXIS? 1150 IF X=0 THEN 1230 1160 REM - POINT ON X-AXIS? 1170 IF Y=0 THEN 1350 1180 REM - COMPUTE POLAR COO 1290 PRINT "ENTER SIDE,SIDE,ANGLE"; 1300 INPUT Y(1),Y(2),X(1) 1310 I=(SIN(X(1))*Y(2))/Y(1) 1314 IF I<1 THEN 1320 1316 X(2)RDINATES, ROUND OFF, PRINT 1190 PRINT "R =";INT(SQR(X^2+Y^2)*100+.5)/100;","; 1200 A=INT(ATN(Y/X)*180/3.14159*100+.5)/100 120=1.5708 \ GO TO 1330 1320 X(2)=ATN(I/SQR(1-I^2)) 1330 X(3)=PI-X(2)-X(1) 1340 Y(3)=Y(1)*SIN(X(3))/SIN(X(1)) 1350 GO TO 1470 2 IF X<0 THEN 1208 1204 IF Y>0 THEN 1210 1206 A=A+360 \ GO TO 1210 1208 A=A+180 1210 PRINT " A =";A \ GO TO 1040 1220 REM -1360 PRINT "ENTER SIDE,SIDE,SIDE"; 1370 INPUT Y(1),Y(2),Y(3) 1380 I=Y(1) \ J=Y(2) 1390 IF IY(3 POINT IS ON Y-AXIS; AT ORIGIN? 1230 IF Y=0 THEN 1320 1240 PRINT "R =";ABS(Y);","; 1250 REM - IS POINT ABOVE OR BELOW ORIGIN?) THEN Y(1)=Y(3) \ Y(3)=I 1410 GO TO 1430 1420 IF J>Y(3) THEN Y(2)=Y(3) \ Y(3)=J 1430 I=(Y(2)^2+Y(3)^2-Y(1)^2)/(2*Y(2)*Y(3))  1260 IF Y<0 THEN 1290 1270 PRINT " A = 90" 1280 GO TO 1040 1290 PRINT " A = 270" 1300 GO TO 1040 1310 REM - POINT IS AT O OXIMATION FORMULAS 1190 L=ABS((1-K)*Z^(1/3)-1+J)/SQR(K*Z^(2/3)+J) 1200 IF R<4 THEN 1240 1210 X=.5/(1+L*(.196854+L*(.115194+L*20 FOR I=1 TO N 1230 A(2,I-1)=A(1,I)*I 1240 NEXT I 1250 PRINT "GUESS APPROXIMATE ROOT "; \ INPUT X 1260 T=0 1270 S=1 \ D=0 (3.44000E-04+L*.019527))))^4 1220 X=INT(X*10000+.5)/10000 1230 GO TO 1270 1240 L=L*(1+.08*L^4/R^3) 1250 GO TO 1210 1260 REM\ C=0 1280 T=T+1 \ REM COUNT ITERATIONS 1290 FOR I=0 TO N 1300 C=C+A(1,I)*S \ REM CALCULATE VALUE OF FUNCTION 1310 D=D+A(2,I - ADJUST IF INVERSE WAS COMPUTED 1270 IF T>=1 THEN 1290 1280 X=1-X 1290 PRINT "TWO TAIL VALUE = ";X 1300 RETURN )*S \ REM AND DERIVATIVE 1320 S=S*X 1330 NEXT I 1340 IF D=0 THEN PRINT "DERIVATIVE = 0 AT X=";X \ GO TO 1470 1350 S=X-C/D \ REM NEW GUESS 1360 IF ABS(X-S)/(ABS(X)+ABS(S))<5.00000E-06 THEN 1430 \ REM FOUND 1370 X=S 1380 IF T<=100 THEN 1270 1390 PRINT TO LEFT OR RIGHT OF ORIGIN? 1370 IF X<0 THEN 1400 1380 PRINT " A = 0" 1390 GO TO 1040 1400 PRINT " A = 180" 1410 GO TO 1040 1420 REM - CONVERT FROM POLAR COORDINATES TO CARTESIAN COORDINATES 1430 REM - ENTER POLAR COORDINATES (MAGNITUDEOFRAY,ANGL1000 A=SYS(7,1) \ REM ROOTS OF POLYNOMIALS: NEWTON 1010 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" TE) 1440 PRINT "R,A"; 1450 INPUT R,A 1460 REM - CONVERT FROM DEGREES TO RADIANS 1470 M=(A-INT(A/360)*360)*3.14159/180 1480 RHEN I9=0 \ GO TO 1050 1030 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "DEGREEM - CALCULATE CARTESIAN COORDINATES, ROUND OFF, PRINT 1490 PRINT "X =";INT(R*COS(M)*100+.5)/100;","; 1500 PRINT " Y =";INT(R*E OF EQUATION? "; \ INPUT #I9,N \ IF I9<>0 THEN PRINT N 1060 IF N<=A8 THEN GO TO 1090 1070 PRINT "DEGREE TOO LARGE "; \ IF I9SIN(M)*100+.5)/100 1510 GO TO 1040 1520 RETURN <>0 THEN CLOSE #I9 1080 GO TO 1490 1090 FOR I=0 TO N 1100 PRINT "COEFFICIENT";I;"? "; \ INPUT #I9,A(1,I) \ IF I9<>0 THEN PRINT A(1,I) 1110 NEXT I 1120 IF I9<>0 THEN CLOSE #I9 \ GO TO 1220 1130 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA "; \ INPUT Y$1000 A=SYS(7,1) \ REM - "STUDENT'S T-DISTRIBUTION" 1010 PRINT "T-VALUE"; 1020 INPUT T 1030 PRINT "DEGREES OF FREEDOM"; 1040 1140 IF SEG$(Y$,1,1)<>"Y" GO TO 1220 1150 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1160 PRINT "DEGREE OF EQ INPUT D2 1050 D1=1 \ T=T^2 1060 X=1 1070 REM - COMPUTE USING INVERSE FOR SMALL T-VALUES 1080 IF T<1 THEN 1130 1090 S=D1 1UATION ";N; \ INPUT Y$ 1170 IF SEG$(Y$,1,1)="Y" THEN PRINT "DEGREE OF EQUATION = "; \ INPUT N 1180 FOR I=0 TO N 1190 PRINT "C100 R=D2 1110 Z=T 1120 GO TO 1160 1130 S=D2 1140 R=D1 1150 Z=1/T 1160 J=2/9/S 1170 K=2/9/R 1180 REM - COMPUTE USING APPROEFFICIENT";I;"= ";A(1,I); \ INPUT Y$ 1200 IF SEG$(Y$,1,1)="Y" THEN PRINT "COEFFICIENT";I;"= "; \ INPUT A(1,I) 1210 NEXT I 12  77 5 6 7 6 4 9.4 44 81 67 5 4 9 6 10.1 31 93 4 6 7 3 2  5 5 5 45 67 90 7 6 576 5 5 45 467 67 4 5 6 7 8 23 56 78  11.6 29 93 5 8 9 7 5 12.6 58 51 7 7 7 7 7 10.9 37 76 4 3 7  6 6 6 4 0 12 34 56 8 7 6 9 8 12 45 8 6 6 6 7 6 34 78  9 0 23.1 46 96 4 4 4 5 5 23 50 77 8 9 0 9 8 21.6 56 95 7  67 6 5 4 4 4 23 45 67 5 4 6 7 8 12 34 56 4 3 2 5 6  6 9 9 9 23 56 95 7 6 2 2 3 1.9 36 54 55 66 5 4 3 26.8 58 168 8 9 0 3 3 29 51 99 87 0 9 9 9 25 50 98 6 5 4 6 6  23 50 80 7 8 9 9 8 4 54 76 8 9 4 3 6 10 31 93 5 6 89 8 1000 A=SYS(7,1) \ PRINT \ REM - ALPHABETISE 1010 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9 8 9 4 77 6 7 6 5 4 22 44 93 6 9 24 22 1 4 56 78 5 6 7 =0 \ GO TO 1050 1030 PRINT "FILE SPECIFICATION "; \ INPUT F1$ 1040 OPEN F1$ FOR INPUT AS FILE #2 \ I9=2 1050 IF I9=0 THEN PRINT "100 ITERATIONS COMPLETED X =";X;"F(X)=";C 1400 PRINT "CONTINUE"; \ INPUT Y$ 1410 IF SEG$(Y$,1,1)="Y" THEN T=0 \ GO TO 127 89 7 -2 45 60 5 4 4 44 4 34 56 56 0 0 0 0 0 23 50 77 6 0 1420 GO TO 1490 1430 PRINT 1440 PRINT "ROOT","ERROR","DERIVATIVE" 1450 PRINT X,C,D 1460 PRINT 1470 PRINT "ANOTHER ROOT66 78 8 8 12 58 51 4 4 4 4 4 22 45 98 14 23 21 32 1 23 45 "; \ INPUT Y$ 1480 IF SEG$(Y$,1,1)="Y" GO TO 1250 1490 RETURN  76 6 5 4 3 2 38 67 68 9 8 7 6 8 23 56 88 6 6 6 6 6 4  23 45 3 3 3 3 3 23 45 67 6 7 8 8 7 6 23 59 6 5 7 8 5  50 , 8 .4 .4 3.1 .6 4.7 1 9 10 1 12.4 60 1 1 1 2 2 3.1 19 71 3  34 56 78 98 76 5 5 54 23 45 67 6 6 6 6 6 67 5 44 3 2 34  4 2 4 1 .66 34 61 3 2 4 6 2 4.7 24 54 3 1 1 1 2 1.7 65 5 6 45 67 90 7 6 6 5 4 4 56 78 6 6 6 68 8 67 567 78 5 5  HEN GO TO 1220 1160 PRINT "TYPE Y TO CHANGE THE DATA, ANYTHING ELSE IF DATA IS OK " 1170 FOR I=1 TO N 1180 PRINT "ITEM ";I;TALUE OF ALL ENDPOINTS ARE EQUAL 1050 PRINT "ABSOLUTE VALUE OF ENDPOINTS"; \ INPUT Z \ PRINT 1060 REM - CALCULATE INCREMENTS OAB(10);A$(I); \ INPUT Y$ 1190 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW ENTRY "; \ INPUT A$(I) 1200 NEXT I 1210 REM - LOOP TO ALPHAF AXES ACCORDING TO CHARACTERS PER AXIS 1070 PRINT "INCREMENT OF X-AXIS =";Z/30 1080 PRINT "INCREMENT OF Y-AXIS =";Z/12 1090 BETISE DATA 1220 FOR I=1 TO N 1230 FOR J=1 TO N-1 1240 REM - THE REPLACEMENT TECHNIQUE CAUSES 1250 REM - LOWER VALUES TO FPRINT 1100 FOR I=1 TO N 1110 REM - CONVERT DEGREES TO RADIANS 1120 D=.0698132*I 1130 X(I)=INT(((FNC(D)*COS(D)/Z+1)*30)+.5) LOAT OR BUBBLE TO THE TOP 1260 A$=A$(J) 1270 B$=A$(J+1) 1280 IF A$=Y(K) THEN 1230 1220 Y(K+H)=Y(K) \ X(EG$(P1$,1,3)="TT:" THEN I9=0 \ GO TO 1390 1360 IF LEN(P1$)=0 THEN GO TO 1340 1370 OPEN P1$ FOR OUTPUT AS FILE #3 \ I9=3 1380K+H)=X(K) \ K=K-H \ IF K>0 THEN 1210 1230 Y(K+H)=B \ X(K+H)=A 1240 NEXT I 1250 H=INT((H-1)/3) 1260 NEXT J 1270 PRINT "WHERE PRINT #I9,N 1390 FOR I=1 TO N 1400 REM - PRINT ALPHABETISED ITEMS 1410 PRINT #I9,A$(I) 1420 NEXT I 1430 PRINT 1440 IF I9 ARE THE RESULTS TO BE PRINTED - TT: OR LP: "; \ INPUT Y$ 1280 IF SEG$(Y$,1,3)="TT:" THEN I9=0 \ GO TO 1320 1290 REM - NEXT PO<>0 THEN CLOSE #I9 1450 PRINT "ANYWHERE ELSE "; \ INPUT Y$ 1460 IF SEG$(Y$,1,1)="Y" THEN GO TO 1340 1470 RETURN INT TO BE PLOTTED STORED IN T 1300 IF SEG$(Y$,1,3)<>"LP:" THEN GO TO 1270 1310 OPEN Y$ FOR OUTPUT AS FILE #3 \ I9=3 1320 T=1 1330 REM - SKIP POINTS OUT OF Y-POSITIVE RANGE 1340 FOR P=0 TO N-1 1350 IF Y(P+1)>=0 THEN 1380 1360 NEXT P 1370 REM - LOOP TO CALL UP EACH Y-INCREMENT FOR LINES OF PRINT 1380 FOR I=0 TO 24 1390 T=T+P 1400 REM - NUMBER OF POINTS TO BE PLOTTED ON EANT "NUMBER OF ITEMS ? "; 1060 INPUT #I9,N 1070 IF N>0 THEN IF N<=A8 THEN 1100 1080 PRINT "NUMBER OF ITEMS MUST BE IN RANGE 1 - ";A8 \ IF I9=0 THEN 1050 1090 GO TO 1470 1100 FOR I=1 TO N 1110 IF I9=0 THEN PRINT "ITEM";I;" ? "; 1120 INPUT #I9,A$(I) 11000 A=SYS(7,1) \ REM - "PLOT OF POLAR EQUATION" \ PRINT 1010 N=90 \ REM - NUMBER OF POINTS TO BE CALCULATED 1020 PR130 NEXT I \ IF I9<>0 THEN CLOSE #I9 1140 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1150 IF SEG$(Y$,1,1)<>"Y" TINT "FORMULA MUST BE ENTERED AT LINE 20, PROCEED"; \ INPUT Y$ 1030 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 2190 1040 REM - ABSOLUTE V CH LINE STORED IN P 1410 P=0 1420 REM - ALL POINTS PLOTTED? 1430 IF T>N THEN 1470 1440 REM - Y-VALUE ON Y-LINE? 1450 IF Y(T J 2100 REM - LABEL X-AXIS 2110 PRINT #I9,"X"; 2120 PRINT #I9, 2130 NEXT I 2140 REM - LABEL Y-AXIS 2150 PRINT #I9,TAB(30);)=I THEN 1530 1460 REM - PRINT X-AXIS? 1470 IF I=12 THEN 1510 1480 REM - PRINT Y-AXIS 1490 PRINT #I9,TAB(30);"*"; 1500 GO T"Y" 2160 IF I9<>0 THEN CLOSE #I9 2170 PRINT "ANYWHERE ELSE "; \ INPUT Y$ 2180 IF SEG$(Y$,1,1)="Y" THEN GO TO 1270 2190 RETUO 2120 1510 S=N+1 1520 GO TO 1970 1530 FOR L=T TO N 1540 REM - NEXT POINT TO BE PLOTTED ON SAME LINE? 1550 IF Y(L)>Y(T) THERN N 1580 1560 P=P+1 1570 NEXT L 1580 IF P=1 THEN 1710 1590 REM - LOOP TO SORT X-COORDINATES WITH EQUAL Y-COORDINATES; 1600 RE10 DIM F$(2),A$(100),X(100),Y(100) \ X8=100 \ Y8=100 \ A8=100 20 DEF FNC(X)=X 30 A=SYS(7,1) \ REM -ENSURE ALL REPLIES INTERPREM - REORDER SMALLEST TO LARGEST 1610 FOR J=1 TO P 1620 FOR L=1 TO P-J 1630 C=X(T+L-1) 1640 A=X(T+L) 1650 IF C<=A THEN 1680 TED AS UPPER CASE 40 PRINT "DEVICE ON WHICH PACKAGE RESIDES "; \ INPUT D9$ 50 Y$=SEG$(D9$,1,4) 60 IF Y$<>" " THEN IF Y$<>"SY: 1660 X(T+L-1)=A 1670 X(T+L)=C 1680 NEXT L 1690 NEXT J 1700 REM - PRINT X-AXIS? 1710 IF I=12 THEN 1950 1720 L=-1 1730 S=0" THEN IF Y$<>"SY1:" THEN IF Y$<>"SY0:" GO TO 40 70 REM SET UP THE ARRAY OF FUNCTION NAMES 80 OPEN D9$+"MISC.DAT" FOR INPUT AS 1740 FOR K=0 TO P-1 1750 REM - MORE THAN ONE POINT TO BE PLOTTED AT SAME POINT ON GRAPH? 1760 IF X(T+K)=L THEN 1900 1770 L= FILE #1 \ INPUT #1,P9 90 PRINT "FUNCTION CODES " 100 FOR I=1 TO P9 110 INPUT #1,F$(0) 120 PRINT I;" = ";F$(0) 130 NEXT I \X(T+K) 1780 REM - PLOT POINT ON Y-AXIS? 1790 IF L=30 THEN 1850 1800 REM - PLOT POINT TO THE LEFT OF Y-AXIS? 1810 IF L<30 THE RESTORE #1 \ INPUT #1,P9 140 PRINT "FUNCTION CODE "; \ INPUT F9 150 IF F9>0 THEN IF F9<=P9 GO TO 180 160 IF F9=99 THEN GO TN 1870 1820 IF S=1 THEN 1870 1830 REM - PRINT Y-AXIS 1840 PRINT #I9,TAB(30);"*"; 1850 S=1 1860 REM - POINT OUTSIDE OF POSITO 9000 170 PRINT "ILLEGAL FUNCTION CODE" \ GO TO 90 180 FOR I=1 TO F9 190 INPUT #1,F$(0),F$(1) 200 NEXT I \ CLOSE #1 210 PRIVE X-RANGE? 1870 IF L>60 THEN 2120 1880 REM - PLOT POINT 1890 PRINT #I9,TAB(L);"+"; 1900 NEXT K 1910 IF S=1 THEN 2120 192INT "LOADING ";F$(0);"...." 220 OVERLAY D9$+F$(1) 230 GOSUB 1000 240 PRINT "MORE DATA "; \ INPUT Y$ 250 IF SEG$(Y$,1,1)="Y" 0 REM - PRINT Y-AXIS 1930 PRINT #I9,TAB(30);"*"; 1940 GO TO 2120 1950 S=T 1960 REM - LOOP TO PRINT LINE OF X-AXIS 1970 FOR THEN GO TO 230 260 PRINT "DO YOU WANT ANOTHER FUNCTION "; \ INPUT Y$ 270 IF SEG$(Y$,1,1)='Y' THEN GO TO 80 280 GO TO 9000 J=0 TO 60 1980 IF X(S)<>J THEN 2080 1990 REM - PLOT POINT ON X-AXIS 2000 PRINT #I9,"+"; 2010 FOR K=S TO T+P-1 2020 IF X(K)=1000 REM 9000 END X(S) THEN 2050 2030 S=K 2040 GO TO 2090 2050 NEXT K 2060 GO TO 2090 2070 REM - PRINT X-AXIS 2080 PRINT #I9,"*"; 2090 NEXT  1230 REM - LOOP TO INPUT X,Y COORDINATES FOR EACH POINT 1240 PRINT "ARE THE COORDINATES IN A FILE "; \ INPUT Y$ 1250 IF SEG$E 1630 FOR P=0 TO N-1 1640 IF Y(P)>=B1 THEN 1670 1650 NEXT P 1660 REM - LOOP TO CALL UP EACH Y-INCREMENT FOR LINES OF PRINT (Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1280 1260 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1270 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 12 1670 FOR I=B1 TO B2 STEP B3 1680 T=T+P \ P=0 1690 REM - COUNT NUMBER OF POINTS TO BE PLOTTED ON EACH LINE IN P 1700 REM - AL80 FOR I=0 TO N-1 1290 IF I>0 THEN 1310 1300 IF I9=0 THEN PRINT "COORDINATES OF POINT ";I+1; \ GO TO 1320 1310 IF I9=0 THEN PL POINTS PLOTTED? 1710 IF T>N-1 THEN 1750 1720 REM - Y-VALUE ON Y-LINE? IF YES, TEST FOR Y 1730 IF Y(T)=I THEN 1790 1740 REMRINT " POINT ";I+1; 1320 INPUT #I9,X(I),Y(I) 1330 REM - ROUND OFF EACH X,Y TO NEAREST INCREMENT ON AXIS 1340 X(I)=INT((X(I)-A1)/A3+.5)*A3 \ Y(I)=INT((Y(I)-B1)/B3+.5)*B3 1350 IF ABS(X(I)+A1)<=1.00000E-05 THEN X(I)=0 \ GO TO 1370 1360 X(I)=X(I)+A1 1371000 A=SYS(7,1) \ PRINT 1010 REM - INPUT INFORMATION TO SET UP AXES 1020 PRINT "IS WIDTH LONG "; \ INPUT Y$ 1030 IF SEG$(Y$0 IF ABS(Y(I)+B1)<=1.00000E-05 THEN Y(I)=0 \ GO TO 1390 1380 Y(I)=Y(I)+B1 1390 NEXT I \ IF I9<>0 THEN CLOSE #I9 1400 PRINT "W,1,1)="Y" THEN L=131 \ GO TO 1050 1040 L=70 1050 PRINT "X-AXIS: LEFTMOST POINT, RIGHTMOST POINT, INCREMENT"; \ INPUT A1,A2,A3 HERE ARE THE RESULTS TO BE PRINTED - TT: OR LP:"; \ INPUT Y$ 1410 IF SEG$(Y$,1,3)="TT:" THEN I9=0 \ GO TO 1440 1420 IF SEG$(Y$ 1060 IF A3<=0 THEN PRINT "ILLEGAL INCREMENT" \ GO TO 1050 1070 A=(A2-A1)/A3 \ IF A<=L THEN 1090 1080 PRINT "X-RANGE TOO LARGE,1,3)<>"LP:" THEN PRINT "INCORRECT REPLY " \ GO TO 1400 1430 OPEN Y$ FOR OUTPUT AS FILE #3 \ I9=3 1440 REM - CALCULATE ADDITIO" \ GO TO 1050 1090 PRINT "Y-AXIS: LOWEST POINT, HIGHEST POINT, INCREMENT"; \ INPUT B1,B2,B3 1100 IF B3<=0 THEN PRINT "ILLEGALNAL X AND Y COORDINATE 1450 Y(N)=INT(B2+.5)+1 1460 X(N)=INT(A2+.5)+1 1470 PRINT \ REM - NOTE WHERE AXES CROSS 1480 PRINT #I INCREMENT " \ GO TO 1090 1110 B=(B2-B1)/B3 1120 REM - Y-AXIS TOO LONG FOR OUTPUT DEVICE? IF YES, CHANGE ENDPOINTS 1130 REM -9,"INTERSECTION OF AXES AT (";A1;",";B1;")" 1490 PRINT \ REM - SORT COORDINATES REORDER Y(1) TO Y(N)SMALLESTTOLARGEST 1500 H= OR INCREASE INCREMENT 1140 IF B<=20 THEN 1160 1150 PRINT "Y-RANGE TOO LARGE" \ GO TO 1090 1160 PRINT "NUMBER OF POINTS TO BE1 \ A=1 1510 H=3*H+1 \ A=A+1 \ IF H=Y(K) THEN 1580 1570 Y(K+H)=Y(K) \ X(K+H)=X(K) \ K=K-H \ IF K>=0 THEN 1560 1580 Y(K+H)=A \ X(K+H)=BTER NUMBER OF POINTS 1200 I=Y8 \ IF X8-Y8<=0 THEN I=X8 1210 IF NA2 THEN 2300 \ REM10 DIM F$(1),A(15,15),B(15,15) \ A8=15 \ B8=15 20 DEF FNC(X)=X 30 A=SYS(7,1) \ REM -ENSURE ALL REPLIES INTERPRETED AS UPPER CA - TEST FOR OUT OF RANGE X-COORDINATE 2060 IF X(T+J)=Z THEN 2090 \ REM - BYPASS DUPLICATE COORDINATES 2070 PRINT #I9,TAB((Z-A1SE 40 PRINT "DEVICE ON WHICH PACKAGE RESIDES"; \ INPUT D9$ 50 Y$=SEG$(D9$,1,4) 60 IF Y$<>" " THEN IF Y$<>"SY:" THEN IF Y$<>"S)/A3);"+"; \ REM - PLOT POINT 2080 Z=X(T+J) 2090 NEXT J 2100 REM - TEST FOR OUT-OF-RANGE X-COORDINATE 2110 IF Z"SY0:" GO TO 40 70 OPEN D9$+"NUM.DAT" FOR INPUT AS FILE #1 \ INPUT #1,P9 80 PRINT "FUNCTION CODES" 90 FOR I=0 2120 IF Z>A2 THEN 2300 2130 REM - PLOT POINT 2140 PRINT #I9,TAB((Z-A1)/A3);"+"; \ GO TO 2300 2150 S=T+L 2160 REM - LOOP T1 TO P9 100 INPUT #1,F$(0) 110 PRINT I;" = ";F$(0) 120 NEXT I \ RESET #1 \ INPUT #1,P9 130 PRINT "FUNCTION CODE "; \ INPUT FO ESTABLISH PRINT FOR FIRST LINE 2170 FOR J=A1 TO A2 STEP A3 2180 REM - POINT TO BE PLOTTED? 2190 IF X(S)<>J THEN 2270 2200 9 140 IF F9>0 THEN IF F9<=P9 THEN P8=1 \ GO TO 160 150 PRINT "ILLEGAL FUNCTION CODE" \ GO TO 80 160 FOR I=1 TO F9 170 INPUT PRINT #I9,"+"; \ REM - PLOT POINT 2210 REM - BYPASS DUPLICATE COORDINATES 2220 FOR K=S TO T+P-1 2230 IF X(K)=X(S) THEN 2250 #1,F$(0),F$(1) 180 NEXT I \ CLOSE #1 190 PRINT "LOADING ";F$(0);"...." 200 OVERLAY D9$+F$(1) 210 GOSUB 1000 220 PRINT "MORE - FIRST LINE?IF YES, X-AXIS MUST BE PLOTTED 1750 IF I=B1 THEN 1780 1760 REM - PLOT X-AXIS 1770 PRINT #I9,"!"; \ GO TO 2300 2240 S=K \ GO TO 2280 2250 NEXT K 2260 GO TO 2280 2270 PRINT #I9,"-"; \ REM - PLOT X-AXIS 2280 NEXT J 2290 PRINT #I9,"X"; \1780 S=N \ GO TO 2170 1790 FOR L=T TO N-1 1800 REM - NEXT POINT PLOTTED ON SAME LINE? 1810 IF Y(L)>Y(T) THEN 1850 1820 P=P+1 REM - LABEL X-AXIS 2300 PRINT #I9, \ REM - ADVANCE OUTPUT DEVICE TO NEXT LINE 2310 NEXT I 2320 PRINT #I9,"Y" \ REM - LABEL  \ REM - COUNT POINTS TO BE PLOTTED ON EACH LINE 1830 NEXT L 1840 REM - PLOT ONE POINT 1850 IF P=1 THEN 1960 1860 REM - LOOPY-AXIS 2330 IF I9<>0 THEN CLOSE #I9 2340 PRINT "ANYWHERE ELSE "; \ INPUT Y$ 2350 IF SEG$(Y$,1,1)="Y" THEN GO TO 1400 2360 R TO SORT X-COORDINATES WITH EQUAL Y-COORDINATES; 1870 REM - REORDER SMALLEST TO LARGEST 1880 FOR J=1 TO P 1890 FOR L=1 TO P-JETURN  1900 D=X(T+L-1) \ B=X(T+L) 1910 IF D<=B THEN 1940 1920 X(T+L-1)=B \ X(T+L)=D 1930 D=Y(T+L-1) \ Y(T+L-1)=Y(T+L) \ Y(T+L)=D 1940 NEXT L 1950 NEXT J 1960 FOR L=0 TO P-1 1970 Z=X(T+L) 1980 IF Z>=A1 THEN 2000 \ REM - TEST FOR OUT OF RANGE X-COORDINATE 1990 NEXT L 2000 IF I=B1 THEN 2150 \ REM - POINT TO BE PLOTTED ON Y-AXIS? 2010 IF Z=A1 THEN 2030 \ REM - POINT TO BE PLOTTED AS EQUATIONS: ROOTS OF QUADRATIC,QUADEQ.BAS POLYNOMIALS: REAL ROOTS - NEWTON'S METHOD,NEWTON.BAS POLYNOMIALS: HALF-INTERVAL SEARCH,HLFINT.BAS POLYNOMIALS: TRIGONOMETRIC,TRGPOL.BAS EQUATIONS: SIMULTANEOUS,SIMUEQ.BAS MATRIX ADD; SUBTRACT; SCALAR MUL1000 A=SYS(7,1) \ REM - "STUDENT'S T-DISTRIBUTION TEST" 1010 PRINT "IS THE DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>T,MATASM.BAS MATRIX MULTIPLY,MATMUL.BAS MATRIX INVERT,MATINV.BAS END "Y" THEN I9=0 \ GO TO 1050 1030 PRINT "FILENAME "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "TEST 1: MEAN=X" 1060 PRINT "TEST 2: MEAN=MEAN,STANDARD DEVIATION=STANDARD DEVIATION" 1070 PRINT "TEST 3: MEAN=MEAN,STANDARD DEVIATION<>STANDARD DEVIATION" 1080 PRINT "WHICH HYPOTHESIS ? "; 1090 INPUT #I9,T \ IF I9<>0 THEN PRINT T 1100 IF T>=1 THEN IF T<=3 GO TO 1120 1110 IF I9=0 GO TO 1080 \ PRINT "ILLEGAL INPUT" \ CLOSE #I9 \ GO TO 1770 1120 PRINT 1130 REM - INPUT 1 OR 2 SAMPLES DE13 MANN-WHITNEY U TEST,MWUTST.BAS MEAN-VARIANCE-STANDARD DEVIATION,MVSDEV.BAS GEOMETRIC MEAN AND DEVIATION,GEOMD.BAS DISTRIBPENDING ON HYPOTHESIS 1140 FOR I=1 TO SGN(T-1)+1 1150 PRINT "SAMPLE ";I;":" 1160 PRINT " NUMBER OF ELEMENTS (MAXIMUM =";S8;")UTION: BINOMIAL,BINDIS.BAS DISTRIBUTION: POISSON,POIDIS.BAS DISTRIBUTION: NORMAL,NORDIS.BAS DISTRIBUTION: CHI-SQUARE,CHIDIS.B ? "; 1170 INPUT #I9,T(I) \ IF I9<>0 THEN PRINT T(I) 1180 IF T(I)>0 THEN IF T(I)<=S8 GO TO 1200 1190 IF I9=0 GO TO 1160 \ CLOAS DISTRIBUTION: CHI-SQUARE TEST,CHITST.BAS DISTRIBUTION: STUDENT'S T,STTDIS.BAS DISTRIBUTION: STUDENT'S T TEST,STTTST.BAS DSE #I9 \ PRINT "ILLEGAL INPUT" \ GO TO 1770 1200 FOR J=1 TO T(I) 1210 PRINT "ELEMENT ";J;" ? "; 1220 INPUT #I9,S(J,I) \ IF I9 DATA "; \ INPUT Y$ 230 IF SEG$(Y$,1,1)="Y" THEN GO TO 210 240 PRINT "DO YOU WANT ANOTHER FUNCTION "; \ INPUT Y$ 250 IF SEG$ISTRIBUTION: F,FDIS.BAS PERMUTATIONS AND COMBINATIONS,PRMCMB.BAS SYSTEM RELIABILITY,SYSREL.BAS END OF FILE (Y$,1,1)='Y' THEN GO TO 70 260 GO TO 9000 1000 REM 9000 END 7 ALPHABETISE,ALPHAB.BAS PARTS OF A TRIANGLE,PRTTRI.BAS CO-ORDINATE CONVERSION,COCNV.BAS CO-ORDINATE PLOT,COPLT.BAS POLAR E15 POLYGON: AREA OF,POLARE.BAS INTERPOLATION: LINEAR,LINI.BAS INTERPOLATION: CURVILINEAR-LAGRANGE,CURVI.BAS INTEGRATION: SIMQUATION PLOT,PLREQ.BAS DAY OF THE WEEK,DOTWK.BAS DAYS BETWEEN TWO DATES,DBDTS END PSON'S RULE,SIMPSO.BAS INTEGRATION: TRAPEZOIDAL RULE,TRAPEZ.BAS INTEGRATION: GAUSSIAN QUADRATURE,GAUSS.BAS DERIVATIVE,DERIV.B  1260 REM -CHECK THE DATA IF INPUT BY FILE 1270 IF I9<>0 THEN CLOSE #I9 1280 PRINT "DO YOU WISH TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1290 IF SEG$(Y$,1,1)<>"Y" GO TO 1430 1300 PRINT "TYPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK " 1310 PRINT "HYPOTHESIS = ";T; 1320 FOR I=1 TO SGN(T-1)+1 1330 PRINT "SAMPLE ";I;", SIZE = ";T(I); \ INPUT Y$ 1340 IF SEG$(Y$,1,1)="Y" THEN 1000 A=SYS(7,1) \ REM LINEAR INTERPOLATION 1010 PRINT "IS THE DATA IN A FILE "; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9=0PRINT "NEW VALUE = "; \ INPUT T(I) 1350 FOR J=1 TO T(I) 1360 PRINT "ELEMENT ";J;" = ";S(J,I); \ INPUT Y$ 1370 IF SEG$(Y$,1,1) \ GO TO 1050 1030 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "ENTER TWO POI="Y" THEN PRINT "NEW VALUE ="; \ INPUT S(J,I) 1380 NEXT J 1390 NEXT I 1400 IF T<>1 GO TO 1430 1410 PRINT "MEAN =";M; \ INPUTNTS ON THE LINE" 1060 FOR I=1 TO 2 1070 PRINT "X,Y OF POINT";I; 1080 INPUT #I9,A(0,I),A(1,I) \ IF I9<>0 THEN PRINT A(0,I);A(1 Y$ 1420 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW VALUE ="; \ INPUT M 1430 REM 1440 FOR I=1 TO SGN(T-1)+1 1450 C(1,I)=0 1460 C(2,I) 1090 NEXT I 1100 PRINT "ENTER X COORDINATES TO BE INTERPOLATED (99999 TO FINISH)" 1110 FOR I=1 TO A8 1120 PRINT "INTERP,I)=0 1470 FOR J=1 TO T(I) 1480 REM - ACCUMULATE SAMPLES 1490 C(1,I)=C(1,I)+S(J,I) 1500 C(2,I)=C(2,I)+S(J,I)^2 1510 NEXT J OLATION POINT";I;" X = "; 1130 INPUT #I9,A(2,I) \ IF I9<>0 THEN PRINT A(2,I) 1140 IF A(2,I)=99999 THEN L=I-1 \ GO TO 1180 115 1520 REM - COMPUTE INTERMEDIATE VALUES 1530 A(I)=C(1,I)/T(I) 1540 C(1,I)=(C(2,I)-C(1,I)^2/T(I))/(T(I)-1) 1550 NEXT I 1560 P0 NEXT I 1160 L=A8 1170 PRINT "ARRAY FULL - WILL INTERPOLATE THE POINTS ENTERED" 1180 IF I9<>0 THEN CLOSE #I9 \ GO TO 1310 1RINT 1570 IF T=2 THEN 1640 1580 IF T=3 THEN 1690 1590 REM - COMPUTE T AND DEGREES OF FREEDOM FOR FIRST HYPOTHESIS 1600 A=(A190 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1200 IF SEG$(Y$,1,1)<>"Y" GO TO 1310 1210 PRINT "TYPE 'Y' TO CHAN(1)-M)*SQR(T(1)/C(1,1)) 1610 B=T(1)-1 1620 GO TO 1730 1630 REM - COMPUTE T AND DEGREES OF FREEDOM FOR SECOND HYPOTHESIS 1640GE DATA, ANYTHING ELSE IF IT'S OK" 1220 FOR I=1 TO 2 1230 PRINT "POINT";I;"= ";A(0,I),A(1,I); \ INPUT Y$ 1240 IF SEG$(Y$,1,1) A=(A(1)-A(2))/SQR(1/T(1)+1/T(2)) 1650 B=T(1)+T(2)-2 1660 A=A/SQR(((T(1)-1)*C(1,1)+(T(2)-1)*C(1,2))/B) 1670 GO TO 1730 1680 <>"Y" GO TO 1260 1250 PRINT "NEW POINT";I;"= "; \ INPUT A(0,I),A(1,I) 1260 NEXT I 1270 FOR I=1 TO L 1280 PRINT "INTERPOLATE REM - COMPUTE T AND DEGREES OF FREEDOM FOR THIRD HYPOTHESIS 1690 A=(A(1)-A(2))/SQR(C(1,1)/T(1)+C(1,2)/T(2)) 1740 PRINT 1750 X =";A(2,I); \ INPUT Y$ 1290 IF SEG$(Y$,1,1)="Y" THEN PRINT "NEW X = "; \ INPUT A(2,I) 1300 NEXT I 1310 FOR I=1 TO L 1320 Y=<>0 THEN PRINT S(J,I) 1230 NEXT J 1240 NEXT I 1250 IF T=1 THEN PRINT "VALUE OF MEAN "; \ INPUT #I9,M \ IF I9<>0 THEN PRINT M PRINT "T-VALUE =";ABS(A) 1760 PRINT "DEGREES OF FREEDOM: SAMPLE 1 =";T(1)-1;" , SAMPLE 2 =";T(2)-1 1770 RETURN  J) 1280 NEXT J 1290 NEXT I 1300 IF I9<>0 THEN CLOSE #I9 1310 C=A \ D=B 1320 NEXT L 1330 IF N=3 THEN PRINT "SCALAR"; \ INPUT S 1340 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1350 IF SEG$(Y$,1,1)<>"Y" GO TO 1560 1360 PRINT "TYPE 'Y' TA(1,1)+(A(1,2)-A(1,1))/(A(0,2)-A(0,1))*(A(2,I)-A(0,1)) 1330 A(3,I)=INT(Y*1000+.5)/1000 1340 NEXT I 1350 PRINT "X-VALUE";TAB(1O CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1370 PRINT "BOUNDS";A","B; \ INPUT Y$ 1380 IF SEG$(Y$,1,1)="Y" GO TO 1060 1390 FOR L0);"Y-VALUE" 1360 PRINT 1370 FOR I=1 TO L 1380 PRINT A(2,I);TAB(10);A(3,I) 1390 NEXT I 1400 PRINT 1410 RETURN =1 TO K 1400 PRINT \ PRINT "MATRIX ";L 1410 FOR I=1 TO A 1420 PRINT "ROW";I 1430 FOR J=1 TO B 1440 PRINT "COLUMN";J;"="; 1000 A=SYS(7,1) \ REM MATRIX ADDITION, SUBTRACTION AND SCALAR MULTIPLICATION 1010 PRINT "1-ADDITION, 2-SUBTRACTION, 3-SCALAR MU1450 IF L=1 THEN PRINT A(I,J); \ GO TO 1470 1460 PRINT B(I,J); 1470 INPUT Y$ \ IF SEG$(Y$,1,1)<>"Y" GO TO 1510 1480 PRINT "NELTIPLICATION" 1020 PRINT "WHICH OPERATION"; \ INPUT N 1030 IF N>0 THEN IF N<4 THEN GO TO 1050 1040 PRINT "INVALID OPERATION"W VALUE"; 1490 IF L=1 THEN INPUT A(I,J) \ GO TO 1510 1500 INPUT B(I,J) 1510 NEXT J 1520 NEXT I 1530 NEXT L 1540 IF N=3 THE \ GO TO 1020 1050 K=2 \ IF N=3 THEN K=1 1060 FOR L=1 TO K 1070 PRINT "IS MATRIX";L;"DATA IN A FILE"; \ INPUT Y$ 1080 IF SEGN PRINT "SCALAR";S; \ INPUT Y$ 1550 IF SEG$(Y$,1,1)="Y" THEN PRINT "SCALAR"; \ INPUT S 1560 FOR I=1 TO A 1570 FOR J=1 TO B 1$(Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1110 1090 PRINT "FILE SPECIFICATION"; \ INPUT Y$ 1100 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 11580 ON N GO TO 1590,1610,1630 1590 A(I,J)=A(I,J)+B(I,J) 1600 GO TO 1640 1610 A(I,J)=A(I,J)-B(I,J) 1620 GO TO 1640 1630 A(I,10 PRINT "MATRIX";L;"BOUNDS? "; \ INPUT #I9,A,B 1120 IF I9<>0 THEN PRINT A,B 1130 I=A8 \ IF L=2 THEN I=B8 1140 IF A>I THEN PRJ)=A(I,J)*S 1640 NEXT J 1650 NEXT I 1660 PRINT "WHERE ARE THE RESULTS TO BE PRINTED/SAVED - TT:, LP:, OR FILESPEC"; 1670 INPINT "MATRIX ";L;"ROW BOUND CANNOT EXCEED ";I \ GO TO 1800 1150 IF B>I THEN PRINT "MATRIX ";L;"COLUMN BOUND CANNOT EXCEED ";I \ UT Y$ \ L=0 \ IF SEG$(Y$,1,3)="TT:" THEN I9=0 \ GO TO 1710 1680 IF LEN(Y$)=0 THEN I9=0 \ GO TO 1710 1690 I9=3 \ OPEN Y$ FOR OUGO TO 1800 1160 IF L=1 GO TO 1200 1170 IF A=C THEN IF B=D GO TO 1200 1180 PRINT "MATRIX BOUNDS MUST BE THE SAME" \ IF I9=0 GOTPUT AS FILE #3 1700 IF SEG$(Y$,1,3)<>"LP:" THEN L=1 \ PRINT #I9,A;",";B 1710 FOR I=1 TO A 1720 FOR J=1 TO B 1730 PRINT #I9, TO 1110 1190 GO TO 1800 1200 FOR I=1 TO A 1210 PRINT "ROW";I 1220 FOR J=1 TO B 1230 PRINT "COLUMN";J;"? "; 1240 IF L<>1 GA(I,J); \ IF L=1 THEN PRINT #I9 1740 NEXT J 1750 PRINT #I9 1760 NEXT I 1770 IF I9<>0 THEN CLOSE #I9 1780 PRINT "ANYWHERE ELO TO 1270 1250 INPUT #I9,A(I,J) \ IF I9<>0 THEN PRINT A(I,J) 1260 GO TO 1280 1270 INPUT #I9,B(I,J) \ IF I9<>0 THEN PRINT B(I,SE "; \ INPUT Y$ 1790 IF SEG$(Y$,1,1)="Y" THEN GO TO 1660 1800 RETURN !  1130 NEXT J% 1140 K%=K%+1 \ B(K%,0)=Y% 1150 FOR I%=1 TO N 1160 B(I%,1)=1 1170 FOR J%=1 TO M 1180 FOR L%=2 TO K% 1190 IF  1440 NEXT K 1450 NEXT J 1460 NEXT I 1470 PRINT #P1, \ PRINT #P1, \ PRINT #P1, 1480 P=0 1490 FOR I=2 TO V+1 1500 P=P+A(I,J%<>B(L%,0) THEN GO TO 1210 1200 B(I%,L%)=A(I%,J%) \ GO TO 1220 1210 NEXT L% 1220 NEXT J% 1230 NEXT I% 1240 V=K%-2 \ V1=V V+2)*(B(0,I)-A(0,I)*B(0,1)/N) 1510 NEXT I 1520 R=B(0,V+2)-B(0,1)^2/N 1530 Z=R-P 1540 L=N-V-1 1550 S1=Z/L 1560 R2=P/R 1570 1250 RETURN  F=(P/V)/S1 1580 PRINT #P1,"MULTIPLE R-SQUARED ";R2 1590 PRINT #P1,"STANDARD ERROR OF ESTIMATE ";SQR(ABS(S1)) 1600 PRINT #P1, \ PRINT #P1, \ PRINT #P1, 1610 REM FIND INVERSE OF X'X MATRIX 1620 FOR J%=1 TO V+1 1630 FOR I%=J% TO V+1 1640 IF A(I%+V+1,J1000 A=SYS(7,1) \ REM - "COMPUTE REG. COEFF.,ANOVA,RESIDUALS" 1010 REM - POPULATE A MATRIX TO BE USED IN CURVE FITTING 1020 F%)<>0 THEN GO TO 1670 1650 NEXT I% 1660 GO TO 2250 1670 T=A(J%+V+1,J%) 1680 IF T=0 THEN GO TO 1740 1690 FOR K%=1 TO V+1  13 , 5 7 26 36 60 78.5 1 29 225 52 74.3 11 56 64 20 104.3 11 31 64 OR I%=0 TO A8 1030 FOR J%=0 TO B8 \ A(I%,J%)=0 \ NEXT J% 1040 NEXT I% 1050 FOR I=1 TO N 1060 FOR K=1 TO V+1 1070 FOR L=1 TO 47 87.6 7 52 36 33 95.9 11 55 81 22 109.2 3 71 289 6 102.7 1 31 V+2 1080 A(K,L)=A(K,L)+B(I,K)*B(I,L) 1090 B(0,K)=A(K,V+2) 1100 NEXT L 1110 NEXT K 1120 B(0,V+2)=B(0,V+2)+B(I,V+2)^2 1130  484 44 72.5 2 54 324 22 93.1 20 47 16 26 115.9 1 40 529 34 83.8 NEXT I 1140 REM SAVE X'X AND CREATE I MATRIX 1150 FOR I%=1 TO V+1 1160 FOR J%=1 TO V+1 \ A(I%+V+1,J%)=A(I%,J%) 1170 A(I%+2*V 11 66 81 12 113.3 10 68 64 12 109.4 +2,J%)=0 \ IF I%=J% THEN A(I%+2*V+2,J%)=1 1180 NEXT J% 1190 NEXT I% 1200 REM - FIT CURVE BY SOLVING THE SYSTEM OF LINEAR EQUA1000 A=SYS(7,1) \ REM - "MULTIPLE LINEAR REGRESSION" 1010 OVERLAY D9$+"RFILE.BAS" 1020 GOSUB 1500 1030 PRINT "TYPE THE NUMBETIONS IN MATRIX A() 1210 FOR I=2 TO V+1 1220 A(0,I)=A(1,I) 1230 NEXT I 1240 FOR I=1 TO V+1 1250 FOR J=I TO V+1 1260 IF A(JR OF THE DEPENDENT VARIABLE "; \ INPUT Y% 1040 IF Y%>0 THEN IF Y%<=M THEN GO TO 1060 1050 PRINT "NO SUCH VARIABLE." \ GO TO 1,I)<>0 THEN 1300 1270 NEXT J 1280 PRINT #P1, \ PRINT #P1, \ PRINT #P1,"NO UNIQUE REGRESSION MODEL" 1290 GO TO 2250 1300 FOR 030 1060 PRINT "TYPE Y IF YOUR MODEL INCLUDES THE VARIABLE, IF NOT TYPE ANYTHING." 1070 K%=1 1080 FOR J%=1 TO M 1090 IF J%=YK=1 TO V+2 1310 B=A(I,K) 1320 A(I,K)=A(J,K) 1330 A(J,K)=B 1340 NEXT K 1350 Z=1/A(I,I) 1360 FOR K=1 TO V+2 1370 A(I,K)=Z*A% THEN GO TO 1130 1100 PRINT "VARIABLE ";J%; \ INPUT Y$ 1110 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1130 1120 K%=K%+1 \ B(K%,0)=J%(I,K) 1380 NEXT K 1390 FOR J=1 TO V+1 1400 IF J=I THEN 1450 1410 Z=-A(J,I) 1420 FOR K=1 TO V+2 1430 A(J,K)=A(J,K)+Z*A(I,K)" T J% 1920 NEXT I% 1930 PRINT #P1,"TABLE OF REGRESSION COEFFICIENTS" \ PRINT #P1, 1940 PRINT #P1,TAB(13);"REGRESSION";TAB(27);1000 A=SYS(7,1) \ REM MATRIX INVERSION 1010 PRINT "IS MATRIX DATA IN A FILE"; \ INPUT Y$ 1020 IF SEG$(Y$,1,1)<>"Y" THEN I9=0 \"STANDARD" 1950 PRINT #P1,"VAR NO.";TAB(13);"COEFFICIENT";TAB(27);"ERROR";TAB(39);"T-VALUE" 1960 PRINT #P1, \ FOR I%=2 TO V+1  GO TO 1050 1030 PRINT "FILE SPECIFICATION"; \ INPUT Y$ 1040 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1050 PRINT "MATRIX BOUND? "; 1970 IF B(V+1,0)=0 THEN PRINT #P1,"X**";I%-1; \ GO TO 1990 1980 PRINT #P1,B(I%,0); 1990 PRINT #P1,TAB(13);A(I%,V+2);TAB(26);A \ INPUT #I9,R 1060 IF I9<>0 THEN PRINT R 1070 FOR I=1 TO R 1080 PRINT "ROW";I 1090 FOR J=1 TO R 1100 PRINT "COLUMN";J;"? "(I%,V+1);TAB(39);A(I%,V) 2000 NEXT I% 2010 PRINT #P1, \ PRINT #P1, \ PRINT #P1,"CONSTANT TERM ";A(1,V+2) 2020 PRINT #P1, \ PR; 1110 INPUT #I9,A(I,J) \ IF I9<>0 THEN PRINT A(I,J) 1120 B(I,J)=0 1130 NEXT J 1140 B(I,I)=1 1150 NEXT I 1160 IF I9<>0 THEINT #P1, \ PRINT #P1, 2030 PRINT #P1,"ANALYSIS OF VARIANCE" \ PRINT #P1, 2040 PRINT #P1,"SOURCE";TAB(14);"DF";TAB(25); 2050 N CLOSE #I9 1170 PRINT "DO YOU WANT TO CHECK/CHANGE THE DATA"; \ INPUT Y$ 1180 IF SEG$(Y$,1,1)<>"Y" GO TO 1310 1190 PRINT "TYPRINT #P1,"SUM OF SQ.";TAB(39);"MEAN SQ.";TAB(52);"F-VALUE" 2060 PRINT #P1, 2070 PRINT #P1,"REGRESSION";TAB(13);V;TAB(24);P;TAPE 'Y' TO CHANGE DATA, ANYTHING ELSE IF IT'S OK" 1200 PRINT "BOUND";R; \ INPUT Y$ 1210 IF SEG$(Y$,1,1)="Y" GO TO 1050 1220 FOB(38);P/V;TAB(52);F 2080 PRINT #P1,"RESIDUAL";TAB(13);L;TAB(24);Z;TAB(38);S1 2090 PRINT #P1,"TOTAL";TAB(13);N-1;TAB(24);R;TAB(R I=1 TO R 1230 PRINT "ROW";I 1240 FOR J=1 TO R 1250 PRINT "COLUMN";J;"="; 1260 PRINT A(I,J); \ INPUT Y$ \ IF SEG$(Y$,1,1)<>38);R/(N-1) 2100 PRINT #P1, \ PRINT #P1, \ PRINT #P1, 2110 PRINT #P1,"TABLE OF RESIDUALS" \ PRINT #P1, 2120 PRINT #P1,"OBS NO"Y" GO TO 1290 1270 PRINT "NEW VALUE"; 1280 INPUT A(I,J) 1290 NEXT J 1300 NEXT I 1310 FOR J=1 TO R 1320 FOR I=J TO R 13301700 A(J%+V+1,K%)=A(J%+V+1,K%)/T 1710 A(J%+2*V+2,K%)=A(J%+2*V+2,K%)/T 1720 NEXT K% 1730 GO TO 1780 1740 FOR K%=1 TO V+1 175.";TAB(13);"Y";TAB(26);"PRED.Y";TAB(39);"RESIDUALS"; 2130 PRINT #P1,TAB(52);"STD.RES." 2140 PRINT #P1, 2150 FOR I%=1 TO N 210 A(J%+V+1,K%)=A(J%+V+1,K%)+A(I%+V+1,K%) \ A(J%+2*V+2,K%)=A(J%+2*V+2,K%)+A(I%+2*V+2,K%) 1760 NEXT K% 1770 GO TO 1670 1780 FOR60 A(I%,0)=0 2170 FOR J%=2 TO V+1 2180 A(I%,0)=A(I%,0)+A(J%,V+2)*B(I%,J%) \ NEXT J% 2190 A(I%,0)=A(I%,0)+A(1,V+2) 2200 NEXT  L%=1 TO V+1 1790 IF L%=J% THEN 1850 1800 T=-A(L%+V+1,J%) 1810 FOR K%=1 TO V+1 1820 A(L%+V+1,K%)=A(L%+V+1,K%)+T*A(J%+V+1,K%)I% 2210 FOR I%=1 TO N 2220 PRINT #P1,I%;TAB(12);B(I%,V+2);TAB(25);A(I%,0);TAB(38);B(I%,V+2)-A(I%,0); 2230 PRINT #P1,TAB(51);( 1830 A(L%+2*V+2,K%)=A(L%+2*V+2,K%)+T*A(J%+2*V+2,K%) 1840 NEXT K% 1850 NEXT L% 1860 NEXT J% 1870 FOR I%=2 TO V+1 1880 FOR B(I%,V+2)-A(I%,0))/SQR(S1) 2240 NEXT I% \ PRINT #P1, \ PRINT #P1, \ PRINT #P1, 2250 RETURN J%=1 TO V+1 1890 IF I%<>J% THEN GO TO 1910 1900 A(I%,V+1)=SQR(ABS(A(I%+2*V+2,J%)*S1)) \ A(I%,V)=A(I%,V+2)/A(I%,V+1) 1910 NEX# J,K)+A(I,K) \ B(J,K)=B(J,K)+B(I,K) \ NEXT K 1460 GO TO 1370 1470 FOR L=1 TO R 1480 IF L=J THEN 1550 1490 T=-A(L,J) 1500 FORNO SUCH VARIABLE." \ GO TO 1060 1090 PRINT "TYPE THE DEGREE OF YOUR MODEL (UP TO";B8-2;")"; \ INPUT D% 1100 IF D%>0 THEN IF D% K=1 TO R 1510 A(L,K)=A(L,K)+T*A(J,K) 1520 IF INT(A(L,K)*1000+.5)=0 THEN A(L,K)=0 1530 B(L,K)=B(L,K)+T*B(J,K) 1540 NEXT K 1<=B8-2 THEN GO TO 1120 1110 PRINT "ENTER AN INTEGER BETWEEN 1 AND ";B8-2\ GO TO 1090 1120 IF D%<=N-2 THEN GO TO 1140 1130 P550 NEXT L 1560 NEXT J 1570 PRINT "WHERE ARE THE RESULTS TO BE PRINTED/SAVED - TT:, LP:, OR FILESPEC"; 1580 INPUT Y$ \ L=0 \ RINT "DEGREE MUST BE LESS THAN OR EQUAL TO ";N-2\ GO TO 1090 1140 B(2,0)=X% \ B(D%+2,0)=Y% 1150 FOR I%=1 TO N 1160 FOR J%=1 TIF SEG$(Y$,1,3)="TT:" THEN I9=0 \ GO TO 1620 1590 IF LEN(Y$)=0 THEN I9=0 \ GO TO 1620 1600 I9=3 \ OPEN Y$ FOR OUTPUT AS FILE #O D%+1 1170 B(I%,J%)=A(I%,X%)^(J%-1) \ NEXT J% 1180 B(I%,D%+2)=A(I%,Y%) 1190 NEXT I% 1200 V=D% \ V1=1 1210 RETURN 3 1610 IF SEG$(Y$,1,3)<>"LP:" THEN L=1 \ PRINT #I9,R;",";R 1620 FOR I=1 TO R 1630 FOR J=1 TO R 1640 PRINT #I9,INT(B(I,J)*1000+.5)/1000;" "; \ IF L=1 THEN PRINT #I9 1650 NEXT J 1660 PRINT #I9 1670 NEXT I 1680 IF I9<>0 THEN CLOSE #I9 1690 PRINT "ANY1500 A=SYS(7,1) \ REM - EDITING AND CREATING FILE 1510 D1=0 1520 FOR I%=0 TO A8 1530 FOR J%=0 TO B8 \ B(I%,J%)=0 \ A(I%,J%)=WHERE ELSE "; \ INPUT Y$ 1700 IF SEG$(Y$,1,1)="Y" THEN GO TO 1570 1710 RETURN 0 \ NEXT J% 1540 NEXT I% 1550 PRINT "IS THE DATA IN A FILE (Y OR N) "; \ INPUT Y$ 1560 IF SEG$(Y$,1,1)<>"Y" THEN I9=0 \ GO TO 1650 1570 PRINT "FILE SPECIFICATION "; \ INPUT Y$ 1580 OPEN Y$ FOR INPUT AS FILE #2 \ I9=2 1590 INPUT #I9,N,M 1600 FOR I%=1 TO N 1610 FOR J%=1 TO M \ INPUT #I9,A(I%,J%) \ NEXT J% 1620 NEXT I% \ CLOSE #I9 1630 PRINT "YOU HAVE ";M;" VARIABLES." 1640 GO TO 1790 1650 PRINT "SAMPLE SIZE(UP TO";A8;"OBS.)"; \ INPUT N 1660 IF N>1 THEN IF N<=A8 THEN GO TO 1680 1670 PRINT "ENTER1000 A=SYS(7,1) \ REM - "POLYNOMIAL REGRESSION" 1010 OVERLAY D9$+"RFILE.BAS" 1020 GOSUB 1500 1030 PRINT "TYPE THE NUMBER OF  AN INTEGER BETWEEN 2 AND";A8 \ GO TO 1650 1680 IF SEG$(F$(1),1,1)="D" THEN M=1 \ GO TO 1730 1690 PRINT "NUMBER OF VARIABLES(I IF A(I,J)<>0 THEN 1370 1340 NEXT I 1350 PRINT "SINGULAR MATRIX" 1360 GO TO 1710 1370 T=A(J,J) 1380 IF T=0 THEN GO TO 1440THE DEPENDENT VARIABLE "; \ INPUT Y% 1040 IF Y%>0 THEN IF Y%<=M THEN GO TO 1060 1050 PRINT "NO SUCH VARIABLE." \ GO TO 1030  1390 FOR K=1 TO R 1400 A(J,K)=A(J,K)/T 1410 B(J,K)=B(J,K)/T 1420 NEXT K 1430 GO TO 1470 1440 FOR K=1 TO R 1450 A(J,K)=A(1060 PRINT "TYPE THE NUMBER OF THE INDEPENDENT VARIABLE "; \ INPUT X% 1070 IF X%>0 THEN IF X%<=M THEN GO TO 1090 1080 PRINT "$ <>"Y" THEN GO TO 2180 2080 IF I9=0 THEN GO TO 2110 2090 PRINT "TO REPLACE THE OLD DATA WITH THE NEW DATA,TYPE THE SAME FILE NCLUDING DEPENDENT VAR. UP TO";B8-1;")"; 1700 INPUT M 1710 IF M>1 THEN IF M<=B8-1 THEN GO TO 1730 1720 PRINT "ENTER AN INTEGNAME." 2100 PRINT "TYPE ANOTHER FILE NAME IF YOU WANT TO CREATE A FILE WITH THE NEW DATA." 2110 PRINT "FILE SPECIFICATION "; \ER BETWEEN 2 AND";B8-1 \ GO TO 1690 1730 FOR I%=1 TO N 1740 PRINT \ PRINT "OBS. ";I% \ PRINT 1750 FOR J%=1 TO M 1760 PRINT INPUT Y$ 2120 OPEN Y$ FOR OUTPUT AS FILE #3 2130 PRINT #3,N;",";M 2140 FOR I%=1 TO N 2150 FOR J%=1 TO M \ PRINT #3,A(I%,J%) "VARIABLE ";J%; 1770 INPUT A(I%,J%) \ NEXT J% 1780 NEXT I% 1790 PRINT "DO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) "; \ IN; \ PRINT #3 \ NEXT J% 2160 PRINT #3 \ NEXT I% 2170 CLOSE #3 2180 RETURN PUT Y$ 1800 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1920 1810 PRINT \ PRINT \ PRINT 1820 T%=8 \ K%=0 1830 PRINT "OBS NO."; 1840 FOR I%=1 TO M \ K%=K%+T% \ PRINT TAB(K%);"VAR";I%; \ NEXT I% 1850 PRINT \ PRINT 1860 FOR I%=1 TO N \ K%=0 1870 PRINT I%; \ FOR J%=1 TO M \ K%=K%+T% \ PRINT TAB(K%);A(I%,J%); \ NEXT J% 1880 PRINT \ NEXT I% \ PRINT \ PRINT \ PRINT 1890 PRINT "DIS1000 A=SYS(7,1) \ REM - "SUMMARY STATISTICS" 1010 OVERLAY D9$+"RFILE.BAS" 1020 GOSUB 1500 1030 PRINT "TYPE Y IF YOU WANT SUMPLAY THE DATA AGAIN (Y OR N) "; \ INPUT Y$ 1900 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1920 1910 GO TO 1810 1920 IF D1=5 THEN RETURMARY STATISTICS OF THE VARIABLE." 1040 PRINT "IF NOT, TYPE ANYTHING." 1050 K%=1 1060 FOR J%=1 TO M 1070 PRINT "VARIABLE ";J%N 1930 PRINT "DO YOU WANT TO EDIT THE DATA (Y OR N) "; \ INPUT Y$ 1940 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 2050 1950 I8=1 1960 ; \ INPUT Y$ 1080 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1100 1090 K%=K%+1 \ B(K%,0)=J% 1100 NEXT J% 1110 FOR I%=1 TO N 1120 B(I%PRINT "OBS NO.(TYPE 0 IF NO MORE CHANGE) "; \ INPUT I% \ IF I%=0 THEN GO TO 2040 1970 IF I%>0 THEN IF I%<=N THEN GO TO 1990 ,1)=1 \ L%=2 1130 FOR J%=1 TO M 1140 IF J%<>B(L%,0) THEN GO TO 1160 1150 B(I%,L%)=A(I%,J%) \ L%=L%+1 \ GO TO 1160 1160 NEXT1980 PRINT "NO SUCH OBSERVATION NO." \ GO TO 1960 1990 PRINT "VAR NO. "; \ INPUT J% 2000 IF J%>0 THEN IF J%<=M THEN GO TO 202 J% 1170 NEXT I% 1180 V=K%-2 \ V1=V 1190 RETURN 0 2010 PRINT "NO SUCH VARIABLE." \ GO TO 1990 2020 PRINT "NEW VALUE "; \ INPUT A(I%,J%) 2030 GO TO 1960 2040 IF I8=1 THEN GO TO 2060 2050 IF I9<>0 THEN GO TO 2180 2060 PRINT "DO YOU WANT TO SAVE THE DATA (Y OR N) "; \ INPUT Y$ 2070 IF SEG$(Y$,1,1)% 0 IF B(V+1,0)=0 THEN PRINT #P1,"DEGREE ";V 1300 T%=8 \ K%=0 1310 PRINT #P1, \ PRINT #P1, \ PRINT #P1, 1320 PRINT #P1,"INPUT DATA" \ PRINT #P1, 1330 PRINT #P1,"OBS NO."; 1340 FOR I%=1 TO V+1 1350 IF B(I%+1,0)=0 THEN GO TO 1370 1360 K%=K%+T% \ PRINT #P1,TAB(K%);"VAR";B(I%+1,0); 1370 NEXT I% 1380 PRINT #P1, \ PRINT #P1, 1390 FOR I%=1 TO N \ K%=0 1400 PRINT #P1,I%; 1410 F4 STATISTICS: SUMMARY,SUMSTA.BAS REGRESSION: MULTIPLE LINEAR,MULREG.BAS REGRESSION: POLYNOMIAL,POLREG.BAS TRANSFORM: DATA,DAOR J%=2 TO V+2 1420 IF B(J%,0)=0 THEN GO TO 1440 1430 K%=K%+T% \ PRINT #P1,TAB(K%);B(I%,J%); 1440 NEXT J% 1450 PRINT #P1, \TRAN.BAS END 1000 A=SYS(7,1) \ REM -SUMMARY STAT AND CORRELATION MATRIX 1010 FOR I%=0 TO A8 1020 FOR J%=0 TO B8 \ A(I%,J%)=0 \ NEXT J% 103 NEXT I% \ PRINT #P1, \ PRINT #P1, \ PRINT #P1, 1460 PRINT #P1,"SUMMARY STATISTICS" \ PRINT #P1, 1470 PRINT #P1,"VAR NO.";TAB(0 NEXT I% 1040 REM SUMMARY STATISTICS \ FOR J%=2 TO V+2 1050 FOR I%=1 TO N \ A(J%-1,1)=A(J%-1,1)+B(I%,J%) 1060 A(J%-112);"NO.OBS.";TAB(24);"MEAN";TAB(36);"VARIANCE"; 1480 PRINT #P1,TAB(48);"STD.DEV.";TAB(60);"COEF.OF VAR.(%)" 1490 PRINT #P1, \,2)=A(J%-1,2)+B(I%,J%)^2 \ NEXT I% 1070 NEXT J% 1080 FOR I%=1 TO V+1 \ A(I%,2)=(A(I%,2)-A(I%,1)^2/N)/(N-1) \ A(I%,1)=A(I%,1)/N T%=12 \ FOR I%=1 TO V+1 \ K%=11 1500 IF B(I%+1,0)=0 THEN GO TO 1530 1510 PRINT #P1,B(I%+1,0);TAB(11);N; 1520 FOR J%=1 TO 4  1090 A(I%,3)=SQR(A(I%,2)) \ A(I%,4)=A(I%,3)/A(I%,1)*100 \ NEXT I% 1100 REM CORRELATION MATRIX 1110 FOR K%=2 TO V+2 1120 FOR\ K%=K%+T% \ PRINT #P1,TAB(K%);A(I%,J%); \ NEXT J% \ PRINT #P1, 1530 NEXT I% \ PRINT #P1, \ PRINT #P1, \ PRINT #P1, \ PRINT #P1 J%=2 TO V+2 \ D=0 1130 FOR I%=1 TO N \ D=D+B(I%,K%)*B(I%,J%) \ NEXT I% 1140 D=D-A(K%-1,1)*A(J%-1,1)*N 1150 A(K%+V,J%-1)=D/(A,"CORRELATION MATRIX" \ PRINT #P1, 1540 T%=8 \ K%=4 \ PRINT #P1,"VAR NO."; 1550 FOR I%=1 TO V+1 1560 IF B(I%+1,0)=0 THEN 1580(K%-1,3)*A(J%-1,3)*(N-1)) 1160 NEXT J% 1170 NEXT K% 1180 PRINT #P1, \ PRINT #P1, \ PRINT #P1, 1190 P$="MULTIPLE" 1200 IF SE 1570 K%=K%+T% \ PRINT #P1,TAB(K%);B(I%+1,0); 1580 NEXT I% \ PRINT #P1, 1590 PRINT #P1, \ FOR I%=1 TO V+1 \ K%=0 1600 IF B(IG$(F$(1),1,1)="P" THEN P$="POLYNOMIAL" 1210 IF SEG$(F$(1),1,1)="S" THEN P$="SUMMARY" 1220 P1$=" REGRESSION" 1230 IF SEG$(F$(1%+1,0)=0 THEN GO TO 1670 1610 PRINT #P1,B(I%+1,0);" "; 1620 FOR J%=1 TO V+1 1630 IF B(J%+1,0)=0 THEN GO TO 1660 1640 A),1,1)="S" THEN P1$=" STATISTICS" 1240 PRINT #P1,"MAP --- ";P$+P1$ 1250 PRINT #P1,"SAMPLE SIZE ";N 1260 PRINT #P1,"NUMBER (I%+V+1,J%)=INT(A(I%+V+1,J%)*10^4+.5)/10^4 1650 PRINT #P1,USING "###.####",A(I%+V+1,J%); 1660 NEXT J% \ PRINT #P1, 1670 NEXT OF VARIABLES ";V1+1 1270 IF SEG$(F$(1),1,1)="S" THEN GO TO 1300 1280 PRINT #P1,"DEPENDENT VARIABLE IS NOW NO. ";B(V+2,0) 129I% 1680 RETURN & NT "FUNCTION CODE (1-4) "; \ INPUT F9 140 IF F9>0 THEN IF F9<=P9 THEN P8=1 \ GO TO 160 150 PRINT "ILLEGAL FUNCTION CODE" \ GO 0 PRINT "TYPE THE NUMBER OF THE VARIABLE TO BE TRANSFORMED "; \ INPUT Y% 1150 IF Y%>0 THEN IF Y%<=M THEN GO TO 1170 1160 PRINTO 80 160 FOR I=1 TO F9 170 INPUT #1,F$(0),F$(1) 180 NEXT I \ CLOSE #1 190 PRINT "LOADING ";F$(0);"...." 200 OVERLAY D9$+F$T "NO SUCH VARIABLE." \ GO TO 1140 1170 FOR I%=1 TO N 1180 A(I%,Y%)=FNC(A(I%,Y%)) 1190 NEXT I% 1200 PRINT "THE DATA HAS BEEN(1) 210 GOSUB 1000 220 IF SEG$(F$(1),1,1)="D" THEN GO TO 320 230 PRINT "WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:) " TRANSFORMED." \ D1=5 1210 GOSUB 1790 \ REM TO RFILE.BAS TO DISPLAY 1220 GO TO 2060 \ REM TO RFILE.BAS TO SAVE 1230 RETURN ; \ INPUT P1$ 240 P1=0 250 IF SEG$(P1$,1,3)="LP:" THEN OPEN P1$ FOR OUTPUT AS FILE #3 \ P1=3 260 OVERLAY D9$+"REG1.BAS" 270  SUBJ: MAP REGRESSION CHAPTER (IPR #54) CONTENTS CHAPTER 5 REGRESSION ANALYSIS 5-1 5.1 INTRODUCGOSUB 1000 280 IF SEG$(F$(1),1,1)="S" THEN GO TO 310 290 OVERLAY D9$+"REG2.BAS" 300 GOSUB 1000 310 IF P1<>0 THEN CLOSE #P1 TION 5.2 SUMMARY STATISTICS 5.2.1 Method Description 5.2.2 Prompts 5.2.3 Example 5.2.3.1 Input data from keybo 320 PRINT "ANOTHER RUN (Y OR N) "; \ INPUT Y$ 330 IF SEG$(Y$,1,1)="Y" THEN GO TO 200 340 PRINT "DO YOU WANT ANOTHER FUNCTIONard 5.2.3.2 Input data from a file 5.3 MULTIPLE LINEAR REGRESSION 5.3.1 Method Description 5.3.2 Prompts 5.3.3  (Y OR N) "; \ INPUT Y$ 350 IF SEG$(Y$,1,1)='Y' THEN GO TO 70 360 GO TO 9000 1000 REM 9000 END  Example 5.4 POLYNOMIAL REGRESSION 5.4.1 Method Description 5.4.2 Prompts 5.4.3 Example 5.5 DATA TRANSFORMATION 1000 A=SYS(7,1) \ REM - "DATA TRANSFORMATION" 1010 PRINT "TRANSFORMATION EQUATION MUST BE DEFINED BEFORE RUNNING THIS PROGRAM.10 COMMON F$(1),A(50,9),B(50,9),N,V,V1,M \ A8=50 \ B8=9 20 DEF FNC(X)=X 30 A=SYS(7,1) \ REM -ENSURE ALL REPLIES INTERPRETED AS" 1020 PRINT "HAVE YOU DEFINED YOUR EQUATION FOR THIS RUN (Y OR N) "; \ INPUT Y$ 1030 IF SEG$(Y$,1,1)<>"Y" THEN GO TO 1050 1 UPPER CASE 40 PRINT "DEVICE ON WHICH PACKAGE RESIDES (SY0: OR SY1:) "; \ INPUT D9$ 50 Y$=SEG$(D9$,1,4) 60 IF Y$<>" " THEN IF040 GO TO 1120 1050 PRINT "TO DEFINE AN EQUATION, USE THE FOLLOWING FORMAT." 1060 PRINT '20 DEF FNC(X)="USER-DEFINED EQUATION" Y$<>"SY:" THEN IF Y$<>"SY1:" THEN IF Y$<>"SY0:" THEN GO TO 40 70 OPEN D9$+"REG.DAT" FOR INPUT AS FILE #1 \ INPUT #1,P9 80 PRI' 1070 PRINT \ PRINT "FOR EXAMPLE," 1080 PRINT "20 DEF FNC(X)=LOG(2*X)+5" 1090 PRINT \ PRINT "TYPE RUN WHEN YOU ARE DONE." NT "FUNCTION CODES" 90 FOR I=1 TO P9 100 INPUT #1,F$(0) 110 PRINT I;" = ";F$(0) 120 NEXT I \ RESET #1 \ INPUT #1,P9 130 PRI 1100 PRINT \ PRINT "WAIT FOR READY TO ENTER !!!" 1110 PRINT \ GO TO 9000 1120 OVERLAY D9$+"RFILE.BAS" 1130 GOSUB 1500 114' uctions on how they should be answered. For each program, an example is given. 5.2 SUMMARY STATISTICS 5.2.1 Method VARIABLE M ? This prompt displays the variable number. Type the value of the variable. 7. DO YOU WANT T Description This program gives summary statistics and the correlation matrix for the input data. The summary statistics inHE DATA TO BE DISPLAYED (Y OR N) ? Type Y to see the input data. Push the no scroll key to prevent the data displclude the arithmetic mean, variance, standard deviation, and coefficient of variation. Limitations per problem: Total numay from advancing. Push the no scroll key to continue the data display. 8. DISPLAY THE DATA AGAIN (Y OR N) ? ber of variables (M<=8) Total number of observations per variable (N<=50) 5.2.2 Prompts The program issues the followi Type Y to see the data again. 9. DO YOU WANT TO EDIT THE DATA (Y OR N) ? Type Y to edit the data. Otherwise, typeng prompts: 1. IS THE DATA IN A FILE (Y OR N) ? If the data for the program is in a file, type Y. Otherwise,  N. 10. OBS NO. (TYPE 0 IF NO MORE CHANGE) ? This prompt appears only if a reply of Y was given to prompt 9. type N. 2. FILE SPECIFICATION ? This prompt appears only if a reply of Y was given to prompt 1. To input the  Type the observation number of the datum to be changed. Type 0 if no more changes are to be made. 11. VAfile specification, use the following format: dev:name.typ where: dev: SY: for drive 0 or SY1: for drive 1. R NO. ? This prompt appears when prompt 10 was not 0. Type the variable number of the datum to be changed.  name your file name of 1 to 6 alphanumeric characters (0-9,A-Z). .ext standard .ext for a data file is ".DAT"  12. NEW VALUE ? Type the new value. 13. DO YOU WANT TO SAVE THE DATA (Y OR N) ? This prompt appe but any .ext of 1 to 3 alphanumeric characters can be specified. 3. YOU HAVE M VARIABLES. This pars if a data file was created or edited on this run. The prompt will not appear when the replies for prompt 1 and rompt displays the number of the variables in the file specified in prompt 2. 4. SAMPLE SIZE (UP TO 50 OBS.) ? 9 were both N. 14. TO REPLACE THE OLD DATA WITH THE NEW DATA, TYPE THE SAME FILE NAME. TYPE ANOTHER FILE NAME 5.5.1 Method Description 5.5.2 Prompts 5.5.3 Example CHAPTER 5 REGRESSION ANALYSIS  This prompt appears only if you are creating a file (N to prompt 1). Type the number of observations per variable. 5 5.1 INTRODUCTION This chapter describes 4 programs: summary statistics, multiple linear regression, polynomial regress. NUMBER OF VARIABLES (INCLUDING DEPENDENT VAR. UP TO 8) ? This prompt appears when you are creating a file (N to prompion, and data transformation. Each program has a description of the method and a description of the prompts issued with instrt 1). Type the number of variables in the problem. 6. OBS. N This prompt displays the observation number. (  IF YOU WANT TO CREATE A FILE WITH THE NEW DATA. This prompt appears only if the replies to prompts 1 and 9 LE 2 ? 109.1 OBS. 7 VARIABLE 1 ? 17 VARIABLE 2 ? 102.7 OBS. 8 VARIABLE 1 ? 22 VARIABLE 2 ? 72.5 OBS.were both Y. 15. FILE SPECIFICATION ? This prompt appears when a Y was given to prompt 13. Use the for 9 VARIABLE 1 ? 18 VARIABLE 2 ? 93.1 OBS. 10 VARIABLE 1 ? 4 VARIABLE 2 ? 115.9 OBS. 11 VARIABLE 1 ? mat given in prompt 2 to specify the file name to save the data. 16. TYPE Y IF YOU WANT SUMMARY STATISTICS OF23 VARIABLE 2 ? 83.8 OBS. 12 VARIABLE 1 ? 9 VARIABLE 2 ? 113.3 OBS. 13 VARIABLE 1 ? 8 VARIABLE 2 ? 109.4 THE VARIABLE. Answer prompt 17 for each variable. 17. VARIABLE M ? Type Y to display summary st DO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) ? Y OBS NO. VAR 1 VAR 2 1 6 78.5 2 15 74.atistics of the variable. 18. WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:) ? Type TT: to display the re3 3 8 104.3 4 8 87.6 5 6 95.9 6 9 109.1 7 17 102.7 8 22 72.5 9 18 93.1 10 sults on the screen or type LP: for line printer output. 19. ANOTHER RUN (Y OR N) ? Type Y to run the s 4 115.9 11 23 83.8 12 9 113.3 13 8 109.4 DISPLAY THE DATA AGAIN (Y OR N) ? N DO YOU WANT ummary statistics program again. 20. DO YOU WANT ANOTHER FUNCTION (Y OR N) ? Type Y to run another program inTO EDIT THE DATA (Y OR N) ? Y OBS NO. (TYPE 0 IF NO MORE CHANGE) ? 6 VAR NO. ? 2 NEW VALUE ? 109.2 OBS NO. (TYPE 0 IF NO MOR the regression package. 5.2.3 Example 5.2.3.1 Input data from keyboard RUN SY1:REG DEVICE ON WHICH PACKAGE E CHANGE) ? 0 DO YOU WANT TO SAVE THE DATA (Y OR N) ? Y FILE SPECIFICATION ? SY1:SAMPL1.DAT TYPE Y IF YOU WANT SUMMARY STATISRESIDES (SY0: OR SY1:) ? SY1: FUNCTION CODES 1 = STATISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR 3 = REGRESSION: POLYTICS OF THE VARIABLE. IF NOT, TYPE ANYTHING. VARIABLE 1 ? Y VARIABLE 2 ? Y WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:NOMIAL 4 = TRANSFORM: DATA FUNCTION CODE (1-4) ? 1 LOADING STATISTICS: SUMMARY.... IS THE DATA IN A FILE (Y OR N) ? N SAMP) ? TT: MAP --- SUMMARY STATISTICS SAMPLE SIZE 13 NUMBER OF VARIABLES 2 INPUT DATA OBS NO. VAR 1 LE SIZE (UP TO 50 OBS.) ? 13 NUMBER OF VARIABLES (INCLUDING DEPENDENT VAR. UP TO 8) ? 2 OBS. 1 VARIABLE 1 ? 6 VARIABL VAR 2 1 6 78.5 2 15 74.3 3 8 104.3 4 8 87.6 5 6 95.9 6 9 109.2 7 17 102.7 E 2 ? 78.5 OBS. 2 VARIABLE 1 ? 15 VARIABLE 2 ? 74.3 OBS. 3 VARIABLE 1 ? 8 VARIABLE 2 ? 104.3 OBS. 4  8 22 72.5 9 18 93.1 10 4 115.9 11 23 83.8 12 9 113.3 13 8 109.4 SUMMARY STATIST VARIABLE 1 ? 8 VARIABLE 2 ? 87.6 OBS. 5 VARIABLE 1 ? 6 VARIABLE 2 ? 95.9 OBS. 6 VARIABLE 1 ? 9 VARIABICS VAR NO. NO.OBS. MEAN VARIANCE STD.DEV. COEF. OF VAR.(%) 1 13 11.7692 41.0256 6.40513 ) CS VAR NO. NO.OBS. MEAN VARIANCE STD.DEV. COEF.OF VAR.(%) 1 13 11.7692 41.0256 6.40. ANOTHER RUN (Y OR N) ? Type Y to run the multiple linear regression program again. 21. DO YOU WANT ANOTHER FUNC513 54.4226 2 13 95.4231 226.313 15.0437 15.7653 CORRELATION MATRIX VAR NO. 1 TION (Y OR N) ? Type Y to run another program in the regression package. 5.3.3 Example RUN SY1:REG DEVI 2 1 1.0000 -0.5347 2 -0.5347 1.0000 ANOTHER RUN (Y OR N) ? N DO YOU WANT ANOTHER FUNCTION (Y OR N) ? CE ON WHICH PACKAGE RESIDES (SY0: OR SY1:) ? SY1: FUNCTION CODES 1 = STATISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR  54.4226 2 13 95.4231 226.313 15.0437 15.7653 CORRELATION MATRIX VAR NO. 1 2 1 N 5.3 MULTIPLE LINEAR REGRESSION 5.3.1 Method description This program gives summary statistics, the correl 1.0000 -0.5347 2 -0.5347 1.0000 ANOTHER RUN (Y OR N) ? N DO YOU WANT ANOTHER FUNCTION (Y OR N) ? N 5.ation matrix, the multiple correlation coefficient, the standard error of estimate, regression coefficients, the analysis of 2.3.2 Input data from a file RUN SY1:REG DEVICE ON WHICH PACKAGE RESIDES (SY0: OR SY1:) ? SY1: FUNCTION CODES 1 = STATvariance table, and the table of residuals. The multiple regression model is of the following form: Y = b + b X + b X +...ISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR 3 = REGRESSION: POLYNOMIAL 4 = TRANSFORM: DATA FUNCTION CODE (1-4) ? 1 L..+ b X where Y is the dependent variable. X ,X ,...,X are the independent variables. b ,b ,...,b are the regrOADING STATISTICS: SUMMARY.... IS THE DATA IN A FILE (Y OR N) ? Y FILE SPECIFICATION? SY1:SAMPL1.DAT YOU HAVE 2 VARIABLES. Dession coefficients. Limitations per problem: Total number of observations per variable (N<=50) Total number of indepeO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) ? N DO YOU WANT TO EDIT THE DATA (Y OR N) ? N TYPE Y IF YOU WANT SUMMARY STATISTIndent variables (n<=7) 5.3.2 Prompts The program issues the following prompts: 1-15. Refer to summary statistics secCS OF THE VARIABLE. IF NOT, TYPE ANYTHING. VARIABLE 1 ? Y VARIABLE 2 ? Y WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:) tion. 16. TYPE THE NUMBER OF THE DEPENDENT VARIABLE ? Type the number of the dependent variable from the input dat? TT: MAP --- SUMMARY STATISTICS SAMPLE SIZE 13 NUMBER OF VARIABLES 2 INPUT DATA OBS NO. VAR 1 VAR a. 17. TYPE Y IF YOUR MODEL INCLUDES THE VARIABLE, IF NOT TYPE ANYTHING. Answer prompt 18 for each independent var 2 1 6 78.5 2 15 74.3 3 8 104.3 4 8 87.6 5 6 95.9 6 9 109.2 7 17 102.7 iable. 18. VARIABLE M ? Type Y to include the variable in the model. 19. WHERE DO YOU WANT THE RESULTS  8 22 72.5 9 18 93.1 10 4 115.9 11 23 83.8 12 9 113.3 13 8 109.4 SUMMARY STATISTIPRINTED (TT: OR LP:) ? Type TT: to display the results on the screen or type LP: for line printer output. 20* ULTS PRINTED (TT: OR LP:) ? TT: MAP --- MULTIPLE REGRESSION SAMPLE SIZE 13 NUMBER OF VARIABLES 4 DEPENDENT VA18598 4 87.6 89.763 -2.16302 -.910154 5 95.9 95.3749 .525078 .220941 6 109.2 105.423 3.77711 1.58933 RIABLE IS NOW NO. 5 INPUT DATA OBS NO. VAR 1 VAR 3 VAR 4 VAR 5 1 7 6 60 78.5 2 1 15 52 74.3 3 1 7 102.7 104.012 -1.31248 -.552262 8 72.5 75.4322 -2.93224 -1.23383 9 93.1 92.2658 .834175 1 8 20 104.3 4 11 8 47 87.6 5 7 6 33 95.9 6 11 9 22 109.2 7 3 17 6 102.7 8 1 22 44 72.5 9 2 18 22 93.1  .351003 10 115.9 115.42 .479591 .201802 11 83.8 81.4502 2.34984 .988761 12 113.3 111.851 1.44913 .609765 13 1 10 21 4 26 115.9 11 1 23 34 83.8 12 11 9 12 113.3 13 10 8 12 109.4 SUMMARY STATISTICS VAR NO. NO.OB09.4 111.209 -1.80908 -.761223 DO YOU WANT ANOTHER RUN (Y OR N) ? N DO YOU WANT ANOTHER FUNCTION (Y OR N) ? S. MEAN VARIANCE STD.DEV. COEF.OF VAR.(%) 1 13 7.46154 34.6026 5.88239 78.8362 3 13 11.7N 5.4 POLYNOMIAL REGRESSION 5.4.1 Method Description This program gives summary statistics, the correlation692 41.0256 6.40513 54.4226 4 13 30 280.167 16.7382 55.7939 5 13 95.4231 226.313 15.0437 15. matrix, the multiple correlation coefficient, the standard error of estimate, regression coefficients, the analysis of vari3 = REGRESSION: POLYNOMIAL 4 = TRANSFORM: DATA FUNCTION CODE (1-4) ? 2 LOADING REGRESSION: MULTIPLE LINEAR.... IS THE DATA 7653 CORRELATION MATRIX VAR NO. 1 3 4 5 1 1.0000 -0.8241 -0.2454 0.7307 3 -0.8241 IN A FILE (Y OR N) ? Y FILE SPECIFICATION ? SY1:SAMPLE.DAT YOU HAVE 5 VARIABLES. DO YOU WANT THE DATA TO BE DISPLAYED (Y OR N1.0000 0.0295 -0.5347 4 -0.2454 0.0295 1.0000 -0.8213 5 0.7307 -0.5347 -0.8213 1.0000 MULTIPLE R-SQ) ? Y OBS NO. VAR 1 VAR 2 VAR 3 VAR 4 VAR 5 1 7 26 6 60 78.5 2 1 29 15 52 74.3 3 11 56 8 20 104.3 4UARED .981283 STANDARD ERROR OF ESTIMATE 2.37655 TABLE OF REGRESSION COEFFICIENTS REGRESSION STANDARD VAR NO. COE 11 31 8 47 87.6 5 7 52 6 33 95.9 6 11 55 9 22 109.2 7 3 71 17 6 102.7 8 1 31 22 44 72.5 9 2 54 18 22 93.1FFICIENT ERROR T-VALUE 1 1.05184 .223675 4.70255 3 -.410054 .199224 -2.05826 4 -.642797 .0445423 -14. 10 21 47 4 26 115.9 11 1 40 23 34 83.8 12 11 66 9 12 113.3 13 10 68 8 12 109.4 DISPLAY THE DATA AGAIN (Y4312 CONSTANT TERM 111.685 ANALYSIS OF VARIANCE SOURCE DF SUM OF SQ. MEAN SQ. F-VALUE REGRESSION 3 2664.93 888. OR N) ? N DO YOU WANT TO EDIT THE DATA (Y OR N) ? N TYPE THE NUMBER OF THE DEPENDENT VARIABLE ? 5 TYPE Y IF YOUR MODEL INCLU309 157.279 RESIDUAL 9 50.8318 5.64798 TOTAL 12 2715.76 226.313 TABLE OF RESIDUALS OBS NO. Y PRED.Y RESIDUALS DES THE VARIABLE, IF NOT TYPE ANYTHING. VARIABLE 1 ? Y VARIABLE 2 ? VARIABLE 3 ? Y VARIABLE 4 ? Y WHERE DO YOU WANT THE RESSTD.RES. 1 78.5 78.0194 .480598 .202225 2 74.3 73.1602 1.13976 .479584 3 104.3 107.119 -2.81854 -1.+ EPENDENT VARIABLE ? Type the number of the dependent variable from the input data. 17. TYPE THE NUMBER OF THE INDE4 115.9 11 23 83.8 12 9 113.3 13 8 109.4 SUMMARY STATISTICS VAR NO. NO.OBS. MEAN VARIANCE STD.DEVPENDENT VARIABLE ? Type the number of the independent variable from the input data. 18. TYPE THE DEGREE OF YOUR MO. COEF.OF VAR.(%) 3 13 11.7692 41.0256 6.40513 54.4226 5 13 95.4231 226.313 15.0437 15.7653 CORRELATION MDEL (UP TO 7) ? Type the degree of your model. 19. WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:) ? ATRIX VAR NO. 3 5 3 1.0000 -0.5347 5 -0.5347 1.0000 MULTIPLE R-SQUARED .297212 STANDARD ERType TT: to display the results on the screen or type LP: for line printer output. 20. ANOTHER RUN (Y OR N) ? ROR OF ESTIMATE 14.5625 TABLE OF REGRESSION COEFFICIENTS REGRESSION STANDARD VAR NO. COEFFICIENT ERROR T-VALUE  Type Y to run the polynomial regression program again. 21. ANOTHER FUNCTION (Y OR N) ? Type Y to run another pr X**1 -1.81166 12.5967 -.14382 X**2 .108814 .999954 .108819 X**3 -3.88806E-03 .0242545 -.160302 CONSTANT TERM 1ogram in the regression package. 5.4.3 Example RUN SY1:REG DEVICE ON WHICH PACKAGE RESIDES (SY0: OR SY1:) ? SY1: 09.642 ANALYSIS OF VARIANCE SOURCE DF SUM OF SQ. MEAN SQ. F-VALUE REGRESSION 3 807.156 269.052 1.26871 RESIDUALFUNCTION CODES 1 = STATISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR 3 = REGRESSION: POLYNOMIAL 4 = TRANSFORM: DATA F 9 1908.6 212.067 TOTAL 12 2715.76 226.313 TABLE OF RESIDUALS OBS NO. Y PRED.Y RESIDUALS STD.RES. 1 78.5 10UNCTION CODE (1-4) ? 3 LOADING REGRESSION: POLYNOMIAL.... IS THE DATA IN A FILE (Y OR N) ? Y FILE SPECIFICATION ? SY1:SAMPLE.1.85 -23.3496 -1.60341 2 74.3 93.8283 -19.5283 -1.341 3 104.3 100.122 4.17773 .286883 4 8DAT YOU HAVE 5 VARIABLES. DO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) ? N DO YOU WANT TO EDIT THE DATA (Y OR N) ? N TYPE T7.6 100.122 -12.5223 -.859897 5 95.9 101.85 -5.94965 -.408559 6 109.2 99.3167 9.88326 .6786ance, and the table of residuals. The polynomial regression model is of the following form: Y = b + b X + b X +....+ b X HE NUMBER OF THE DEPENDENT VARIABLE ? 5 TYPE THE NUMBER OF THE INDEPENDENT VARIABLE ? 3 TYPE THE DEGREE OF YOUR MODEL (UP TO 7 where Y is the dependent variable X is the independent variable b ,b ,b ,...,b are the regression coefficients Li) ? 3 WHERE DO YOU WANT THE RESULTS PRINTED (TT: OR LP:) ? TT: MAP --- POLYNOMIAL REGRESSION SAMPLE SIZE 13 NUMmitations per problem: Total number of observations per variable (N<=50) Degree of a model (n<=7) 5.4.2 Prompts TBER OF VARIABLES 2 DEPENDENT VARIABLE IS NOW NO. 5 INPUT DATA OBS NO. VAR 1 VAR 2 1 6 78.5 2 he program issues the following prompts: 1-15. Refer to the summary statistics section. 16. TYPE THE NUMBER OF THE D15 74.3 3 8 104.3 4 8 87.6 5 6 95.9 6 9 109.2 7 17 102.7 8 22 72.5 9 18 93.1 10 , 10 115.9 103.888 12.0123 .824881 11 83.8 78.2308 5.56924 .382437 12 113.3 99.3167 13.9833 .960223 13 109.4 100.1rs only if an N was given to prompt 1. Follow the above instructions to define an equation. 3-17. Refer to prompts 122 9.27773 .637097 ANOTHER RUN (Y OR N) ? N ANOTHER FUNCTION (Y OR N) ? N 5.5 DATA TRANSFORMATION 5-15 in the summary statistics section. These prompts appear when a Y was given to prompt 1. 18. THE DATA HAS BEEN .5.1 Method Description This program transforms input data with the user defined equation. The transformed data can be saveTRANSFORMED. The input data has been transformed by the equation defined at the BASIC program line 20. 19. d so that it can be used as input data to any programs in the regression package. The user defined equation can be of any foDO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) ? Type Y to display the transformed data. 20. DISPLAY THE DATA AGAINrm as long as it conforms to MINC BASIC syntax. Here are some examples: FNC(X) = A*X^B + C FNC(X) = A*LOG(B*X) + C FNC (Y OR N) ? Type Y to display the transformed data again. 21. DO YOU WANT TO SAVE THE DATA (Y OR N) ? Typ(X) = A*EXP(B*X) + C FNC(X) = A*SIN(B*X) + C FNC(X) = A*COS(B*X) + C The user defined equation must be entered before rue Y to save the transformed data. 22. TO REPLACE THE OLD DATA WITH THE NEW DATA, TYPE THE SAME FILE NAME. TYPE ANOnning this program. To define an equation, use the following format: 20 DEF FNC(X)="USER DEFINED EQUATION" The statement labeTHER FILE NAME IF YOU WANT TO CREATE A FILE WITH THE NEW DATA. This prompt appears only if you used an old filel must be 20. Do not use another label. For example, 20 DEF FNC(X)=LOG(2*X)+5 Limitations per problem: Total number. 23. FILE SPECIFICATION ? Type the name of the file to be created using the format dev:name.typ. 24.  of observations per variable (N<=50) Number of variables that can be transformed (V=1) 5.5.2 Prompts The program ANOTHER RUN (Y OR N) ? Type Y to run the data transformation program again. 25. ANOTHER FUNCTION (Y OR N) ?  issues the following prompts: 1. TRANSFORMATION EQUATION MUST BE DEFINED BEFORE RUNNING THIS PROGRAM.  Type Y to run another program in the regression package. 5.5.3 Example RUN SY1:REG DEVICE ON WHICH PACKAGE R HAVE YOU DEFINED YOUR EQUATION FOR THIS RUN (Y OR N) ? Type Y if the equation has been defined. Otherwise type N. ESIDES (SY0: OR SY1:) ? SY1: FUNCTION CODES 1 = STATISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR 3 = REGRESSION: POLYN 2. TO DEFINE AN EQUATION,USE THE FOLLOWING FORMAT. 20 DEF FNC(X)="USER DEFINED EQUATION" FOR EXAMPLE, OMIAL 4 = TRANSFORM: DATA FUNCTION CODE (1-4) ? 4 LOADING TRANSFORM: DATA.... TRANSFORMATION EQUATION MUST BE DEFINED BEFOR78 7 102.7 91.1892 11.5108 .790438 8 72.5 81.0517 -8.55172 -.587242 9 93.1 89.6129 3.48705 .239454  20 DEF FNC(X)=LOG(2*X)+5 TYPE RUN WHEN YOU ARE DONE. WAIT FOR READY TO ENTER !!! This prompt appea- AR 1 VAR 2 VAR 3 VAR 4 VAR 5 1 49 26 6 60 78.5 2 1 29 15 52 74.3 3 121 56 8 20 104.3 4 121 31 8 47 87.6 5 49 52 6 33 95.9 6 121 55 9 22 109.2 7 9 71 17 6 102.7 8 1 31 22 44 72.5 9 4 54 18 22 93.1 10 441 47 4 26 115.9 E RUNNING THIS PROGRAM. HAVE YOU DEFINED YOUR EQUATION FOR THIS RUN (Y OR N) ? N TO DEFINE AN EQUATION, USE THE FOLLOWING FORM 11 1 40 23 34 83.8 12 121 66 9 12 113.3 13 100 68 8 12 109.4 DISPLAY THE DATA AGAIN (Y OR N) ? N DO YOU WANT TOAT. 20 DEF FNC(X)="USER DEFINED FUNCTION" FOR EXAMPLE, 20 DEF FNC(X)=LOG(2*X)+5 TYPE RUN WHEN YOU ARE DONE. WAIT FOR  SAVE THE DATA (Y OR N) ? Y TO REPLACE THE OLD DATA WITH THE NEW DATA, TYPE THE SAME FILE NAME. TYPE ANOTHER FILE NAME IF YOU READY TO ENTER !!! READY 20 DEF FNC(X)=X^2 RUN REG 02-JAN-79 04:40:40 DEVICE ON WHICH PACKAGE RESIDES (SY0: OR SYWANT TO CREATE A FILE WITH THE NEW DATA. FILE SPECIFICATION ? SY1:SAMPL3.DAT ANOTHER RUN (Y OR N) ? N ANOTHER FUNCTION (Y OR 1:) ? SY1: FUNCTION CODES 1 = STATISTICS: SUMMARY 2 = REGRESSION: MULTIPLE LINEAR 3 = REGRESSION: POLYNOMIAL 4 = TRANSFN) ? N ORM: DATA FUNCTION CODE (1-4) ? 4 LOADING TRANSFORM: DATA.... TRANSFORMATION FUNCTION MUST BE DEFINED BEFORE RUNNING THIS PROGRAM. HAVE YOU DEFINED YOUR FUNCTION FOR THIS RUN (Y OR N) ? Y IS THE DATA IN A FILE (Y OR N) ? Y FILE SPECIFICATION ? SY1:SAMPLE.DAT YOU HAVE 5 VARIABLES. DO YOU WANT THE DATA DISPLAYED ? Y OBS NO. VAR 1 VAR 2 VAR 3 VAR 4 VAR 5 1 7 26 6 60 78.5 2 1 29 15 52 74.3 3 11 56 8 20 104.3 4 11 31 8 47 87.6 5 7 52 6 33 95.9 6 11 55 9 22 109.2 7 3 7ΊΊΊ͠͠ΊʊD@}.+ʊ|&"4ΊΊP1 17 6 102.7 8 1 31 22 44 72.5 9 2 54 18 22 93.1 10 21 47 4 26 115.9 11 1 40 23 34 83.8 12 11 66 9 12 113.3 13 NONAMEDIR  10 68 8 12 109.4 DISPLAY THE DATA AGAIN (Y OR N) ? N DO YOU WANT TO EDIT THE DATA (Y OR N) ? N TYPE THE NUMBER OF THE VARIABLE TO BE TRANSFORMED ? 1 THE DATA HAS BEEN TRANSFORMED. DO YOU WANT THE DATA TO BE DISPLAYED (Y OR N) ? Y OBS NO. V. / 0 1 2 3 4 5 6 7 8 9 : ; < = > ? @ A B C D E F G H I J K L