a = \[Alpha]*(r + s)/2 /. {r -> 0.00037, s -> 0.0013};
b = \[Alpha]*(\[Alpha]*r*s - 1) /. {r -> 0.00037, s -> 0.0013};
Subscript[\[Beta], 1] = Sqrt[Sqrt[a^2 - b] - a];
Subscript[\[Beta], 2] = Sqrt[Sqrt[a^2 - b] + a];
Subscript[m,
1] = (\[Alpha]*s + Subscript[\[Beta], 1]^2)/Subscript[\[Beta], 1];
Subscript[m,
2] = (\[Alpha]*s - Subscript[\[Beta], 2]^2)/Subscript[\[Beta], 2];
T1 = {{1, 0, 1, 0}, {Subscript[m, 1]*Subscript[\[Beta], 1], 0,
Subscript[m, 2]*Subscript[\[Beta], 2],
0}, {0, (Subscript[\[Beta], 1] - Subscript[m, 1]),
0, (Subscript[\[Beta], 2] + Subscript[m, 2])}, {-Subscript[m, 1]*
Subscript[\[Beta], 1]*t, Subscript[\[Beta],
1], -Subscript[m, 2]*Subscript[\[Beta], 2]*t, Subscript[\[Beta],
2]}} /. {t -> 0.2154};
T2 = Inverse[{{Cosh[Subscript[\[Beta], 1]/2],
Sinh[Subscript[\[Beta], 1]/2], Cos[Subscript[\[Beta], 2]/2],
Sin[Subscript[\[Beta], 2]/2]}, {Subscript[m, 1]*Subscript[\[Beta],
1]*Cosh[Subscript[\[Beta], 1]/2],
Subscript[m, 1]*Subscript[\[Beta], 1]*
Sinh[Subscript[\[Beta], 1]/2],
Subscript[m, 2]*Subscript[\[Beta], 2]*
Cos[Subscript[\[Beta], 2]/2],
Subscript[m, 2]*Subscript[\[Beta], 2]*
Sin[Subscript[\[Beta], 2]/2]}, {(Subscript[\[Beta], 1] -
Subscript[m, 1])*
Sinh[Subscript[\[Beta], 1]/2], (Subscript[\[Beta], 1] -
Subscript[m, 1])*
Cosh[Subscript[\[Beta], 1]/
2], -(Subscript[\[Beta], 2] + Subscript[m, 2])*
Sin[Subscript[\[Beta], 2]/2], (Subscript[\[Beta], 2] + Subscript[
m, 2])*Cos[Subscript[\[Beta], 2]/2]}, {Subscript[\[Beta], 1]*
Sinh[Subscript[\[Beta], 1]/2],
Subscript[\[Beta], 1]*
Cosh[Subscript[\[Beta], 1]/2], -Subscript[\[Beta], 2]*
Sin[Subscript[\[Beta], 2]
/2], Subscript[\[Beta], 2]*Cos[Subscript[\[Beta], 2]/2]}}]
