Продолжим не дожидаясь Меркулова.
Это текст программы.
\!\(\*
RowBox[{\(a = 0.0072973525628`30;\),
"\[IndentingNewLine]", \(mm = 105.6583715`30;\), \( (*MeV*) \), "\
\[IndentingNewLine]", \(mme = \(-9.2847647043`30\)*10^\((\(-24\))\);\), \( (*\
\ J/T\ *) \), "\[IndentingNewLine]", \(mmp = \
1.41060679736`30*10^\((\(-26\))\);\), \( (*\ J/T\ *) \), \
"\[IndentingNewLine]", \(H = 1.1904982`30;\), "\[IndentingNewLine]", \(h8 =
8.4778765`30;\), "\[IndentingNewLine]",
StyleBox[\(Ea = mm\/\(3*\((\((\(a*mme\)\/mmp)\)\^2 - 1)\)\);\),
FontSize->18], "\[IndentingNewLine]", \(m = 1;\), "\[IndentingNewLine]",
StyleBox[\(E1 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((
1 - n + h8\/H)\)\))\);\),
FontSize->16],
"\[IndentingNewLine]", \(m = 2;\), "\[IndentingNewLine]",
StyleBox[\(E2 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((1 - n +
h8\/H)\)\))\);\),
FontSize->
16], "\[IndentingNewLine]", \(m = 3;\), "\[IndentingNewLine]",
StyleBox[\(E3 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((1 -
n + h8\/H)\)\))\);\),
FontSize->16],
"\[IndentingNewLine]", \(m = 4;\), "\[IndentingNewLine]",
StyleBox[\(E4 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((1 -
n + h8\/H)\)\))\);\),
FontSize->16],
"\[IndentingNewLine]", \(m = 5;\), "\[IndentingNewLine]",
StyleBox[\(E5 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((1 - n +
h8\/H)\)\))\);\),
FontSize->16],
"\[IndentingNewLine]", \(m = 6;\), "\[IndentingNewLine]",
StyleBox[\(E6 = 1\/2*\((m + \(Ea*H\^m\)\/h8*\(∏\+\(n = 1\)\%m\((1 -
n + h8\/H)\)\))\);\),
FontSize->16],
StyleBox["\[IndentingNewLine]",
FontSize->16], \( (*\ \ *) \), "\[IndentingNewLine]", \(n = 1;\), "\
\[IndentingNewLine]", \(M1 = \((5 -
2*n)\)*E3 - \((2 + 2*
n)\)*E1;\), "\[IndentingNewLine]", \(Print["\<Ea = \>",
Ea];\), "\[IndentingNewLine]", \(Print["\<E1 = \>", E1];\), "\
\[IndentingNewLine]", \(Print["\<E2 = \>", E2];\), "\[IndentingNewLine]", \
\(Print["\<E3 = \>", E3];\), "\[IndentingNewLine]", \(Print[
"\<E4 = \>", E4];\), "\[IndentingNewLine]", \(Print[
"\<E5 = \>", E5];\), "\[IndentingNewLine]", \(Print["\<E6 = \>",
E6];\), "\[IndentingNewLine]", \(Print["\<\>"];\), "\[IndentingNewLine]", \
\(n = 1;\), "\[IndentingNewLine]", \(M2 = \((5 - 2*n)\)*E5 - M1;\),
"\[IndentingNewLine]", \(Print["\<Muon = \>", 105.6583715];\),
"\[IndentingNewLine]", \(Print[M1];\), "\[IndentingNewLine]", \
\(Print["\<Taon = \>", 1776.82];\),
"\[IndentingNewLine]", \(Print[M2];\), "\[IndentingNewLine]", \(n \
= 1;\), "\[IndentingNewLine]", \(M3 = \((5 - 2*n)\)*E6 - M2;\), "\
\[IndentingNewLine]", \(Print["\<XyempoH = 10 GeV по Баруту\>"];\),
"\[IndentingNewLine]", \(Print[M3];\)}]\)
Я продолжил для таона и для следующего за таоном мифического "хренотрона". (Согласно подъё6kaм CASTRO).
Это заключительная часть читабельного кода:
n = 1;
M1 = (5 - 2*n)*E3 - (2 + 2*n)*E1;
Print["Ea = ", Ea];
Print["E1 = ", E1];
Print["E2 = ", E2];
Print["E3 = ", E3];
Print["E4 = ", E4];
Print["E5 = ", E5];
Print["E6 = ", E6];
Print[""];
n = 1;
M2 = (5 - 2*n)*E5 - M1;
Print["Muon = ", 105.6583715];
Print[M1];
Print["Taon = ", 1776.82];
Print[M2];
n = 1;
M3 = (5 - 2*n)*E6 - M2;
Print["XyempoH = 10 GeV по Баруту"];
Print[M3];
Впереди этого кода стоит вычисление соответственно: Ea,E1,E2,E3,E4,E5 и E6.
Так же здесь зашита подгонка величины H к массе мюона, чтобы вычислить величину h8.
В данном конкретном случае для
H = 1.1904982`30;
соответствует
h8 = 8.4778765`30;
Результат будет такой:
Taon Macca = 1777.82 МэВ Экспериментальное
1840.76 МэВ
Хренотрон по Баруту = 10 ГэВ
3064.76 МэВ