D:\_jknuba_25_26_autem\_jg_ERE_2025\j_Artc_KG_2025.wxmx
--> kill(all);

\[\]\[\tag{%o0} \ensuremath{\mathrm{done}}\]

--> Q:1.5;del:1.6;i:0.007;n:0.014;t:15;C:400;h_s:del·r_s;

\[\]\[\tag{%o1} 1.5\]

\[\]\[\tag{%o2} 1.6\]

\[\]\[\tag{%o3} 0.007\]

\[\]\[\tag{%o4} 0.014\]

\[\]\[\tag{%o5} 15\]

\[\]\[\tag{%o6} 400\]

\[\]\[\tag{%o7} 1.6 {r_s}\]

--> K_0:Q/sqrt(i);

\[\]\[\tag{%o8} 17.928429140015904\]

y = -0,2392x3 + 0,7665x2 - 0,0778x + 0,0073
--> K_s(di):=−0.2392·((di)··3) + 0.7665·((di)··2) − 0.0778·(di) + 0.0073;

\[\]\[\tag{%o9} \mathop{K\_ s}\left( \ensuremath{\mathrm{di}}\right) \mathop{:=}\mathop{-}0.2392 {{\ensuremath{\mathrm{di}}}^{3}}\mathop{+}0.7665 {{\ensuremath{\mathrm{di}}}^{2}}\mathop{+}\mathop{-}0.0778 \ensuremath{\mathrm{di}}\mathop{+}0.0073\]

--> K_s(del);

\[\]\[\tag{%o10} 0.8652968\]

--> d16:[2.694, 0.608];

\[\]\[\tag{%o11} \left[ 2.694\mathop{,}0.608\right] \]

1 - iterations j=4 r1:400 D1:800
--> j:4;

\[\]\[\tag{%o12} 4\]

--> ri:[0.1,0.2,0.3,0.4,0.5,0.6,0.7];

\[\]\[\tag{%o13} \left[ 0.1\mathop{,}0.2\mathop{,}0.3\mathop{,}0.4\mathop{,}0.5\mathop{,}0.6\mathop{,}0.7\right] \]

--> r:ri[j];

\[\]\[\tag{%o14} 0.4\]

--> wi:d16[1]·r··2;

\[\]\[\tag{%o15} 0.4310400000000001\]

--> Ri:d16[2]·r;

\[\]\[\tag{%o16} 0.2432\]

Большаков с.305 Таб.5 R=0,24 n=0.014
--> R24:[37,33.4,30.4,27.5,25.6];

\[\]\[\tag{%o17} \left[ 37\mathop{,}33.4\mathop{,}30.4\mathop{,}27.5\mathop{,}25.6\right] \]

--> W:R24[4]; Wtr1:K_s(del)·W;

\[\]\[\tag{%o18} 27.5\]

\[\]\[\tag{%o19} 23.795662\]

--> Ktr1:Wtr1·wi;

\[\]\[\tag{%o20} 10.256882148480003\]

--> if K_0 > Ktr1 then j:j+2 else j:j−2; display(j);

\[\]\[\tag{%o21} 6 \]\[j\mathop{=}6\]

\[\]\[\tag{%o22} \ensuremath{\mathrm{done}}\]

2 itrations
--> r:ri[j];

\[\]\[\tag{%o23} 0.6\]

--> d16[1];wi:d16[1]·r··2;

\[\]\[\tag{%o24} 2.694\]

\[\]\[\tag{%o25} 0.9698399999999999\]

--> Ri:d16[2]·r;

\[\]\[\tag{%o26} 0.36479999999999996\]

Большаков с.305 Таб.5 R=(0,24)0.36 n=0.014
--> R36:[47.8,43.4,39.6,36.5,33.8];

\[\]\[\tag{%o27} \left[ 47.8\mathop{,}43.4\mathop{,}39.6\mathop{,}36.5\mathop{,}33.8\right] \]

--> W:R36[4]; Wtr1:K_s(del)·W;

\[\]\[\tag{%o28} 36.5\]

\[\]\[\tag{%o29} 31.5833332\]

--> Ktr1:Wtr1·wi;

\[\]\[\tag{%o30} 30.630779870687995\]

--> if K_0 > Ktr1 then j:j+2 else j:j−1; display(j);

\[\]\[\tag{%o31} 5 \]\[j\mathop{=}5\]

\[\]\[\tag{%o32} \ensuremath{\mathrm{done}}\]

3 Iterations
--> r:ri[j];

\[\]\[\tag{%o33} 0.5\]

--> d16[1];wi:d16[1]·r··2;

\[\]\[\tag{%o34} 2.694\]

\[\]\[\tag{%o35} 0.6735\]

--> Ri:d16[2]·r;

\[\]\[\tag{%o36} 0.304\]

Большаков с.305 Таб.5 R=(0,24)(0.36)(0.3) n=0.014
--> R30:[42.6,38.6,35.2,32.3,29.9];

\[\]\[\tag{%o37} \left[ 42.6\mathop{,}38.6\mathop{,}35.2\mathop{,}32.3\mathop{,}29.9\right] \]

--> W:R30[4]; Wtr1:K_s(del)·W;

\[\]\[\tag{%o40} 32.3\]

\[\]\[\tag{%o41} 27.949086639999997\]

--> Ktr1:Wtr1·wi;

\[\]\[\tag{%o42} 18.823709852039997\]

Relations
--> Rel:(K_0 − Ktr1)/(K_0/100);

\[\]\[\tag{%o44} \mathop{-}4.993637228516823\]

--> if Rel < 0 then Rel:Rel·(−1);

\[\]\[\tag{%o45} 4.993637228516823\]

--> if Rel < 7 then display ("stop") else display ("continue");

\[\]\[\mbox{}\\"stop"\mathop{=}"stop"\]

\[\]\[\tag{%o46} \ensuremath{\mathrm{done}}\]

if K_0 > Ktr1 then j:j+2 else j:j-1; display(j);
Exactly

Created with wxMaxima.

The source of this Maxima session can be downloaded here.