10 CLS 20 DEFDBL A-X 25 DIM Y(10) 27 Y = 0 40 PRINT "Input 3x3 matrix row by row" 41 INPUT ; " 1,1"; W11 42 INPUT ; " 1,2"; W21 43 INPUT " 1,3"; W31 44 INPUT ; " 2,1"; W12 45 INPUT ; " 2,2"; W22 46 INPUT " 2,3"; W32 47 INPUT ; " 3,1"; W13 48 INPUT ; " 3,2"; W23 49 INPUT " 3,3"; W33 140 PRINT W11, W21, W31 150 PRINT W12, W22, W32 160 PRINT W13, W23, W33 161 PLAY "mbo4cccl2el4o3bgeafdfal2gl4" 170 FOR Z = 1 TO 8000 180 NEXT Z 190 Q11 = W11 200 Q12 = W12 210 LE = 0 220 Q13 = W13 230 Q21 = W21 240 Q22 = W22 250 Q23 = W23 260 Q31 = W31 270 Q32 = W32 280 Q33 = W33 290 K = 0 300 AU11 = Q11 310 A12 = Q12 320 A13 = Q13 330 A21 = Q21 340 AU22 = Q22 350 A23 = Q23 360 A31 = Q31 370 A32 = Q32 380 AU33 = Q33 390 A11 = AU11 - AU33 400 A22 = AU22 - AU33 410 A33 = AU33 - AU33 420 R11 = ((A11 ^ 2) + (A12 ^ 2) + (A13 ^ 2)) ^ (1 / 2) 430 A11 = A11 / R11 440 A12 = A12 / R11 450 A13 = A13 / R11 460 R12 = (A11 * A21) + (A12 * A22) + (A13 * A23) 470 R13 = (A11 * A31) + (A11 * A32) + (A11 * A33) 480 A21 = A21 - (A11 * R12) 490 A22 = A22 - (A12 * R12) 500 A23 = A23 - (A13 * R12) 510 A31 = A31 - (A11 * R13) 520 A32 = A32 - (A12 * R13) 530 A33 = A33 - (A13 * R13) 540 R22 = ((A21 ^ 2) + (A22 ^ 2) + (A23 ^ 2)) ^ (1 / 2) 550 A21 = A21 / R22 560 A22 = A22 / R22 570 A23 = A23 / R22 580 R23 = ((A21 * A31) + (A22 * A32) + (A23 * A33)) 590 A31 = A31 - (A21 * R23) 600 A32 = A32 - (A22 * R23) 610 A33 = A33 - (A23 * R23) 620 R33 = ((A31 ^ 2) + (A32 ^ 2) + (A33 ^ 2)) ^ (1 / 2) 625 IF ABS(R33) < .0000000000000005# THEN K = 100 626 IF ABS(R33) < .0000000000000005# THEN GOTO 840 630 A31 = A31 / R33 640 A32 = A32 / R33 650 A33 = A33 / R33 660 B11 = ((R11 * A11) + (R12 * A12) + (R13 * A13)) + AU33 670 B21 = (R11 * A21) + (R12 * A22) + (R13 * A23) 680 B31 = (R11 * A31) + (R12 * A32) + (R13 * A33) 690 B12 = (0) + (R22 * A12) + (R23 * A13) 700 B22 = ((0) + (R22 * A22) + (R23 * A23)) + AU33 710 B32 = (0) + (R22 * A32) + (R23 * A33) 720 B13 = (0) + (0) + (R33 * A13) 730 B23 = R33 * A23 740 B33 = (R33 * A33) + AU33 750 Q11 = B11 760 Q12 = B12 770 Q13 = B13 780 Q21 = B21 790 Q22 = B22 800 Q23 = B23 810 Q31 = B31 820 Q32 = B32 830 Q33 = B33 840 K = K + 1 850 IF K < 99 THEN GOTO 300 ELSE GOTO 851 851 IF (CINT(B11)) <> (CINT(B33)) THEN GOTO 856 852 FF = B11 853 B11 = B22 854 B22 = FF 856 IF (CINT(B11)) <> (CINT(B22)) THEN GOTO 860 857 FFF = B11 .858: : : : : : : : : : : B11 = B33 859 B33 = FFF 860 PRINT "EIGENVALUES ARE" 861 IF (CINT(B22)) = (CINT(B33)) THEN B22 = (B22 + B33) / 2 862 IF (CINT(B22)) = (CINT(B33)) THEN B33 = B22 870 PRINT B11, B22, B33 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 1.250 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 14.30 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" 2241 OPEN "E.TXT" FOR OUTPUT AS 2 2242 PRINT #2, W11, W21, W31 2243 PRINT #2, W12, W22, W32 2244 PRINT #2, W13, W23, W33 2245 CLOSE 2 2250 FOR ZZ = 1 TO 5000000 2260 NEXT ZZ 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