[Swark]

Может кто-то хочет в MatLab повторить расчеты. Вот текст программки для этого:

Mp = 1.6726e-27; % kg
Ms = 1.989e+30; % kg
h = 6.626e-34; % J*s
c = 2.99792e+8; % m/s
Me = 2*h/c/c;

K1 = (Mp/Me)^(1/7);
K2 = (Ms/16/Mp)^(1/17);
Rm = 4.0e+8; % m
Rp = 0.8418e-15; % m
Rs = 1.39e+9; % m
K3 = ((Rs/Rp)^3/16)^(1/17);
K4 = K3^(1/3);
K5 = Rp*K4^6;
RH = 0.529e-10; % m
K6 = K5*RH/Rp;
K7 = Rs/(16^(1/3))/K4;
MR = 25.17;
MM = 1943.6;
MD = 8.20;
KK = RH / Rp;
Dp = 2.04e+17; % kg/m^3

for j=1:36
i = j - 8;
Ra(j) = Rp*MR^(i-1);
Ma(j) = Mp*MM^(i-1);
De(j) = Dp*MD^(1-i);
RAT(j) = Ra(j)*KK;
DeAT(j) = Ma(j)/(4/3)/3.1415/RAT(j)^3;

end