4.4 Блок-схема алгоритма решения задачи

Рисунок 4.1 – Блок-схема к заданию №4

4.5 Текст программы

program imr4;

var a0, a1, a2, h1, p0, h, c, k1, k2, k3, k4, p:real;

i:integer;

f1,f2:text;

begin

assign(f1,'dano.txt');

assign(f2,'otvet.txt');

reset(f1);

rewrite(f2);

read(f1,a0,a1,a2,h1,p0,h,c);

writeln(f2,'p0= ',p0:6:3);

for i:=1 to 10 do

begin

k1:=h*c*(a0*sqr(p0)+a1*p0+a2);

k2:=h*c*(a0*sqr(p0+k1/2)+a1*(p0+k1/2)+a2);

k3:=h*c*(a0*sqr(p0+k2/2)+a1*(p0+k2/2)+a2);

k4:=h*c*(a0*sqr(p0+k3)+a1*(p0+k3)+a2);

p:=p0+(k1+2*k2+2*k3+k4)/6;

writeln(f2,'p',i:1,'= ',p:6:3);

p0:=p;

end;

close(f1);

close(f2);

end.

4.6 Результаты вычисления

p₀ = 60.000 МПа

p₁ = 60.369 МПа

p₂ = 60,740 МПа

p₃ = 61,113 МПа

p₄ = 61,490 МПа

p₅ = 61,868 МПа

p₆ = 62,250 МПа

p₇ = 62,634 МПа

p₈ = 63,020 МПа

p₉ = 63,410 МПа

p₁₀ = 63,801 МПа

4.7 График

Рисунок 4.2 - График зависимости давления газовой скважины от глубины

26