X_b := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, t, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 2*t, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_a := [ [ 1, t, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, t, 2*t^2, 0, 0 ], [ 0, 0, 0, 1, t, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 2*t ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_ab := [ [ 1, 0, 2*t, 0, 0, 0, 0 ], [ 0, 1, 0, t, 0, 2*t^2, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, t, 0 ], [ 0, 0, 0, 0, 1, 0, t ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_2ab := [ [ 1, 0, 0, 2*t, 0, 0, 2*t^2 ], [ 0, 1, 0, 0, 2*t, 0, 0 ], [ 0, 0, 1, 0, 0, t, 0 ], [ 0, 0, 0, 1, 0, 0, 2*t ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_3ab := [ [ 1, 0, 0, 0, 2*t, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, t ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_3a2b := [ [ 1, 0, 0, 0, 0, 2*t, 0 ], [ 0, 1, 0, 0, 0, 0, t ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_mb := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, t, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 2*t, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ]; X_ma := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ t, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 2*t, 1, 0, 0, 0 ], [ 0, 0, 2*t^2, 2*t, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 2*t, 1 ] ]; X_mab := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 2*t, 0, 1, 0, 0, 0, 0 ], [ 0, 2*t, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 2*t^2, 0, 2*t, 0, 1, 0 ], [ 0, 0, 0, 0, t, 0, 1 ] ]; X_m2ab := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ t, 0, 0, 1, 0, 0, 0 ], [ 0, 2*t, 0, 0, 1, 0, 0 ], [ 0, 0, t, 0, 0, 1, 0 ], [ 2*t^2, 0, 0, t, 0, 0, 1 ] ]; X_m3ab := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 2*t, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, t, 0, 0, 0, 1 ] ]; X_m3a2b := [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 2*t, 0, 0, 0, 0, 1, 0 ], [ 0, t, 0, 0, 0, 0, 1 ] ]; xS:=simplify(subs({tt = t^(3^m), vt = v^(3^m), ut = u^(3^m)}, [ [ 1, tt, 2*ut, 2*vt+tt*ut, 2*u+2*tt*vt+2*t*tt^3, 2*v+2*ut*vt+2*tt*ut^2+2*t^2*tt^3, tt*v+2*u*ut+2*vt^2+2*tt*ut*vt+2*t*tt^3*ut+t^2*tt^4 ], [ 0, 1, t, ut+t*tt, 2*vt+2*t*tt^2, t*vt+2*ut^2+t*tt*ut, v+t*u+2*ut*vt+2*t*tt*vt+2*t*tt^2*ut+2*t^2*tt^3 ], [ 0, 0, 1, tt, 2*tt^2, vt+tt*ut, u+2*tt*vt+2*tt^2*ut+t*tt^3 ], [ 0, 0, 0, 1, tt, ut, 2*vt+tt*ut ], [ 0, 0, 0, 0, 1, 2*t, ut+t*tt ], [ 0, 0, 0, 0, 0, 1, 2*tt ], [ 0, 0, 0, 0, 0, 0, 1 ] ])); xSp:=simplify(subs({tt = t^(3^m), vt = v^(3^m), ut = u^(3^m)}, [ [ 1, 0, 0, 0, 0, 0, 0 ], [ tt, 1, 0, 0, 0, 0, 0 ], [ 2*ut+t*tt, t, 1, 0, 0, 0, 0 ], [ vt+tt*ut, 2*ut, 2*tt, 1, 0, 0, 0 ], [ 2*u+2*tt*vt+tt^2*ut+t*tt^3, 2*vt+tt*ut, 2*tt^2, 2*tt, 1, 0, 0 ], [ 2*v+t*u+2*ut*vt+t*tt*vt+2*t*tt^2*ut+t^2*tt^3, t*vt+2*ut^2+2*t*tt*ut, vt+t*tt^2, 2*ut+t*tt, 2*t, 1, 0 ], [ tt*v+2*u*ut+2*vt^2+tt*ut*vt+t*tt^3*ut+t^2*tt^4, v+2*ut*vt+tt*ut^2+t^2*tt^3, u+2*tt*vt+2*t*tt^3, vt+tt*ut, ut, 2*tt, 1 ] ])); n_b := simplify(subs(tt = t^(3^m), [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 0, t, 0, 0, 0, 0 ], [ 0, 2*t^(-1), 0, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 2*t, 0 ], [ 0, 0, 0, 0, t^(-1), 0, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ])); n_a := simplify(subs(tt = t^(3^m), [ [ 0, t, 0, 0, 0, 0, 0 ], [ 2*t^(-1), 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 2*t^2, 0, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 0, 2*t^(-2), 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, 2*t ], [ 0, 0, 0, 0, 0, t^(-1), 0 ] ])); n_ab := simplify(subs(tt = t^(3^m), [ [ 0, 0, 2*t, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 2*t^2, 0 ], [ t^(-1), 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, t ], [ 0, 2*t^(-2), 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 2*t^(-1), 0, 0 ] ])); n_2ab := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, 0, 0, 2*t^2 ], [ 0, 0, 0, 0, 2*t, 0, 0 ], [ 0, 0, 0, 0, 0, t, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, t^(-1), 0, 0, 0, 0, 0 ], [ 0, 0, 2*t^(-1), 0, 0, 0, 0 ], [ 2*t^(-2), 0, 0, 0, 0, 0, 0 ] ])); n_3ab := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, 2*t, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, t ], [ 0, 0, 0, 1, 0, 0, 0 ], [ t^(-1), 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 2*t^(-1), 0, 0, 0, 0 ] ])); n_3a2b := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, 0, 2*t, 0 ], [ 0, 0, 0, 0, 0, 0, t ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ t^(-1), 0, 0, 0, 0, 0, 0 ], [ 0, 2*t^(-1), 0, 0, 0, 0, 0 ] ])); n_mb := simplify(subs(tt = t^(3^m), [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 0, 2*t^(-1), 0, 0, 0, 0 ], [ 0, t, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 0, t^(-1), 0 ], [ 0, 0, 0, 0, 2*t, 0, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ])); n_ma := simplify(subs(tt = t^(3^m),[ [ 0, 2*t^(-1), 0, 0, 0, 0, 0 ], [ t, 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 2*t^(-2), 0, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 0, 2*t^2, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, t^(-1) ], [ 0, 0, 0, 0, 0, 2*t, 0 ] ])); n_mab := simplify(subs(tt = t^(3^m), [ [ 0, 0, t^(-1), 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 2*t^(-2), 0 ], [ 2*t, 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, 2*t^(-1) ], [ 0, 2*t^2, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, t, 0, 0 ] ])); n_m2ab := simplify(subs(tt = t^(3^m),[ [ 0, 0, 0, 0, 0, 0, 2*t^(-2) ], [ 0, 0, 0, 0, t^(-1), 0, 0 ], [ 0, 0, 0, 0, 0, 2*t^(-1), 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 2*t, 0, 0, 0, 0, 0 ], [ 0, 0, t, 0, 0, 0, 0 ], [ 2*t^2, 0, 0, 0, 0, 0, 0 ] ])); n_m3ab := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, t^(-1), 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 0, 2*t^(-1) ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 2*t, 0, 0, 0, 0, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, t, 0, 0, 0, 0 ] ])); n_m3a2b := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, 0, t^(-1), 0 ], [ 0, 0, 0, 0, 0, 0, 2*t^(-1) ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 2*t, 0, 0, 0, 0, 0, 0 ], [ 0, t, 0, 0, 0, 0, 0 ] ])); h_b := simplify(subs(tt = t^(3^m), [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, -2*t, 0, 0, 0, 0, 0 ], [ 0, 0, -2*t^(-1), 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, -2*t, 0, 0 ], [ 0, 0, 0, 0, 0, -2*t^(-1), 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ])); h_a := simplify(subs(tt = t^(3^m), [ [ -2*t, 0, 0, 0, 0, 0, 0 ], [ 0, -2*t^(-1), 0, 0, 0, 0, 0 ], [ 0, 0, 4*t^2, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, 4*t^(-2), 0, 0 ], [ 0, 0, 0, 0, 0, -2*t, 0 ], [ 0, 0, 0, 0, 0, 0, -2*t^(-1) ] ])); h_ab := simplify(subs(tt = t^(3^m), [ [ -2*t, 0, 0, 0, 0, 0, 0 ], [ 0, 4*t^2, 0, 0, 0, 0, 0 ], [ 0, 0, -2/t, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, -2*t, 0, 0 ], [ 0, 0, 0, 0, 0, 4/t^2, 0 ], [ 0, 0, 0, 0, 0, 0, -2/t ] ])); h_2ab := simplify(subs(tt = t^(3^m), [ [ 4*t^2, 0, 0, 0, 0, 0, 0 ], [ 0, -2*t, 0, 0, 0, 0, 0 ], [ 0, 0, -2*t, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, -2/t, 0, 0 ], [ 0, 0, 0, 0, 0, -2/t, 0 ], [ 0, 0, 0, 0, 0, 0, 4/t^2 ] ])); h_3ab := simplify(subs(tt = t^(3^m), [ [ -2*t, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, -2*t, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, -2/t, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, -2/t ] ])); h_3a2b := simplify(subs(tt = t^(3^m), [ [ -2*t, 0, 0, 0, 0, 0, 0 ], [ 0, -2*t, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, -2/t, 0 ], [ 0, 0, 0, 0, 0, 0, -2/t ] ])); h_mb := simplify(subs(tt = t^(3^m), [ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, -2/t, 0, 0, 0, 0, 0 ], [ 0, 0, -2*t, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, -2/t, 0, 0 ], [ 0, 0, 0, 0, 0, -2*t, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ])); h_ma := simplify(subs(tt = t^(3^m), [ [ -2/t, 0, 0, 0, 0, 0, 0 ], [ 0, -2*t, 0, 0, 0, 0, 0 ], [ 0, 0, 4/t^2, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, 4*t^2, 0, 0 ], [ 0, 0, 0, 0, 0, -2/t, 0 ], [ 0, 0, 0, 0, 0, 0, -2*t ] ])); h_mab := simplify(subs(tt = t^(3^m), [ [ -2/t, 0, 0, 0, 0, 0, 0 ], [ 0, 4/t^2, 0, 0, 0, 0, 0 ], [ 0, 0, -2*t, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, -2/t, 0, 0 ], [ 0, 0, 0, 0, 0, 4*t^2, 0 ], [ 0, 0, 0, 0, 0, 0, -2*t ] ])); h_m2ab := simplify(subs(tt = t^(3^m), [ [ 4/t^2, 0, 0, 0, 0, 0, 0 ], [ 0, -2/t, 0, 0, 0, 0, 0 ], [ 0, 0, -2/t, 0, 0, 0, 0 ], [ 0, 0, 0, 4, 0, 0, 0 ], [ 0, 0, 0, 0, -2*t, 0, 0 ], [ 0, 0, 0, 0, 0, -2*t, 0 ], [ 0, 0, 0, 0, 0, 0, 4*t^2 ] ])); h_m3ab := simplify(subs(tt = t^(3^m), [ [ -2/t, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, -2/t, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, -2*t, 0, 0 ], [ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, -2*t ] ])); h_m3a2b := simplify(subs(tt = t^(3^m), [ [ -2/t, 0, 0, 0, 0, 0, 0 ], [ 0, -2/t, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ], [ 0, 0, 0, 0, 0, -2*t, 0 ], [ 0, 0, 0, 0, 0, 0, -2*t ] ])); h := simplify(subs(tt = t^(3^m), [ [ tt, 0, 0, 0, 0, 0, 0 ], [ 0, t/tt, 0, 0, 0, 0, 0 ], [ 0, 0, tt^2/t, 0, 0, 0, 0 ], [ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, t/tt^2, 0, 0 ], [ 0, 0, 0, 0, 0, tt/t, 0 ], [ 0, 0, 0, 0, 0, 0, 1/tt ] ])); n := simplify(subs(tt = t^(3^m), [ [ 0, 0, 0, 0, 0, 0, 2 ], [ 0, 0, 0, 0, 0, 2, 0 ], [ 0, 0, 0, 0, 2, 0, 0 ], [ 0, 0, 0, 2, 0, 0, 0 ], [ 0, 0, 2, 0, 0, 0, 0 ], [ 0, 2, 0, 0, 0, 0, 0 ], [ 2, 0, 0, 0, 0, 0, 0 ] ]));