10 CLS
20 DEFDBL A-X
30 LLL = 0
40 DIM Y(10)
50 DIM CW(4)
60 Y = 0
70 LLL = 0
80 RANDOMIZE (TIMER / 100)
90 PRINT "Grey Hicks Speaker Dix Waters 4-04-05"
100 W11 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
110 W12 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
120 W13 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
130 W21 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
140 W22 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
150 W23 = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
160 CW(1) = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
170 CW(2) = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
180 CW(3) = INT(RND * 100) * ((-1) ^ (INT(RND * 100)))
190 W31 = CW(1)
200 W32 = CW(2)
210 W33 = CW(3)
220 PRINT W11, W21, W31, LLL
230 PRINT W12, W22, W32
240 PRINT W13, W23, W33
250 PLAY "mbo4cccl2el4o3bgeafdfaL2gL4"
260 FOR Z = 1 TO 8000
270 NEXT Z
280 Q11 = W11
290 Q12 = W12
300 LE = 0
310 Q13 = W13
320 Q21 = W21
330 Q22 = W22
340 Q23 = W23
350 Q31 = W31
360 Q32 = W32
370 Q33 = W33
380 K = 0
390 AU11 = Q11
400 A12 = Q12
410 A13 = Q13
420 A21 = Q21
430 AU22 = Q22
440 A23 = Q23
450 A31 = Q31
460 A32 = Q32
470 AU33 = Q33
480 A11 = AU11 - AU33
490 A22 = AU22 - AU33
500 A33 = AU33 - AU33
510 R11 = ((A11 ^ 2) + (A12 ^ 2) + (A13 ^ 2)) ^ (1 / 2)
520 A11 = A11 / R11
530 A12 = A12 / R11
540 A13 = A13 / R11
550 R12 = (A11 * A21) + (A12 * A22) + (A13 * A23)
560 R13 = (A11 * A31) + (A11 * A32) + (A11 * A33)
570 A21 = A21 - (A11 * R12)
580 A22 = A22 - (A12 * R12)
590 A23 = A23 - (A13 * R12)
600 A31 = A31 - (A11 * R13)
610 A32 = A32 - (A12 * R13)
620 A33 = A33 - (A13 * R13)
630 R22 = ((A21 ^ 2) + (A22 ^ 2) + (A23 ^ 2)) ^ (1 / 2)
640 A21 = A21 / R22
650 A22 = A22 / R22
660 A23 = A23 / R22
670 R23 = ((A21 * A31) + (A22 * A32) + (A23 * A33))
680 A31 = A31 - (A21 * R23)
690 A32 = A32 - (A22 * R23)
700 A33 = A33 - (A23 * R23)
710 R33 = ((A31 ^ 2) + (A32 ^ 2) + (A33 ^ 2)) ^ (1 / 2)
715 IF ABS(R33) < .0000000000000005# THEN K = 100
716 IF ABS(R33) < .0000000000000005# THEN GOTO 930
720 A31 = A31 / R33
730 A32 = A32 / R33
740 A33 = A33 / R33
750 B11 = ((R11 * A11) + (R12 * A12) + (R13 * A13)) + AU33
760 B21 = (R11 * A21) + (R12 * A22) + (R13 * A23)
770 B31 = (R11 * A31) + (R12 * A32) + (R13 * A33)
780 B12 = (0) + (R22 * A12) + (R23 * A13)
790 B22 = ((0) + (R22 * A22) + (R23 * A23)) + AU33
800 B32 = (0) + (R22 * A32) + (R23 * A33)
810 B13 = (0) + (0) + (R33 * A13)
820 B23 = R33 * A23
830 B33 = (R33 * A33) + AU33
840 Q11 = B11
850 Q12 = B12
860 Q13 = B13
870 Q21 = B21
880 Q22 = B22
890 Q23 = B23
900 Q31 = B31
910 Q32 = B32
920 Q33 = B33
930 K = K + 1
940 IF K < 99 THEN GOTO 390 ELSE GOTO 950
950 IF (CINT(B11)) <> (CINT(B33)) THEN GOTO 990
960 FF = B11
970 B11 = B22
980 B22 = FF
990 IF (CINT(B11)) <> (CINT(B22)) THEN GOTO 1030
1000 FFF = B11
1010 B11 = B33
1020 B33 = FFF
1030 PRINT "EIGENVALUES ARE"
1040 IF (CINT(B22)) = (CINT(B33)) THEN B22 = (B22 + B33) / 2
1050 IF (CINT(B22)) = (CINT(B33)) THEN B33 = B22
1060 PRINT B11, B22, B33
1070 IF LLL > 2 THEN GOTO 1130
1080 IF (CINT(B22)) = (CINT(B33)) THEN GOTO 1130
1090 LLL = LLL + 1
1100 CW(LLL - 1) = CW(LLL - 1) * -1
1110 CW(LLL) = CW(LLL) * -1
1120 GOTO 190
1130 XXX = 0
1140 PLAY "mbo3dedego4co3bo4cdc"
1150 PLAY "mbo3bo4cdl2cl4"
1160 FOR ZZ = 1 TO 8000
1170 NEXT ZZ
1180 H11 = W11 - B11
1190 H12 = W21
1200 H13 = W31 + 0
1210 H21 = W12
1220 H22 = W22 - B11
1230 H23 = W32 + 0
1240 H31 = W13
1250 H32 = W23
1260 H33 = W33 + ((-1) * B11) + (1 / 3)
1270 N11 = H11 * ((H22 * H33) - (H23 * H32))
1280 N12 = (-1 * H12) * ((H21 * H33) - (H23 * H31))
1290 N13 = H13 * ((H21 * H32) - (H22 * H31))
1300 N15 = N11 + N12 + N13
1310 G11 = ((-1) ^ 2) * ((H22 * H33) - (H23 * H32))
1320 G12 = ((-1) ^ 3) * ((H21 * H33) - (H23 * H31))
1330 G13 = ((-1) ^ 4) * ((H21 * H32) - (H22 * H31))
1340 G21 = ((-1) ^ 3) * ((H12 * H33) - (H13 * H32))
1350 G22 = ((-1) ^ 4) * ((H11 * H33) - (H13 * H31))
1360 G23 = ((-1) ^ 5) * ((H11 * H32) - (H12 * H31))
1370 G31 = ((-1) ^ 4) * ((H12 * H23) - (H13 * H22))
1380 G32 = ((-1) ^ 5) * ((H11 * H23) - (H13 * H21))
1390 G33 = ((-1) ^ 6) * ((H11 * H22) - (H12 * H21))
1400 F11 = G11 / N15
1410 F21 = G12 / N15
1420 F31 = G13 / N15
1430 F12 = G21 / N15
1440 F22 = G22 / N15
1450 F32 = G23 / N15
1460 F13 = G31 / N15
1470 F23 = G32 / N15
1480 F33 = G33 / N15
1490 Y = Y + 1
1500 O(Y) = (F13) * (1 / 3)
1510 Y = Y + 1
1520 O(Y) = (F23) * (1 / 3)
1530 Y = Y + 1
1540 O(Y) = (F33) * (1 / 3)
1550 LE = LE + 1
1560 IF LE > 1.5 THEN GOTO 1590 ELSE GOTO 1570
1570 B11 = B22
1580 GOTO 1180
1590 B11 = B33
1600 IF LE > 2.5 THEN GOTO 1610 ELSE GOTO 1180
1610 PRINT "EIGENVECTORS BY COLUMN AND EIGENVALUE"
1620 PRINT O(1), O(4), O(7)
1630 PRINT O(2), O(5), O(8)
1640 PRINT O(3), O(6), O(9)
1650 PRINT "FROBENIUS NORM OF EIGENVECTORS"
1660 M = 0
1670 FOR Y = 1 TO 9
1680 M = M + ((O(Y)) ^ 2)
1690 NEXT Y
1700 PRINT M ^ (1 / 2)
1710 CCC = M ^ (1 / 2)
1720  H11 = O(1)
1730 H12 = O(2)
1740 H13 = O(3)
1750 H21 = O(4)
1760 H22 = O(5)
1770 H23 = O(6)
1780 H31 = O(7)
1790 H32 = O(8)
1800 H33 = (O(9)) + 1
1810 N11 = H11 * ((H22 * H33) - (H23 * H32))
1820 N12 = (-1 * H12) * ((H21 * H33) - (H23 * H31))
1830 N13 = H13 * ((H21 * H32) - (H22 * H31))
1840 N15 = N11 + N12 + N13
1850 G11 = ((-1) ^ 2) * ((H22 * H33) - (H23 * H32))
1860 G12 = ((-1) ^ 3) * ((H21 * H33) - (H23 * H31))
1870 G13 = ((-1) ^ 4) * ((H21 * H32) - (H22 * H31))
1880 G21 = ((-1) ^ 3) * ((H12 * H33) - (H13 * H32))
1890 G22 = ((-1) ^ 4) * ((H11 * H33) - (H13 * H31))
1900 G23 = ((-1) ^ 5) * ((H11 * H32) - (H12 * H31))
1910 G31 = ((-1) ^ 4) * ((H12 * H23) - (H13 * H22))
1920 G32 = ((-1) ^ 5) * ((H11 * H23) - (H13 * H21))
1930 G33 = ((-1) ^ 6) * ((H11 * H22) - (H12 * H21))
1940 F11 = G11 / N15
1950 F21 = G12 / N15
1960 F31 = G13 / N15
1970 F12 = G21 / N15
1980 F22 = G22 / N15
1990 F32 = G23 / N15
2000 F13 = G31 / N15
2010 F23 = G32 / N15
2020 F33 = G33 / N15
2030 PRINT F13; F23; F33
2040 W = CINT(F13 * 10)
2050 IF W > -9.000001E-06 THEN GOSUB 2410
2060 IF W < 0 THEN GOSUB 2640
2070 W = CINT(F23 * 10)
2080 IF W > -9.000001E-06 THEN GOSUB 2410
2090 IF W < 0 THEN GOSUB 2640
2100 W = CINT(F33 * 10)
2110 IF W > -9.000001E-06 THEN GOSUB 2410
2120 IF W < 0 THEN GOSUB 2640
2130 UU = ((F11 ^ 2) + (F21 ^ 2) + (F31 ^ 2) + (F12 ^ 2) + (F22 ^ 2) + (F32 ^ 2) + (F13 ^ 2) + (F23 ^ 2) + (F33 ^ 2)) ^ (1 / 2)
2140 PRINT "FROBENIUS NORM OF INVERSE TRANSVERSE EIGENVECTORS +1 TO 3,3 & /4.5 &/1.81"
2150 PRINT UU, UU / 4.5, UU / 1.81
2160 PRINT " 1                                 2"
2170 PRINT " 3 4                     5 6 7 8 910   "; W11; W21; W31; LLL
2180 PRINT "1112                    131415161718   "; W12; W22; W32
2190 PRINT "192021222324252627282930313233343536   "; W13; W23; W33
2200 PRINT "373839404142434445464748495051525354"
2210 PRINT "555671727374757677787980818283848586"
2220 PRINT "8788103"
2230 PRINT "    5758596061626364656667686970"
2240 PRINT "    8990919293949596979899100101102"
2250 FOR ZZ = 1 TO 7000000
2260 NEXT ZZ
2270 OPEN "E.TXT" FOR OUTPUT AS 1
2280 PRINT #1, W11, W21, W31
2290 PRINT #1, W12, W22, W32
2300 PRINT #1, W13, W23, W33
2310 CLOSE 1
2315 GOTO 2390
2320 OPEN "SPTRIX.TXT" FOR APPEND AS 1
2330 PRINT #1, W11, W21, W31
2340 PRINT #1, W12, W22, W32
2350 PRINT #1, W13, W23, W33
2360 PRINT #1, CCC
2370 PRINT #1, UU
2380 CLOSE 1
2390 SYSTEM
2400 INPUT W
2410 W = CINT(W)
2420 ZZZ = 0
2430 Y = 0
2440 X = W
2450 F = 0
2460 A = 0
2470 C = 0
2480 IF X = 0 THEN Y = 0
2490 IF X = 0 THEN GOTO 2610
2500 IF X = 0 THEN GOTO 2610
2510 IF 2 ^ A > X THEN GOTO 2560
2520 C = 10 ^ A
2530 E = 2 ^ A
2540 A = A + 1
2550 GOTO 2510
2560 Y = Y + C
2570 A = 0
2580 F = F + E
2590 X = W - F
2600 GOTO 2500
2610 PRINT Y
2620 RETURN
2630 INPUT W
2640 X = ABS(W)
2650 ZZZ = 0
2660 Y = 0
2670 A = 0
2680 C = 0
2690 F = 0
2700 W = ABS(W)
2710 IF X = 0 THEN Y = 0
2720 IF X = 0 THEN GOTO 2960
2730 IF 2 ^ A < X THEN GOTO 2820
2740 EG = 2 ^ A
2750 AR = EG - 1
2760 A = 0
2770 X = AR
2780 C = 0
2790 Y = 0
2800 F = 0
2810 GOTO 2840
2820 A = A + 1
2830 GOTO 2730
2840 IF X = 0 THEN GOTO 2950
2850 IF 2 ^ A > X THEN GOTO 2900
2860 C = 10 ^ A
2870 E = 2 ^ A
2880 A = A + 1
2890 GOTO 2850
2900 Y = Y + C
2910 A = 0
2920 F = F + E
2930 X = AR - F
2940 GOTO 2840
2950 IF ZZZ > 0 THEN GOTO 2980
2960 BO = 1111111! - Y
2970 ZZZ = 1
2980 IF ZZZ = 2 THEN GOTO 3060
2990 X = EG - W
3000 AR = EG - W
3010 F = 0
3020 Y = 0
3030 A = 0
3040 ZZZ = 2
3050 GOTO 2840
3060 BO = BO + Y
3070 A$ = STR$(BO)
3080 PRINT "111111111"; MID$(A$, 2, 7)
3090 RETURN