Мкр явная схема

program yavnaya;

const

alpha=22; lamda=0.16; r0=0.021; p=1200; c=1380; n=21; tau=3200;

a=lamda/(p*c);hi=r0/n; hj=hi*hi/(2*a); m=round(tau/hj+1); t0=95; tcp=15;

var

i,j:integer; r:real;

t:array[0..n,0..m] of real;

begin

for i:=0 to n do

t[i,0]:=t0;

for j:=0 to m-1 do

begin

for i:=1 to n-1 do

begin

r:=i*hi;

t[i,j+1]:=a*hj/hi/hi*(t[i+1,j]-2*t[i,j]+t[i-1,j])+a*hj/2/hi/r*(t[i+1,j]-t[i-1,j])+t[i,j];

end;

t[0,j+1]:=t[1,j+1];

t[n,j+1]:=(tcp*alpha*hi/lamda+t[n-1,j+1])/(1+alpha*hi/lamda);

end;

for i:=0 to n do

begin

write (' | ',t[i,0]:2:3,' |');

write (' | ',t[i,round(m/100)]:2:3,' |');

write (' | ',t[i,round(m/2)]:2:3,' |');

writeln (' | ',t[i,m]:2:3,' |');

end;

end.

Мкр неявная схема

program neyavnay;

const

r0=0.021; t0=95; tcp=15; lambda=0.16; ro=1380; c=1200; alfa=22; a=lambda/ro/c;

tau=3200; n=21; m=619; hi=r0/n; hj=tau/m;

var

i,j,tg:integer;

r,X,Y,Z:real;

t:array [0..n,0..m] of Real;

alpha,betta:array[1..30] of real;

begin

for i:=0 to n do

t[i,0]:=t0;

for j:=0 to m-1 do

begin

for i:=1 to n-1 do

begin

r:=i*hi;

X:=a*hj/hi*(1/hi-1/2/r);

Y:=1+2*a*hj/sqr(hi);

Z:=a*hj/hi*(1/hi+1/2/r);

alpha[1]:=1;

betta[1]:=0;

alpha[i+1]:=Z/(Y-X*alpha[i]);

betta[i+1]:=(X*betta[i]+t[i,j])/(Y-X*alpha[i]);

end;

t[n,j+1]:=(tcp*alfa*hi/lambda+betta[n])/((1+alfa*hi/lambda)*(1-alpha[n]/(1+alfa*hi/lambda)));

for i:=n-1 downto 0 do

begin

t[i,j+1]:=alpha[i+1]*t[i+1,j+1]+betta[i+1];

t[0,j+1]:=t[1,j+1];

end;

end;

for i:=0 to n do

begin

write (' | ',t[i,0]:2:3,' |');

write (' | ',t[i,round(m/100)]:2:3,' |');

write (' | ',t[i,round(m/2)]:2:3,' |');

writeln(' | ',t[i,m]:2:3,' |');

end;

end.