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 ] ]));