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
LOCATE 14, 1
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 = 1 TO 9 STEP .07
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
LOCATE 17, 1
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)