2 Ost
Ну, давайте вот здесь сначала разберемся с мезонами. Вот прога для самого простого случая с ud-мезонами
to4 = 50; pa3 = 100;
s = 0.04168723670021172737215360293164504721700486232504;
ew = 30.1567307802595803785344923081668874541094417852214;
ev = 0.51099891000000;
V = 1;
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
ui := x^(x^(Y^x));
di := x^(x^V);
(*Meson*)
p[1] := di^ui;
p[2] := ui^di;
(*Calculate*)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
x1 = 0.00100000000000000000000000000000000000000000000000000;
x2 = 0.12000000000000000000000000000000000000000000000000000;
b1 = 0.00000100000000000000000000000000000000000000000000000;
b7 =\[InvisibleSpace]0.00010000000000000000000000000000000000000000000000000;
Do[
pb = b7 - b1; g = pb/7;
b2 = b1 + g; b6 = b7 - g;
Y = b1;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax1 = ko;
Y = b1 + g;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax2 = ko;
Y = b1 + 2*g;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax3 = ko;
Y = b7 - 2*g;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax5 = ko;
Y = b7 - g;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax6 = ko;
Y = b7;
ko = 0;
For[i = 1, i < 3,
res = FindMinimum[p[i], {x, x1, x2}];
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ax7 = ko;
If[And[ax1 > ax2, ax2 > ax3], {b1 = b2}];
If[And[ax7 > ax6, ax6 > ax5], {b7 = b6}],
{pa3}];
o = ax3*ev;
Print["======================="]
Print["Pi +- = 139.57018"];
Print["Macca = ", o, " MeV"];
Print["========================"];
Print["V = ", V];
Print["A = ", Y];
Здесь вычисляется минимум массы пи-мезона (+-)
Один из параметров равен единице, а другой контрпараметр имеет строго определённое значение, которое соответствует минимуму массы.