Вставлю свои 5 копеек.
Вот программа для собственно минимума по протону
to4 = 28;
s = 0.0416872367002117273721536029316450472170;
ew = 30.15673078025958037853449230816688745410;
w = 0.0260499593922058169319422986287300000000;
ev = 0.510998910000000000000000000000000000000;
h = 0.8734230000000000000000000000000000000000;
No = 0.798850000000000000000000000000000000000;
N2 = 0.799100000000000000000000000000000000000;
Print["Ecmb CMEщение!"];
t = 0;
o = 0.0734230000000000000000000000000000000000;
r = 0.1015274300000000000000000000000000000000;
Print[(No + N2)/2];
(* _______________________ *)
x := N[((1 + x0) + s)^(1/(1 - y0)), to4];
(* _______________________ *)
x0 = 0; y0 = 0;
ko := 0;
(* Quark *)
di := x^(x^B);
ai := x^(x^(A^x));
ui := x^(x^(Y^x));
(* Proton *)
p[1] := ai^(ui^(1/di));
p[2] := ai^(di^(1/ui));
p[3] := (ui^(1/di))^ai;
p[4] := (di^(1/ui))^ai;
p[5] := (ui^di)^(1/ai);
p[6] := (di^ui)^(1/ai);
p[7] := ai^(1/(ui^di));
p[8] := ai^(1/(di^ui));
p[9] := ai^(1/(ui^(1/di)));
p[10] := ai^(1/(di^(1/ui)));
p[11] := (ui^(1/di))^(1/ai);
p[12] := (di^(1/ui))^(1/ai);
(* End Proton *)
(* ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ *)
kvo = 50;
kv1 = kvo + 1;
cnucok := Array[g, kvo];
m1 = 0.100000000000000000000000000000000000;
m2 = 0.700000000000000000000000000000000000;
dd = (N2 - No)/kvo;
For[j = 1, j < kv1,
ko = 0;
No = No + dd;
A = No + r;
B = No + o;
Y = No + t;
For[i = 1, i < 13,
res = FindMinimum[p[i], {x, m1, m2}, AccuracyGoal -> 28, PrecisionGoal \
-> 22, WorkingPrecision -> to4];
(*res = FindMinimum[p[i], {x, m1, m2}];*)
x0 = x /. Last[res];
y0 = First[res];
k1 = p[i]/ew; ko = ko + k1;
i++];
ko = ko*ev;
g[j] = ko/2;
j++];
ListPlot[cnucok, PlotJoined -> True, PlotStyle -> Hue[.6]];
Print["Pro = 938.272036"];
Этот минимум равен 0.7989750
и соответствует значению массы 802.275 MeV
Можно хихикать, что не совпадает. Посмотрим, как у этих ребят в желтых штанах совпадёт SCC с 4 GeV