restart; assume(t::real, t1::real, t2::real); z[1] := t = conjugate(a1)*a4+conjugate(b1)*b4+conjugate(d1)*d4; z[2] := 0 = conjugate(a5)*a4+conjugate(b5)*b4+conjugate(d5)*d4; z[3] := 0 = conjugate(a1)*a8+conjugate(b1)*b8+conjugate(d1)*d8; z[4] := 0 = conjugate(a5)*a8+conjugate(b5)*b8+conjugate(d5)*d8; z[5] := t = conjugate(a4)*a3+conjugate(b4)*b3+conjugate(f4)*f3; z[6] := 0 = conjugate(a8)*a3+conjugate(b8)*b3+conjugate(f8)*f3; z[7] := 0 = conjugate(a4)*a7+conjugate(b4)*b7+conjugate(f4)*f7; z[8] := 0 = conjugate(a8)*a7+conjugate(b8)*b7+conjugate(f8)*f7; z[9] := t = conjugate(a3)*a2+conjugate(e3)*e2+conjugate(f3)*f2; z[10] := 0 = conjugate(a7)*a2+conjugate(e7)*e2+conjugate(f7)*f2; z[11] := 0 = conjugate(a3)*a6+conjugate(e3)*e6+conjugate(f3)*f6; z[12] := 0 = conjugate(a7)*a6+conjugate(e7)*e6+conjugate(f7)*f6; z[13] := t = conjugate(a2)*a1+conjugate(e2)*e1+conjugate(d2)*d1; z[14] := 0 = conjugate(a6)*a1+conjugate(e6)*e1+conjugate(d6)*d1; z[15] := 0 = conjugate(a2)*a5+conjugate(e2)*e5+conjugate(d2)*d5; z[16] := 0 = conjugate(a6)*a5+conjugate(e6)*e5+conjugate(d6)*d5; z[17] := t1 = conjugate(a3)*a1+conjugate(b3)*b1+conjugate(e3)*e1; z[18] := 0 = conjugate(a7)*a1+conjugate(b7)*b1+conjugate(e7)*e1; z[19] := 0 = conjugate(a3)*a5+conjugate(b3)*b5+conjugate(e3)*e5; z[20] := 0 = conjugate(a7)*a5+conjugate(b7)*b5+conjugate(e7)*e5; z[21] := t2 = conjugate(a4)*a2+conjugate(d4)*d2+conjugate(f4)*f2; z[22] := 0 = conjugate(a8)*a2+conjugate(d8)*d2+conjugate(f8)*f2; z[23] := 0 = conjugate(a4)*a6+conjugate(d4)*d6+conjugate(f4)*f6; z[24] := 0 = conjugate(a8)*a6+conjugate(d8)*d6+conjugate(f8)*f6; " "; "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++"; " "; "Solving of system of equations from Blocks 1 -- 4"; solution1 := eliminate({z[1], z[2], z[3], z[4]}, {a1, a5, b1, b5})[1]; solution1 := factor(simplify(expand(expand(expand(solution1))))); solution2 := eliminate({z[5], z[6], z[7], z[8]}, {a3, a7, b3, b7})[1]; solution2 := factor(simplify(expand(expand(expand(solution2))))); solution3 := eliminate({z[9], z[10], z[11], z[12]}, {a2, a6, e2, e6})[1]; solution3 := factor(simplify(expand(expand(expand(solution3))))); solution3 := factor(simplify(expand(expand(expand(subs(solution2, solution3)))))); solution4 := eliminate({z[13], z[14], z[15], z[16]}, {d2, d6, e1, e5})[1]; solution4 := factor(simplify(expand(expand(expand(solution4))))); solution4 := factor(simplify(expand(expand(expand(subs(solution1, solution3, solution4)))))); " "; "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++"; " "; "Checking of the results"; for i to 16 do factor(simplify(expand(expand(expand(expand(subs(solution1, solution2, solution3, solution4, z[i]))))))); end do; i := 'i'; " "; "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++"; " "; " Result for Block 5 "; for i from 17 to 20 do eq1[i-16] := factor(simplify(expand(expand(expand(expand(subs(solution1, solution2, solution3, solution4, z[i]))))))): end do: |