Matrix math operations.
Gallery
Source Code
1000 CLS 1010 PRINT "NUMBER OF OBSERVATIONS="; 1020 INPUT NOB 1030 PRINT NOB 1040 PRINT "NO. INDEP. VARIABLES ="; 1050 INPUT NV 1060 PRINT NV 1070 PRINT 1080 DIM X(NOB,NV+1) 1085 DIM Y(NOB,1) 1090 FOR I=1 TO NOB 1095 LET X(I,1)=1 1097 CLS 1100 PRINT ,,"OBSERVATION ",I,,, 1102 PRINT " DEP. VAR. = ", 1103 INPUT Y(I,1) 1104 PRINT Y(I,1) 1105 FOR J=2 TO NV+1 1110 PRINT " VAR ";J;" = ", 1120 INPUT X(I,J) 1130 PRINT X(I,J) 1140 NEXT J 1150 PRINT 1160 PRINT "CORRECT?" 1170 INPUT A$ 1180 IF A$(1)<>"Y" THEN GOTO 1097 1190 NEXT I 1200 REM %T%R%A%N%S%O%S%E% %X 1220 LET ROWA=NV+1 1225 LET COLA=NOB 1230 DIM A(ROWA,COLA) 1240 FOR I=1 TO COLA 1250 FOR J=1 TO ROWA 1260 LET A(J,I)=X(I,J) 1280 NEXT J 1290 NEXT I 1300 LET ROWB=NOB 1310 LET COLB=1 1320 DIM B(ROWB,COLB) 1330 FOR I=1 TO ROWB 1340 LET B(I,1)=Y(I,1) 1350 NEXT I 1360 GOSUB 7300 1370 DIM T(NV+1,1) 1380 FOR I=1 TO NV+1 1390 LET T(I,1)=C(I,1) 1400 NEXT I 1410 LET ROWB=NOB 1420 LET COLB=NV+1 1430 DIM B(ROWB,COLB) 1440 FOR I=1 TO ROWB 1450 FOR J=1 TO COLB 1460 LET B(I,J)=X(I,J) 1470 NEXT J 1480 NEXT I 1490 GOSUB 7300 1500 DIM A(NV+1,NV+1) 1510 LET ROWA=NV+1 1520 LET COLA=NV+1 1530 FOR I=1 TO ROWA 1540 FOR J=1 TO COLA 1550 LET A(I,J)=C(I,J) 1560 NEXT J 1570 NEXT I 1580 LET ROWB=NV+1 1590 LET COLB=1 1600 DIM B(ROWB,COLB) 1610 FOR I=1 TO ROWB 1630 LET B(I,1)=T(I,1) 1640 NEXT I 1650 GOSUB 7500 1660 FOR I=1 TO ROWA 1670 FOR J=1 TO COLA 1680 LET A(I,J)=C(I,J) 1690 NEXT J 1700 NEXT I 1710 GOSUB 7300 1800 FOR I=1 TO NV+1 1810 PRINT I,C(I,1) 1820 NEXT I 1830 FOR I=1 TO NOB 1832 LET P=0 1835 FOR J=1 TO NV+1 1840 LET P=P+C(J,1)*X(I,J) 1845 NEXT J 1850 PRINT Y(I,1),INT (P*100+.5)/100 1860 NEXT I 1900 STOP 7000 REM %T%H%E% %S%U%B%R%O%U%T%I%N%E%S% %I%N% %T%H%E% %7%0%0%0% %B%L%O%C%K% %O%F% %T%H%I%S% %P%R%O%G%R%A%M% %A%R%E% %T%H%E% %M%A%T%R%I%X% %O%P%E%R%A%T%I%O%N%S%.% %7%1%0%0% %=% %S%C%A%L%A%R% %M%U%L%T%I%P%L%I%C%A%T%I%O%N%,% %7%2%0%0% %=% %M%A%T%R%I%X% %A%D%D%I%T%I%O%N% %A%N%D% %S%U%B%T%R%A%C%T%I%O%N%,% %7%3%0%0% %=% %M%A%T%R%I%X% %M%U%L%T%I%P%L%I%C%A%T%I%O%N%,% %7%4%0%0% %=% %M%A%T%R%I%X% %T%R%A%N%S%P%O%S%I%T%I%O%N%,% %7%5%0%0% %=% %M%A%T%R%I%X% %I%N%V%E%R%S%I%O%N%. 7100 REM %S%C%A%L%A%R% %M%U%L%T%I%P%L%I%C%A%T%I%O%N%:%I%N%P%U%T%S%:% %S%C%A%L%A%R%,% %A% %(%M%A%T%R%I%X%)%,% %R%O%W%A%,% %C%O%L%A%.% %O%U%T%P%U%T%:% %C% %(%M%A%T%R%I%X%)% %W%H%E%R%E% %C%=%S%C%A%L%A%R%*%A 7110 DIM C(ROWA,COLA) 7120 FOR I=1 TO ROWA 7130 FOR J=1 TO COLA 7140 LET C(I,J)=SCALAR*A(I,J) 7150 NEXT J 7160 NEXT I 7170 RETURN 7200 REM %M%A%T%R%I%X% %A%D%D%I%T%I%O%N% %A%N%D% %S%U%B%T%R%A%C%T%I%O%N%.% %I%N%P%U%T%S%:% %A% %(%M%A%T%R%I%X%)%,% %B% %(%M%A%T%R%I%X%)%,% %R%O%W%A%,% %C%O%L%A%,% %S%I%G%N% %(%+%1% %F%O%R% %A%D%D%,% %-%1% %F%O%R% %S%U%B%T%.%)%.% %O%U%T%P%U%T%S%:% %C% %(%M%A%T%R%I%X%)% %W%H%E%R%E% %C%=%A%+% %O%R% %-%B%. 7210 DIM C(ROWA,COLA) 7220 FOR I=1 TO ROWA 7230 FOR J=1 TO COLA 7240 LET C(I,J)=A(I,J)+SIGN*B(I,J) 7250 NEXT J 7260 NEXT I 7270 RETURN 7300 REM %M%A%T%R%I%X% %M%U%L%T%I%P%L%I%C%A%T%I%O%N%.% %I%N%P%U%T%S% %A% %(%M%A%T%R%I%X%)%,% %B% %(%M%A%T%R%I%X%)%,% %R%O%W%A%,% %C%O%L%A%,% %C%O%L%B%.% %O%U%T%P%U%T%S%:% %C% %(%M%A%T%R%I%X%)% %W%H%E%R%E% %C%=%A%*%B%. 7310 DIM C(ROWA,COLB) 7320 FOR I=1 TO ROWA 7330 FOR J=1 TO COLB 7340 FOR K=1 TO COLA 7350 LET C(I,J)=C(I,J)+A(I,K)*B(K,J) 7360 NEXT K 7370 NEXT J 7380 NEXT I 7390 RETURN 7400 REM %M%A%T%R%I%X% %T%R%A%N%S%P%O%S%I%T%I%O%N%.% %I%N%P%U%T%S%:% %A% %(%M%A%T%R%I%X%)%,% %R%O%W%A%,% %C%O%L%A%.% %O%U%T%P%U%T%S%:% %C% %(%M%A%T%R%I%X%)% %W%H%E%R%E% %C%=%A%"%. 7410 DIM C(ROWA,COLA) 7420 FOR I=1 TO ROWA 7430 FOR J=1 TO COLA 7440 LET C(I,J)=A(J,I) 7450 NEXT J 7460 NEXT I 7470 RETURN 7500 REM %M%A%T%R%I%X% %I%N%V%E%R%S%I%O%N%.% %I%N%P%U%T%S%:% %A% %(%M%A%T%R%I%X%)%,% %R%O%W%A%.% %O%U%T%P%U%T%S%:% %C% %(%M%A%T%R%I%X%)% %W%H%E%R%E% %C%=%A%*%*%-%1%. 7502 LET E$="" 7505 DIM C(ROWA,ROWA) 7510 FOR J=1 TO ROWA 7515 LET C(J,J)=1 7520 NEXT J 7525 FOR J=1 TO ROWA 7530 FOR I=J TO ROWA 7535 IF A(I,J)<>0 THEN GOTO 7555 7540 NEXT I 7545 LET E$="SINGULAR MATRIX" 7550 RETURN 7555 FOR K=1 TO ROWA 7560 LET TEMP=A(J,K) 7565 LET A(J,K)=A(I,K) 7570 LET A(I,K)=TEMP 7575 LET TEMP=C(J,K) 7580 LET C(J,K)=C(I,K) 7585 LET C(I,K)=TEMP 7590 NEXT K 7595 LET TEMP=1/A(J,J) 7600 FOR K=1 TO ROWA 7605 LET A(J,K)=TEMP*A(J,K) 7610 LET C(J,K)=TEMP*C(J,K) 7615 NEXT K 7620 FOR L=1 TO ROWA 7625 IF L=J THEN GOTO 7655 7630 LET TEMP=-A(L,J) 7635 FOR K=1 TO ROWA 7640 LET A(L,K)=A(L,K)+TEMP*A(J,K) 7645 LET C(L,K)=C(L,K)+TEMP*C(J,K) 7650 NEXT K 7655 NEXT L 7660 NEXT J 7665 RETURN 8000 SAVE "REGRESSIO%N" 8010 STOP 9000 REM %D%E%L%E%T%E%S% %L%I%N%E%S% %O%F% %P%R%O%G 9005 CLS 9007 PRINT "FIRST LINE TO BE DELETED="; 9010 INPUT FIRSTLINE 9012 PRINT FIRSTLINE 9015 PRINT "LAST LINE TO BE DELETED="; 9020 INPUT LASTLINE 9022 PRINT LASTLINE 9030 PRINT 9032 PRINT "CORRECT?" 9034 INPUT Z$ 9036 IF Z$(1)<>"Y" THEN GOTO 9000 9040 LET PROGBYTE=16509 9045 LET DFILE=PEEK 16396+256*PEEK 16397 9050 IF (256*PEEK PROGBYTE+PEEK (PROGBYTE+1))=FIRSTLINE THEN GOTO 9080 9060 LET PROGBYTE=PROGBYTE+4+PEEK (PROGBYTE+2)+256*PEEK (PROGBYTE+3) 9065 IF PROGBYTE>=DFILE THEN GOTO 9200 9070 GOTO 9050 9080 LET FIRSTLINE=PROGBYTE 9090 IF (256*PEEK PROGBYTE+PEEK (PROGBYTE+1))>LASTLINE THEN GOTO 9120 9100 LET PROGBYTE=PROGBYTE+4+PEEK (PROGBYTE+2)+256*PEEK (PROGBYTE+3) 9105 IF PROGBYTE>=DFILE THEN GOTO 9115 9110 GOTO 9090 9115 LET PROGBYTE=DFILE 9120 LET LASTLINE=PROGBYTE 9130 POKE (FIRSTLINE+2),(LASTLINE-FIRSTLINE-4)-256*INT ((LASTLINE-FIRSTLINE-4)/256) 9140 POKE (FIRSTLINE+3),INT ((LASTLINE-FIRSTLINE-4)/256) 9150 STOP 9200 PRINT "ERROR: FIRST LINE ";FIRSTLINE;" NOT FOUND." 9210 STOP