al: mas1);
var, j: integer;, par1, par11, par2: tipR;
begin:=100; i:=1 to N - 1 doj:=1 to M - 1 doh [i, j] <> 0 then (h [i, j +1] <> 0 ) and (h [i +1, j] <> 0) then:=dx * dx * (kh ^ [i, j +1] + kh ^ [i, j]) / (kh ^ [i + 1, j] + kh ^ [i, j]) / dy / dy;
par1:=pi * pi / (2 * N * N * (1 + par));:=pi * pi / (2 * M * M * (1 +1 / par));
if par1
for i:=2 to 4 do al [i]:=1-power (par11, (i - 1) / 3);; skv;
{5} [9,13]:=q5_1/86.4;
{75} [8,5]:=q75_1/86.4;
{10} [5,15]:=q10_1/86.4;
{8} [9,26]:=q8_1/86.4;
{16}[6,29]:=q16_1/86.4;;(kh,M*N*SizeOf(tipR));(b1,M*N*SizeOf(tipR));(pp,M*N*SizeOf(tipR));(dp,M*N*SizeOf(tipR));(V,M*N*SizeOf(tipR));(g1,M*N*SizeOf(tipR));(f1,M*N*SizeOf(tipR));:=0;i:=1 to N doj:=1 to M doh [i, j] <> 0 thenj <= 21 then kh^[i,j]:=k2*h[i,j]kh^[i,j]:=k1*h[i,j];[i,j]:=giv1[i,j];:=ZapN+h[i,j]*dx*dy*kp*al1*Ppl/pat/z(Ppl/1e6)*z(pat/1e6);^[i,j]:=0;[i,j]:=0;;^[i,j]:=p[i,j];[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;;;(al);:=0;Qdob:=0;delta:=0;:=0;:=0;;:=count+1;
Продовження додатка А
if odd (count) then
begiN {obratnoe yporyadochivanie} i:=2 to N - 1 doj:=M - 1 downto 2 do:=0; c:=0; b:=0; f:=0; a:=0; d :=0; R:=0; h [i, j] <> 0 thenh [i, j - 1]=0 then g:=0g:=(dx / dy) * pp ^ [i, j] * (kh ^ [i, j - 1] + kh ^ [i, j]) / mu / z (pp ^ [i, j] / 1e6) / 2; h [i - 1, j]=0 then c :=0c:=(dy / dx) * pp ^ [i, j] * (kh ^ [i - 1, j] + kh ^ [i, j]) / mu / z (pp ^ [i, j]/1e6) / 2; h [i +1, j]=0 then b:=0b:=(dy / dx) * pp ^ [i, j] * (kh ^ [i +1, j] + kh ^ [i, j]) / mu / z (pp ^ [i, j] / 1e6) / 2; h [i, j +1]=0 then f:=0f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);:=2*pat*qm[i,j]-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);:=d-(g*pp^[i,j- 1] + c * pp ^ [i - 1, j] + a * pp ^ [i, j] + b * pp ^ [i +1, j] + f * pp ^ [i, j +1]) ; odd (count) then k:=k +1; k> 4 then k:=1; (i, j, g, c, a, b, f, R, g1 ^ [i - 1, j], f1 ^ [i, j +1], al [k], V ^ [i - 1, j], V ^ [i, j +1], ggg, fff, v1); ^ [i, j]:=v1 ; ^ [i, j]:=ggg; ^ [i, j]:=fff;;;:=true; i:=N - 1 downto 2 doj:=2 to M - 1 doh [i, j] < ;> 0 then ^ [i, j]:=V ^ [i, j]-g1 ^ [i, j] * dp ^ [i, j - 1]-f1 ^ [i, j] * dp ^ [ i +1, j]; abs (dp ^ [i, j])> 1e - 4 then flag:=false; ^ [i, j]:=pp ^ [i, j] + dp ^ [i, j ];; i:=1 to N doj:=1 to M do begin g1 ^ [i, j]:=0; f1 ^ [i, j]:=0 end; i:=2 to N - 1 doj: =2 to M - 1 do:=0; c:=0; b:=0; f:=0; a:=0; d:=0; R:=0; h [i, j] <> ; 0 thenh [i, j - 1]=0 then g:=0g:=(dx / dy) * pp ^ [i, j] * (kh ^ [i, j - 1] + kh ^ [i, j ]) / mu / z (pp ^ [i, j] / 1e6) / 2; h [i - 1, j]=0 then c:=0c:=(dy / dx) * pp ^ [i, j] * (kh ^ [i - 1, j] + kh ^ [i, j]) / mu / z (pp ^ [i, j] / 1e6) / 2; h [i +1, j]=0 then b :=0b:=(dy / dx) * pp ^ [i, j] * (kh ^ [i +1, j] + kh ^ [i, j]) / mu / z (pp ^ [i, j]/1e6) / 2; h [i, j +1]=0 then f:=0f:=(dx/dy)*pp^[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(pp^[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*2*al1*kp*h[i,j]/z(pp^[i,j]/1e6);:=2*pat*qm[i,j]-dx*dy/dt*2*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6);:=d-(g*pp^[i,j- 1] + c * pp ^ [i - 1, j] + a * pp ^ [i, j] + b * pp ^ [i +1, j] + f * pp ^ [i, j +1]) ; odd (count) then k:=k +1; k> 4 then k:=1; (i, j, g, c, a, b, f, R, b1 ^ [i, j - 1], f1 ^ [i - 1, j], al [k], V ^ [i - 1, j], V ^ [i, j - 1], bbb, fff, v1); ^ [i, j]:=v1 ; ^ [i, j]:=bbb; ^ [i, j]:=fff;...
