1 lll = 0
2  CLS
3  SCREEN 7
5  WINDOW (-20, 35)-(20, -20)
10 REM DEFINE DX/DT=D1=F(T,X,Y,Z) IN LINE 140
15 REM DEFINE DY/DT=D2=F(T,X,Y,Z) IN LINE 150
17 REM DEFINE DZ/DT=D3=F(T,X,Y,Z) IN LINE 155
18 PRINT "D1 = -Y - Z"
19 PRINT "D2 = x + Y / 5"
20 PRINT "D3 = 1 / 5 + Z * (x - 5.7)"
24 r = 1
45 d = 550
47 T1 = 1
48 T2 = 200
50 x = 0
51 Y = 0
52 Z = 0
53 N = 20000
110 FOR T9 = T1 TO T2 STEP d
130 FOR T = T9 TO T9 + d STEP d / N
135 COLOR (RND(r) * 10)
140 D1 = -Y - Z
150 D2 = x + Y / 5
155 D3 = 1 / 5 + Z * (x - 5.7)
160 x = x + D1 * d / N
170 Y = Y + D2 * d / N
175 Z = Z + D3 * d / N
176 PSET (x, Y + Z)
177 FOR w = 1 TO 1000
178 NEXT w
180 NEXT T
190 NEXT T9
CLS
SCREEN 7
WINDOW (-1.75, 1.4)-(1.75, -1.4)
COLOR (15)
PRINT "x = SIN(m) * COS(N)"
PRINT "Y = ((SIN(m)) * (((-1) + (2 / 3.141592654)) * ((SIN(N))))) + COS(m)"
FOR j = 0 TO 15 STEP .5
COLOR (j / 1.5)
FOR K = 0 TO 50 STEP .5
N = (K * 2) * (3.141592654# / 50)
m = j * (3.141596 / 15)
x = SIN(m) * COS(N)
Y = ((SIN(m)) * (((-1) + (2 / 3.141592654#)) * ((SIN(N))))) + COS(m)
PSET (x, Y)
FOR wo1 = 1 TO 1000
NEXT wo1
NEXT K
NEXT j
FOR i = 1 TO 10000000
NEXT i
CLS
SCREEN 1
WINDOW (-22, 15)-(22, -15)
CLS
RANDOMIZE (TIMER / 100)
FOR x = -10 TO 10 STEP .1
Y = (100 - (x ^ 2)) ^ (1 / 2)
PSET (x, Y)
Y = -((100 - (x ^ 2)) ^ (1 / 2))
PSET (x, Y)
NEXT x
a = RND(a)
a = 300 * a
PRINT "degrees =", a
PRINT "perpendicular = (-b / TAN(an)) + (10 * (1 / (COS((3.14159265 / 2) - (an)))))"
an = a * (3.14159265# / 180)
FOR b = -13 TO 13 STEP .1
c = b * TAN(an)
PSET (b, c)
p = (-b / TAN(an)) + (10 * (1 / (COS((3.14159265# / 2) - (an)))))
PSET (b, p)
NEXT b
FOR i = 1 TO 10000000 STEP 1
NEXT i
CLS
SCREEN 7
WINDOW (-1.05, 1.05)-(1.05, -1.05)
DIM a(25)
DIM b(25)
DIM c(138)
DIM e(138)
DIM O(258)
DIM N(258)
DIM w(1025)
DIM u(1025)
DIM q1(1025)
DIM q2(1025)
a(1) = -.5
a(2) = -.5
a(3) = -.5
a(4) = .5
a(5) = .5
b(1) = -.5
b(2) = -.5
b(3) = .5
b(4) = .5
b(5) = -.5
begin:
CLS
FOR K = 2 TO 5
COLOR (K)
LINE (a(K - 1), b(K - 1))-(a(K), b(K))
FOR r = 1 TO 300000
NEXT r
NEXT K
FOR r = 1 TO 300000
NEXT r
CLS
PRINT "itration 2 to 17 by 4"
FOR m = 1 TO 300000
NEXT m
FOR g = 2 TO 5
e(g) = (-((a(g)) * (-1)) / 2) - .5
c(g) = (b(g) / 2) - .5
c(1) = c(2)
e(1) = e(2)
LINE ((c(g - 1)), (e(g - 1)))-(c(g), e(g))
FOR r = 1 TO 300000
NEXT r
NEXT g
FOR g = 6 TO 9
COLOR (g)
c(g) = ((a(g - 4)) / 2) - .5
e(g) = (b(g - 4) / 2) + .5
c(1) = c(2)
e(1) = e(2)
LINE ((c(g - 1)), (e(g - 1)))-(c(g), e(g))
FOR r = 1 TO 300000
NEXT r
NEXT g
FOR g = 10 TO 13
COLOR (g - 5)
c(g) = ((a(g - 8)) / 2) + .5
e(g) = (b(g - 8) / 2) + .5
c(1) = c(2)
e(1) = e(2)
LINE ((c(g - 1)), (e(g - 1)))-(c(g), e(g))
FOR r = 1 TO 300000
NEXT r
NEXT g
FOR g = 14 TO 17
COLOR (g - 10)
e(g) = (((a(g - 12)) * (-1)) / 2) - .5
c(g) = -(b(g - 12) / 2) + .5
c(1) = c(2)
e(1) = e(2)
LINE ((c(g - 1)), (e(g - 1)))-(c(g), e(g))
FOR r = 1 TO 300000
NEXT r
NEXT g
FOR r = 1 TO 300000
NEXT r
CLS
PRINT "iteration 2 to 65 by 15"
FOR m = 1 TO 300000
NEXT m
FOR g = 2 TO 17
COLOR (g / 2)
N(g) = (-((c(g)) * (-1)) / 2) - .5
O(g) = (e(g) / 2) - .5
O(1) = O(2)
N(1) = N(2)
LINE ((O(g - 1)), (N(g - 1)))-(O(g), N(g))
NEXT g
FOR g = 18 TO 33
COLOR ((g / 2) - 8)
O(g) = ((c(g - 16)) / 2) - .5
N(g) = (e(g - 16) / 2) + .5
O(1) = O(2)
N(1) = N(2)
LINE ((O(g - 1)), (N(g - 1)))-(O(g), N(g))
NEXT g
FOR g = 34 TO 49
COLOR ((g / 4) - 7)
O(g) = ((c(g - 32)) / 2) + .5
N(g) = (e(g - 32) / 2) + .5
O(1) = O(2)
N(1) = N(2)
LINE ((O(g - 1)), (N(g - 1)))-(O(g), N(g))
NEXT g
FOR g = 50 TO 65
COLOR ((g / 5) - 8)
N(g) = (((c(g - 48)) * (-1)) / 2) - .5
O(g) = -(e(g - 48) / 2) + .5
O(1) = O(2)
N(1) = N(2)
LINE ((O(g - 1)), (N(g - 1)))-(O(g), N(g))
NEXT g
FOR r = 1 TO 300000
NEXT r
CLS
PRINT "iteration 2 to 257 by 63"
FOR m = 1 TO 300000
NEXT m
FOR g = 2 TO 65
COLOR (RND(r) * 10)
u(g) = (-((O(g)) * (-1)) / 2) - .5
w(g) = (N(g) / 2) - .5
w(1) = w(2)
u(1) = u(2)
LINE ((w(g - 1)), (u(g - 1)))-(w(g), u(g))
NEXT g
FOR g = 66 TO 129
COLOR (RND(r) * 10)
w(g) = ((O(g - 64)) / 2) - .5
u(g) = (N(g - 64) / 2) + .5
w(1) = w(2)
u(1) = u(2)
LINE ((w(g - 1)), (u(g - 1)))-(w(g), u(g))
NEXT g
FOR g = 130 TO 193
COLOR (RND(r) * 10)
w(g) = ((O(g - 128)) / 2) + .5
u(g) = (N(g - 128) / 2) + .5
w(1) = w(2)
u(1) = u(2)
LINE ((w(g - 1)), (u(g - 1)))-(w(g), u(g))
NEXT g
FOR g = 194 TO 257
COLOR (RND(r) * 10)
u(g) = (((O(g - 192)) * (-1)) / 2) - .5
w(g) = -(N(g - 192) / 2) + .5
w(1) = w(2)
u(1) = u(2)
LINE ((w(g - 1)), (u(g - 1)))-(w(g), u(g))
FOR r = 1 TO 20
NEXT r
NEXT g
FOR r = 1 TO 300000
NEXT r
CLS
PRINT "iteration 2 to 1025 by 255"
FOR m = 1 TO 300000
NEXT m
FOR g = 2 TO 257
COLOR (RND(r) * 10)
q2(g) = (-((w(g)) * (-1)) / 2) - .5
q1(g) = (u(g) / 2) - .5
q1(1) = q1(2)
q2(1) = q2(2)
LINE ((q1(g - 1)), (q2(g - 1)))-(q1(g), q2(g))
NEXT g
FOR g = 258 TO 513
COLOR (RND(r) * 10)
q1(g) = ((w(g - 256)) / 2) - .5
q2(g) = (u(g - 256) / 2) + .5
q1(1) = q1(2)
q2(1) = q2(2)
LINE ((q1(g - 1)), (q2(g - 1)))-(q1(g), q2(g))
NEXT g
FOR g = 514 TO 769
COLOR (RND(r) * 10)
q1(g) = ((w(g - 512)) / 2) + .5
q2(g) = (u(g - 512) / 2) + .5
q1(1) = q1(2)
q2(1) = q2(2)
LINE ((q1(g - 1)), (q2(g - 1)))-(q1(g), q2(g))
NEXT g
FOR g = 770 TO 1025
COLOR (RND(r) * 10)
q2(g) = (((w(g - 768)) * (-1)) / 2) - .5
q1(g) = -(u(g - 768) / 2) + .5
q1(1) = q1(2)
q2(1) = q2(2)
LINE ((q1(g - 1)), (q2(g - 1)))-(q1(g), q2(g))
FOR r = 1 TO 2000
NEXT r
NEXT g
FOR r = 1 TO 1000000
NEXT r
CLS
SCREEN 1
WINDOW (-10, 10)-(10, -10)
100 CLS
DEFDBL B-C, G, L-N, X
DIM b1(4, 4)
DIM b5(5, 5)
RANDOMIZE (TIMER / 100)
a = 1
b = RND * 10
c = 2
d = RND * 10
e = 3
f = RND * 10
g = 4
h = RND * 10
i = 5
j = RND * 10
PRINT b, d, f, h, j
LINE (-10, 0)-(10, 0)
LINE (0, 10)-(0, -10)
PSET (a, b)
PSET (c, d)
PSET (e, f)
PSET (g, h)
PSET (i, j)
FOR amm = 1 TO 1000000
NEXT amm
b5(1, 1) = 5
b5(1, 2) = a + c + e + g + i
b5(1, 3) = (a ^ 2) + (c ^ 2) + (e ^ 2) + (g ^ 2) + (i ^ 2)
b5(1, 4) = (a ^ 3) + (c ^ 3) + (e ^ 3) + (g ^ 3) + (i ^ 3)
b5(1, 5) = (a ^ 4) + (c ^ 4) + (e ^ 4) + (g ^ 4) + (i ^ 4)
b5(2, 1) = b5(1, 2)
b5(2, 2) = b5(1, 3)
b5(2, 3) = (a ^ 3) + (c ^ 3) + (e ^ 3) + (g ^ 3) + (i ^ 3)
b5(2, 4) = (a ^ 4) + (c ^ 4) + (e ^ 4) + (g ^ 4) + (i ^ 4)
b5(2, 5) = (a ^ 5) + (c ^ 5) + (e ^ 5) + (g ^ 5) + (i ^ 5)
b5(3, 1) = b5(1, 3)
b5(3, 2) = b5(2, 3)
b5(3, 3) = (a ^ 4) + (c ^ 4) + (e ^ 4) + (g ^ 4) + (i ^ 4)
b5(3, 4) = (a ^ 5) + (c ^ 5) + (e ^ 5) + (g ^ 5) + (i ^ 5)
b5(3, 5) = (a ^ 6) + (c ^ 6) + (e ^ 6) + (g ^ 6) + (i ^ 6)
b5(4, 1) = b5(1, 4)
b5(4, 2) = b5(2, 4)
b5(4, 3) = b5(3, 4)
b5(4, 4) = (a ^ 6) + (c ^ 6) + (e ^ 6) + (g ^ 6) + (i ^ 6)
b5(4, 5) = (a ^ 7) + (c ^ 7) + (e ^ 7) + (g ^ 7) + (i ^ 7)
b5(5, 1) = b5(1, 5)
b5(5, 2) = b5(2, 5)
b5(5, 3) = b5(3, 5)
b5(5, 4) = b5(4, 5)
b5(5, 5) = (a ^ 8) + (c ^ 8) + (e ^ 8) + (g ^ 8) + (i ^ 8)
g1(1, 1) = b + d + f + h + j
g1(1, 2) = (a * b) + (c * d) + (e * f) + (g * h) + (i * j)
g1(1, 3) = (b * (a ^ 2)) + (d * (c ^ 2)) + (f * (e ^ 2)) + (h * (g ^ 2)) + (j * (i ^ 2))
g1(1, 4) = (b * (a ^ 3)) + (d * (c ^ 3)) + (f * (e ^ 3)) + (h * (g ^ 3)) + (j * (i ^ 3))
g1(1, 5) = (b * (a ^ 4)) + (d * (c ^ 4)) + (f * (e ^ 4)) + (h * (g ^ 4)) + (j * (i ^ 4))
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 2 TO 5 STEP 1
b1((i1) - 1, (i2) - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
h51 = nc
m51 = b5(1, 1) * h51
FOR i1 = 2 TO 5 STEP 1
b1((i1) - 1, 1) = b5(i1, 1)
NEXT i1
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 3 TO 5 STEP 1
b1((i1) - 1, (i2) - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
h52 = nc
m52 = b5(1, 2) * h52
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 1 TO 2 STEP 1
b1((i1) - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 4 TO 5 STEP 1
b1((i1) - 1, (i2) - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
h53 = nc
        m53 = b5(1, 3) * h53
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 1 TO 3 STEP 1
b1((i1) - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 2 TO 5 STEP 1
b1((i1) - 1, 4) = b5(i1, 5)
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
h54 = nc
m54 = b5(1, 4) * h54
FOR i1 = 2 TO 5 STEP 1
FOR i2 = 1 TO 4 STEP 1
b1((i1) - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
h55 = nc
m55 = b5(1, 5) * h55
mt = m51 - m52 + m53 - m54 + m55
FOR i2 = 2 TO 5 STEP 1
b1(1, (i2) - 1) = b5(1, i2)
NEXT i2
FOR i1 = 3 TO 5
FOR i2 = 2 TO 5
b1(i1 - 1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c21 = -nc
b1(1, 1) = b5(1, 1)
b1(1, 2) = b5(1, 3)
b1(1, 3) = b5(1, 4)
b1(1, 4) = b5(1, 5)
FOR i1 = 3 TO 5
b1((i1) - 1, 1) = b5(i1, 1)
NEXT i1
FOR i1 = 3 TO 5
FOR i2 = 3 TO 5
b1((i1) - 1, (i2) - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c22 = nc
b1(1, 1) = b5(1, 1)
b1(1, 2) = b5(1, 2)
FOR i1 = 3 TO 5 STEP 1
FOR i2 = 1 TO 2 STEP 1
b1((i1) - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c23 = -nc
FOR i2 = 1 TO 3 STEP 1
b1(1, i2) = b5(1, i2)
NEXT i2
FOR i1 = 3 TO 5 STEP 1
FOR i2 = 1 TO 3 STEP 1
b1(i1 - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 3 TO 5 STEP 1
b1(i1 - 1, 4) = b5(i1, 5)
NEXT i1
b1(1, 4) = b5(1, 5)
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c24 = nc
FOR i2 = 1 TO 4
q = 1
b1(q, i2) = b5(q, i2)
NEXT i2
FOR i1 = 3 TO 5 STEP 1
FOR i2 = 1 TO 4 STEP 1
b1((i1) - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c25 = -nc
FOR i1 = 1 TO 2 STEP 1
FOR i2 = 2 TO 5 STEP 1
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 4 TO 5 STEP 1
FOR i2 = 2 TO 5 STEP 1
b1(i1 - 1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c31 = nc
b1(1, 1) = b5(1, 1)
b1(2, 1) = b5(2, 1)
FOR i1 = 1 TO 2
FOR i2 = 3 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 4 TO 5
FOR i2 = 3 TO 5
b1(i1 - 1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
b1(3, 1) = b5(4, 1)
b1(4, 1) = b5(5, 1)
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c32 = -nc
b1(1, 1) = b5(1, 1)
b1(2, 1) = b5(2, 1)
b1(1, 2) = b5(1, 2)
b1(2, 2) = b5(2, 2)
FOR i1 = 4 TO 5
FOR i2 = 1 TO 2
b1(i1 - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 1 TO 2
FOR i2 = 4 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 4 TO 5
FOR i2 = 4 TO 5
b1(i1 - 1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c33 = nc
FOR i1 = 1 TO 2
FOR i2 = 1 TO 3
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 4 TO 5
FOR i2 = 1 TO 3
b1(i1 - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 1 TO 2
b1(i1, 4) = b5(i1, 5)
NEXT i1
b1(3, 4) = b5(4, 5)
b1(4, 4) = b5(5, 5)
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c34 = -nc
FOR i1 = 1 TO 2
FOR i2 = 1 TO 4
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 4 TO 5
FOR i2 = 1 TO 4
b1(i1 - 1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c35 = nc
FOR i1 = 1 TO 3
FOR i2 = 2 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i2 = 2 TO 5
b1(4, i2 - 1) = b5(5, i2)
NEXT i2
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c41 = -nc
b1(1, 1) = b5(1, 1)
b1(2, 1) = b5(2, 1)
b1(3, 1) = b5(3, 1)
b1(4, 1) = b5(5, 1)
FOR i1 = 1 TO 3
FOR i2 = 3 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i2 = 3 TO 5
b1(4, i2 - 1) = b5(5, i2)
NEXT i2
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c42 = nc
FOR i1 = 1 TO 3
FOR i2 = 1 TO 2
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i2 = 1 TO 3
b1(4, i2) = b5(5, i2)
NEXT i2
FOR i1 = 1 TO 3
FOR i2 = 4 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
b1(4, 4) = b5(5, 5)
b1(4, 3) = b5(5, 4)
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c43 = -nc
FOR i1 = 1 TO 3
FOR i2 = 1 TO 3
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i2 = 1 TO 3
b1(4, i2) = b5(5, i2)
NEXT i2
FOR i1 = 1 TO 3
b1(i1, 4) = b5(i1, 5)
NEXT i1
b1(4, 4) = b5(5, 5)
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c44 = nc
FOR i1 = 1 TO 3
FOR i2 = 1 TO 4
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i2 = 1 TO 4
b1(4, i2) = b5(5, i2)
NEXT i2
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c45 = -nc
FOR i1 = 1 TO 4
FOR i2 = 2 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c51 = nc
FOR i1 = 1 TO 4
b1(i1, 1) = b5(i1, 1)
NEXT i1
FOR i1 = 1 TO 4
FOR i2 = 3 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c52 = -nc
FOR i1 = 1 TO 4
FOR i2 = 1 TO 2
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 1 TO 4
FOR i2 = 4 TO 5
b1(i1, i2 - 1) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c53 = nc
FOR i1 = 1 TO 4
FOR i2 = 1 TO 3
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
FOR i1 = 1 TO 4
b1(i1, 4) = b5(i1, 5)
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c54 = -nc
FOR i1 = 1 TO 4
FOR i2 = 1 TO 4
b1(i1, i2) = b5(i1, i2)
NEXT i2
NEXT i1
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
c55 = nc
l(1, 1) = ((h51 / mt) * g1(1, 1)) + ((-h52 / mt) * g1(1, 2)) + ((h53 / mt) * g1(1, 3)) + ((-h54 / mt) * g1(1, 4)) + ((h55 / mt) * g1(1, 5))
l(2, 1) = ((c21 / mt) * g1(1, 1)) + ((c22 / mt) * g1(1, 2)) + ((c23 / mt) * g1(1, 3)) + ((c24 / mt) * g1(1, 4)) + ((c25 / mt) * g1(1, 5))
l(3, 1) = ((c31 / mt) * g1(1, 1)) + ((c32 / mt) * g1(1, 2)) + ((c33 / mt) * g1(1, 3)) + ((c34 / mt) * g1(1, 4)) + ((c35 / mt) * g1(1, 5))
l(4, 1) = ((c41 / mt) * g1(1, 1)) + ((c42 / mt) * g1(1, 2)) + ((c43 / mt) * g1(1, 3)) + ((c44 / mt) * g1(1, 4)) + ((c45 / mt) * g1(1, 5))
l(5, 1) = ((c51 / mt) * g1(1, 1)) + ((c52 / mt) * g1(1, 2)) + ((c53 / mt) * g1(1, 3)) + ((c54 / mt) * g1(1, 4)) + ((c55 / mt) * g1(1, 5))
FOR x = -10 TO 10 STEP .05
Y = ((l(5, 1)) * (x ^ 4)) + ((l(4, 1)) * (x ^ 3)) + ((l(3, 1)) * (x ^ 2)) + (l(2, 1) * x) + l(1, 1)
PSET (x, Y)
NEXT x
PRINT "Equation of line fitting these five plots is"
PRINT l(5, 1); "x^4"; l(4, 1); "x^3"; l(3, 1); "x^2"; l(2, 1); "x"; l(1, 1)
FOR s = 1 TO 2000000
NEXT s
CLS
SCREEN 7
WINDOW (-20, 50)-(20, -25)
REM DEFINE DX/DT=D1=F(T,X,Y,Z)
REM DEFINE DY/DT=D2=F(T,X,Y,Z)
REM DEFINE DZ/DT=D3=F(T,X,Y,Z)
PRINT "D1 = 10 * (Y - x)"
PRINT "D2 = 28 * x - x * Z"
PRINT "D3 = x * Y - (8 / 3) * Z"
d = 100
T1 = 1
T2 = 2
x = 5
Y = 5
Z = 5
N = 19000
FOR T9 = T1 TO T2 STEP d
FOR T = T9 TO T9 + d STEP d / N
COLOR (RND(r) * 10)
D1 = 10 * (Y - x)
D2 = 28 * x - x * Z
D3 = x * Y - (8 / 3) * Z
x = x + D1 * d / N
Y = Y + D2 * d / N
Z = Z + D3 * d / N
PSET (x, Y)
FOR so1 = 1 TO 1000
NEXT so1
NEXT T
NEXT T9
CLS
SCREEN 7
COLOR (7)
WINDOW (-1000, 1000)-(1000, -1000)
DEFDBL A-Z
PRINT "THIS PROGRAM TAKES A RANDOM 3 X 3 MATRIX AND CALCULATES"
PRINT "THE EIGENVALUES AND EIGENVECTORS FOR THE MATRIX. IT THEN"
PRINT "PLOTS ON A GRAPH THE CHARACTERISTIC   "
PRINT "EQUATION. THEN IT"
PRINT "TAKES SAMPLES AND REPLOTS AN EQUATION IN TERMS OF X OVER THE"
PRINT "CHARACTERISTIC PLOT."
FOR r = 1 TO 3000000
NEXT r
RANDOMIZE (TIMER - 100)
SCREEN 2
SCREEN 0
PRINT "matrix is"
Y = 0
w11 = RND * 10
w12 = RND * 10
w13 = RND * 10
w21 = RND * 10
w22 = RND * 10
w23 = RND * 10
w31 = RND * 10
w32 = RND * 10
w33 = RND * 10
PRINT w11, w21, w31
PRINT w12, w22, w32
PRINT w13, w23, w33
FOR m = 1 TO 100000
NEXT m
q11 = w11
q12 = w12
l = 0
q13 = w13
q21 = w21
q22 = w22
q23 = w23
q31 = w31
q32 = w32
q33 = w33
K = 0
abegin:
au11 = q11
a12 = q12
a13 = q13
a21 = q21
au22 = q22
a23 = q23
a31 = q31
a32 = q32
au33 = q33
a11 = au11 - au33
a22 = au22 - au33
a33 = au33 - au33
r11 = ((a11 ^ 2) + (a12 ^ 2) + (a13 ^ 2)) ^ (1 / 2)
a11 = a11 / r11
a12 = a12 / r11
a13 = a13 / r11
r12 = (a11 * a21) + (a12 * a22) + (a13 * a23)
r13 = (a11 * a31) + (a11 * a32) + (a11 * a33)
a21 = a21 - (a11 * r12)
a22 = a22 - (a12 * r12)
a23 = a23 - (a13 * r12)
a31 = a31 - (a11 * r13)
a32 = a32 - (a12 * r13)
a33 = a33 - (a13 * r13)
r22 = ((a21 ^ 2) + (a22 ^ 2) + (a23 ^ 2)) ^ (1 / 2)
a21 = a21 / r22
a22 = a22 / r22
a23 = a23 / r22
r23 = ((a21 * a31) + (a22 * a32) + (a23 * a33))
a31 = a31 - (a21 * r23)
a32 = a32 - (a22 * r23)
a33 = a33 - (a23 * r23)
R33 = ((a31 ^ 2) + (a32 ^ 2) + (a33 ^ 2)) ^ (1 / 2)
IF ABS(R33) < .0000000000000005# THEN K = 100
IF ABS(R33) < .0000000000000005# THEN GOTO aaaaa
a31 = a31 / R33
a32 = a32 / R33
a33 = a33 / R33
B11 = ((r11 * a11) + (r12 * a12) + (r13 * a13)) + au33
b21 = (r11 * a21) + (r12 * a22) + (r13 * a23)
b31 = (r11 * a31) + (r12 * a32) + (r13 * a33)
b12 = (0) + (r22 * a12) + (r23 * a13)
B22 = ((0) + (r22 * a22) + (r23 * a23)) + au33
b32 = (0) + (r22 * a32) + (r23 * a33)
b13 = (0) + (0) + (R33 * a13)
b23 = R33 * a23
B33 = (R33 * a33) + au33
q11 = B11
q12 = b12
q13 = b13
q21 = b21
q22 = B22
q23 = b23
q31 = b31
q32 = b32
q33 = B33
660 K = K + 1
aaaaa:
670 IF K < 99 THEN GOTO abegin ELSE GOTO hello
hello:
PRINT "eigenvalues are"
PRINT B11, B22, B33
FOR m = 1 TO 2000000
NEXT m
PRINT "eigenvectors are"
eddy:
h11 = w11 - B11
h12 = w21
h13 = w31 + 0
h21 = w12
h22 = w22 - B11
h23 = w32 + 0
h31 = w13
h32 = w23
h33 = w33 + ((-1) * B11) + (1 / 3)
n11 = h11 * ((h22 * h33) - (h23 * h32))
n12 = (-1 * h12) * ((h21 * h33) - (h23 * h31))
n13 = h13 * ((h21 * h32) - (h22 * h31))
n15 = n11 + n12 + n13
g11 = ((-1) ^ 2) * ((h22 * h33) - (h23 * h32))
g12 = ((-1) ^ 3) * ((h21 * h33) - (h23 * h31))
g13 = ((-1) ^ 4) * ((h21 * h32) - (h22 * h31))
g21 = ((-1) ^ 3) * ((h12 * h33) - (h13 * h32))
g22 = ((-1) ^ 4) * ((h11 * h33) - (h13 * h31))
g23 = ((-1) ^ 5) * ((h11 * h32) - (h12 * h31))
g31 = ((-1) ^ 4) * ((h12 * h23) - (h13 * h22))
g32 = ((-1) ^ 5) * ((h11 * h23) - (h13 * h21))
g33 = ((-1) ^ 6) * ((h11 * h22) - (h12 * h21))
f11 = g11 / n15
f21 = g12 / n15
f31 = g13 / n15
f12 = g21 / n15
f22 = g22 / n15
f32 = g23 / n15
F13 = g31 / n15
F23 = g32 / n15
F33 = g33 / n15
F13 = F13
F23 = F23
F33 = F33
Y = Y + 1
O(Y) = (F13) * (1 / 3)
Y = Y + 1
O(Y) = (F23) * (1 / 3)
Y = Y + 1
O(Y) = (F33) * (1 / 3)
l = l + 1
IF l > 1.5 THEN GOTO try ELSE GOTO oil
oil:
B11 = B22
GOTO eddy
try:
B11 = B33
IF l > 2.5 THEN GOTO caow ELSE GOTO eddy
caow:
PRINT "EIGENVECTORS BY COLUMN AND EIGENVALUE"
PRINT O(1), O(4), O(7)
PRINT O(2), O(5), O(8)
PRINT O(3), O(6), O(9)
FOR m = 1 TO 5000000
NEXT m
CLS
FOR x = -25 TO 25 STEP .2
h11 = w11 - x
h12 = w21
h13 = w31 + 0
h21 = w12
h22 = w22 - x
h23 = w32 + 0
h31 = w13
h32 = w23
h33 = w33 - x
n11 = h11 * ((h22 * h33) - (h23 * h32))
n12 = (-1 * h12) * ((h21 * h33) - (h23 * h31))
n13 = h13 * ((h21 * h32) - (h22 * h31))
n15 = n11 + n12 + n13
SCREEN 7
COLOR (7)
WINDOW (-20, 35)-(20, -20)
PSET (x, n15)
NEXT x
LINE (0, 1000)-(0, -1000)
LINE (-100, 0)-(100, 0)
FOR m = 1 TO 2000
NEXT m
FOR m = 1 TO 1000
NEXT m
FOR x = 5 TO 25 STEP 5
h11 = w11 - x
h12 = w21
h13 = w31 + 0
h21 = w12
h22 = w22 - x
h23 = w32 + 0
h31 = w13
h32 = w23
h33 = w33 - x
n11 = h11 * ((h22 * h33) - (h23 * h32))
n12 = (-1 * h12) * ((h21 * h33) - (h23 * h31))
n13 = h13 * ((h21 * h32) - (h22 * h31))
n15 = n11 + n12 + n13
v(1, x / 5) = n15
NEXT x
a = 5
b = v(1, 1)
c = 10
d = v(1, 2)
e = 15
f = v(1, 3)
g = 20
h = v(1, 4)
i = 25
j = v(1, 5)
b0(1, 1) = 1
b0(1, 2) = 1
b0(1, 3) = 1
b0(1, 4) = 1
b0(1, 5) = 1
b0(2, 1) = a
b0(2, 2) = c
b0(2, 3) = e
b0(2, 4) = g
b0(2, 5) = i
b0(3, 1) = a ^ 2
b0(3, 2) = c ^ 2
b0(3, 3) = e ^ 2
b0(3, 4) = g ^ 2
b0(3, 5) = i ^ 2
b0(4, 1) = a ^ 3
b0(4, 2) = c ^ 3
b0(4, 3) = e ^ 3
b0(4, 4) = g ^ 3
b0(4, 5) = i ^ 3
b1(1, 1) = 5
b1(1, 2) = a + c + e + g + i
b1(1, 3) = (a ^ 2) + (c ^ 2) + (e ^ 2) + (g ^ 2) + (i ^ 2)
b1(1, 4) = (a ^ 3) + (c ^ 3) + (e ^ 3) + (g ^ 3) + (i ^ 3)
b1(2, 1) = b1(1, 2)
b1(2, 2) = b1(1, 3)
b1(2, 3) = (a ^ 3) + (c ^ 3) + (e ^ 3) + (g ^ 3) + (i ^ 3)
b1(2, 4) = (a ^ 4) + (c ^ 4) + (e ^ 4) + (g ^ 4) + (i ^ 4)
b1(3, 1) = b1(1, 3)
b1(3, 2) = b1(2, 3)
b1(3, 3) = (a ^ 4) + (c ^ 4) + (e ^ 4) + (g ^ 4) + (i ^ 4)
b1(3, 4) = (a ^ 5) + (c ^ 5) + (e ^ 5) + (g ^ 5) + (i ^ 5)
b1(4, 1) = b1(1, 4)
b1(4, 2) = b1(2, 4)
b1(4, 3) = b1(3, 4)
b1(4, 4) = (a ^ 6) + (c ^ 6) + (e ^ 6) + (g ^ 6) + (i ^ 6)
g1(1, 1) = b + d + f + h + j
g1(1, 2) = (a * b) + (c * d) + (e * f) + (g * h) + (i * j)
g1(1, 3) = (b * (a ^ 2)) + (d * (c ^ 2)) + (f * (e ^ 2)) + (h * (g ^ 2)) + (j * (i ^ 2))
g1(1, 4) = (b * (a ^ 3)) + (d * (c ^ 3)) + (f * (e ^ 3)) + (h * (g ^ 3)) + (j * (i ^ 3))
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
nx1 = b1(1, 1) * (n11 - n12 + n13)
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
nx2 = b1(1, 2) * (n21 - n22 + n23)
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx3 = b1(1, 3) * (n31 - n32 + n33)
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
nx4 = b1(1, 4) * (n41 - n42 + n43)
nc = nx1 - nx2 + nx3 - nx4
n11 = (b1(2, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n12 = (b1(2, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n13 = (b1(2, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
c(1, 1) = (n11 - n12 + n13) / nc
n21 = (b1(2, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n22 = (b1(2, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n23 = (b1(2, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
c(1, 2) = -(n21 - n22 + n23) / nc
n31 = (b1(2, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n32 = (b1(2, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n33 = (b1(2, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
c(1, 3) = (n31 - n32 + n33) / nc
n41 = (b1(2, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n42 = (b1(2, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n43 = (b1(2, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
c(1, 4) = -(n41 - n42 + n43) / nc
n211 = (b1(1, 2) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n212 = (b1(1, 3) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n213 = (b1(1, 4) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
c(2, 1) = -(n211 - n212 + n213) / nc
n221 = (b1(1, 1) * ((b1(3, 3) * b1(4, 4)) - (b1(3, 4) * b1(4, 3))))
n222 = (b1(1, 3) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n223 = (b1(1, 4) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
c(2, 2) = (n221 - n222 + n223) / nc
n231 = (b1(1, 1) * ((b1(3, 2) * b1(4, 4)) - (b1(3, 4) * b1(4, 2))))
n232 = (b1(1, 2) * ((b1(3, 1) * b1(4, 4)) - (b1(3, 4) * b1(4, 1))))
n233 = (b1(1, 4) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
c(2, 3) = -(n231 - n232 + n233) / nc
n241 = (b1(1, 1) * ((b1(3, 2) * b1(4, 3)) - (b1(3, 3) * b1(4, 2))))
n242 = (b1(1, 2) * ((b1(3, 1) * b1(4, 3)) - (b1(3, 3) * b1(4, 1))))
n243 = (b1(1, 3) * ((b1(3, 1) * b1(4, 2)) - (b1(3, 2) * b1(4, 1))))
c(2, 4) = (n241 - n242 + n243) / nc
n311 = (b1(1, 2) * ((b1(2, 3) * b1(4, 4)) - (b1(2, 4) * b1(4, 3))))
n312 = (b1(1, 3) * ((b1(2, 2) * b1(4, 4)) - (b1(2, 4) * b1(4, 2))))
n313 = (b1(1, 4) * ((b1(2, 2) * b1(4, 3)) - (b1(2, 3) * b1(4, 2))))
c(3, 1) = (n311 - n312 + n313) / nc
n321 = (b1(1, 1) * ((b1(2, 3) * b1(4, 4)) - (b1(2, 4) * b1(4, 3))))
n322 = (b1(1, 3) * ((b1(2, 1) * b1(4, 4)) - (b1(2, 4) * b1(4, 1))))
n323 = (b1(1, 4) * ((b1(2, 1) * b1(4, 3)) - (b1(2, 3) * b1(4, 1))))
c(3, 2) = -(n321 - n322 + n323) / nc
n331 = (b1(1, 1) * ((b1(2, 2) * b1(4, 4)) - (b1(2, 4) * b1(4, 2))))
n332 = (b1(1, 2) * ((b1(2, 1) * b1(4, 4)) - (b1(2, 4) * b1(4, 1))))
n333 = (b1(1, 4) * ((b1(2, 1) * b1(4, 2)) - (b1(2, 2) * b1(4, 1))))
c(3, 3) = (n331 - n332 + n333) / nc
n341 = (b1(1, 1) * ((b1(2, 2) * b1(4, 3)) - (b1(2, 3) * b1(4, 2))))
n342 = (b1(1, 2) * ((b1(2, 1) * b1(4, 3)) - (b1(2, 3) * b1(4, 1))))
n343 = (b1(1, 3) * ((b1(2, 1) * b1(4, 2)) - (b1(2, 2) * b1(4, 1))))
c(3, 4) = -(n341 - n342 + n343) / nc
n411 = (b1(1, 2) * ((b1(2, 3) * b1(3, 4)) - (b1(2, 4) * b1(3, 3))))
n412 = (b1(1, 3) * ((b1(2, 2) * b1(3, 4)) - (b1(2, 4) * b1(3, 2))))
n413 = (b1(1, 4) * ((b1(2, 2) * b1(3, 3)) - (b1(2, 3) * b1(3, 2))))
c(4, 1) = -(n411 - n412 + n413) / nc
n421 = (b1(1, 1) * ((b1(2, 3) * b1(3, 4)) - (b1(2, 4) * b1(3, 3))))
n422 = (b1(1, 3) * ((b1(2, 1) * b1(3, 4)) - (b1(2, 4) * b1(3, 1))))
n423 = (b1(1, 4) * ((b1(2, 1) * b1(3, 3)) - (b1(2, 3) * b1(3, 1))))
c(4, 2) = (n421 - n422 + n423) / nc
n431 = (b1(1, 1) * ((b1(2, 2) * b1(3, 4)) - (b1(2, 4) * b1(3, 2))))
n432 = (b1(1, 2) * ((b1(2, 1) * b1(3, 4)) - (b1(2, 4) * b1(3, 1))))
n433 = (b1(1, 4) * ((b1(2, 1) * b1(3, 2)) - (b1(2, 2) * b1(3, 1))))
c(4, 3) = -(n431 - n432 + n433) / nc
n441 = (b1(1, 1) * ((b1(2, 2) * b1(3, 3)) - (b1(2, 3) * b1(3, 2))))
n442 = (b1(1, 2) * ((b1(2, 1) * b1(3, 3)) - (b1(2, 3) * b1(3, 1))))
n443 = (b1(1, 3) * ((b1(2, 1) * b1(3, 2)) - (b1(2, 2) * b1(3, 1))))
c(4, 4) = (n441 - n442 + n443) / nc
b3(1, 1) = ((b1(2, 2) * b1(3, 3)) - (b1(2, 3) * b1(3, 2))) / nc
b3(1, 2) = ((b1(2, 1) * b1(3, 3)) - (b1(2, 3) * b1(3, 1))) / -nc
b3(1, 3) = ((b1(2, 1) * b1(3, 2)) - (b1(2, 2) * b1(3, 1))) / nc
b3(2, 1) = ((b1(1, 2) * b1(3, 3)) - (b1(1, 3) * b1(3, 2))) / -nc
b3(2, 2) = ((b1(1, 1) * b1(3, 3)) - (b1(1, 3) * b1(3, 1))) / nc
b3(2, 3) = ((b1(1, 1) * b1(3, 2)) - (b1(1, 2) * b1(3, 1))) / -nc
b3(3, 1) = ((b1(1, 2) * b1(2, 3)) - (b1(1, 3) * b1(2, 2))) / nc
b3(3, 2) = ((b1(1, 1) * b1(2, 3)) - (b1(1, 3) * b1(2, 1))) / -nc
b3(3, 3) = ((b1(1, 1) * b1(2, 2)) - (b1(1, 2) * b1(2, 1))) / nc
l(1, 1) = (c(1, 1) * g1(1, 1)) + (c(1, 2) * g1(1, 2)) + (c(1, 3) * g1(1, 3)) + (c(1, 4) * g1(1, 4))
l(2, 1) = (c(2, 1) * g1(1, 1)) + (c(2, 2) * g1(1, 2)) + (c(2, 3) * g1(1, 3)) + (c(2, 4) * g1(1, 4))
l(3, 1) = (c(3, 1) * g1(1, 1)) + (c(3, 2) * g1(1, 2)) + (c(3, 3) * g1(1, 3)) + (c(3, 4) * g1(1, 4))
l(4, 1) = (c(4, 1) * g1(1, 1)) + (c(4, 2) * g1(1, 2)) + (c(4, 3) * g1(1, 3)) + (c(4, 4) * g1(1, 4))
FOR m = -25 TO 25 STEP .05
COLOR (RND(r) * 10)
O = ((l(4, 1)) * (m ^ 3)) + ((l(3, 1)) * (m ^ 2)) + (l(2, 1) * m) + l(1, 1)
PSET (m, O)
NEXT m
COLOR (7)
PRINT l(4, 1), "x^3+", l(3, 1); "x^2+"; l(2, 1); "x+"; l(1, 1)
PRINT "again", w11, w21, w31, w12, w22, w32, w13, w23, w33
FOR see = 1 TO 1000000
NEXT see
lll = lll + 1
SYSTEM
bbl: