Eсть одна вещица, которую я просто обязан здесь тиснуть.
Сказать, что она меня удивила — ничего не сказать.
Вот код:
to4 = 30; wv = 28;
A = N[GoldenRatio - 1, 30];
K = 1.0`30*10*9;
x := N[(1 + x0)^(1/(1 - y0)), to4];
qu := N[(x^x)^(A^x), to4];
x1 = 0.2`30;
x2 = 0.4`30;
Plot[qu, {x, x1, x2}, ImageSize -> {300, 300}];
res = FindMinimum[qu, {x, x1, x2}, AccuracyGoal -> Automatic,
PrecisionGoal -> Automatic, WorkingPrecision -> wv];
x0 = x /. Last[res];
y0 = First[res];
Print["x0 = ", x0];
Print["y0 = ", y0];
Print["qu = ", qu];
Print["2.085 GeV"];
Код наипростейший. Опишу словами.
Значит, берётся вот эта самая формула масс, только при s = 0. Далее, берём формулу векторного бозона (x^x)^(A^x), в ней ищем координаты минимума данной функции (x0, y0), результаты вставляем в общую формулу масс. Всё.
Далее, выбираем значение A (это ищем руками). То есть, подгоняем. Пробуем Золотое сечение.
Результат знаете какой получается?
2.08774...
Только без GeV
А в реале, по данным PDG для W-бозона, среднее значение = 2.085 GeV +- 0.042
Далее можно посмотреть и для Z-бозона. Всё ровно то же, только вместо A = 0.618033... , вместо Золотого сечения берём 2/3.
Результат: 2.47832...
B реале, по данным PDG для Z-бозона, среднее значение = 2.4952 GeV +- 0.0023
... Сижу oxyeвший... . А как быть с GeV'ом? GeV-то здесь откуда взялся?
То есть, GeV — это ширина распада W-бозона. Но у меня-то здесь он с какой сырости?
Ню каращё, пусть. Пусть это случайное совпадение. Совпадение 2 раза. С W-бозоном и с Z-бозоном. 2 раза. Пусть. Х с ним.
Теперь возвращаемся к формуле
Только теперь под "М" мы будем подразумевать не массу, а ширину распада. Хехе, но она же и время жизни.
Соответствие между ГэВ и секундами следующее:
1 GeV = 6.582119569 *10^(-25) секунд.
И что мы с этого поимеем?
