var T: array0..20 of real;
A: array0..20 of real;
t0: real = 10; {начальная температура}
d: real = 0.005; {расстояние между центрами}
r: real = 0.002; {ширина ребра}
c: real = 0.2; {длина |||}
e: real = 0.9; {эпсилон}
s: real = 5.67e-8; {сигма}
h: real = 200; {лямбда}
p: real = 2700; {плотность}
m: real = p * r * d * c; {масса}
ce: real = 900; {теплоемкость}
f: real = r * c; {площадь соприкосновения}
fp: real = c * d; {верхняя площадь}
o: real = 0.01 / (ce * m);
begin
h := h * f / d;
for i := 0 to 20 do
begin
Ti := 0;
Ai := 0;
end;
for i := 1 to 200 do
begin
A1 := T1 + o * ((2 * T2 + (t0 + 273) - 3 * T1) * (h / 2) - e * fp * s * Power(T1 + 273, 4));
writeln('A1: ', A1, ' ');
for j := 2 to 19 do
begin
Aj := Tj + o * ((Tj+1 + Tj-1 - 2 * Tj) * h - e * s * fp * Power(Tj + 273, 4));
writeln('A', j, ': ', Aj, ' ');
end;
A20 := T20 + o * ((2 * T19 + (t0 + 273) - 3 * T20) * (h / 2) - e * fp * s * Power(T20 + 273, 4));
writeln('A20: ', A20, ' ');
for j := 1 to 20 do
begin
Tj := Aj;
end;
writeln(' ');
end;
end.
A= [0]*21
t0=10 #начальная температура
d=0.005 # расстояние между центрами
r=0.002 #ширина ребра
c=0.2 #длина |||
e=0.9 #эпсилон
s=5.67*(10) **(-8)#сигма
h=200 #лямбда
p=2700 #плотность
m=p*r*d*c #масса
ce=900 #теплоемкость
f=r*c #площать соприкосновения
fp=c*d #верхняя площадь
o=0.01/(ce*m)
h=(h*f)/d
for _ in range(1,200):
A[1]=T[1]+o*((2*T[2]+(t0+273)-3*T[1])*(h/(2))-e*fp*s*((T[1]+273) **(4)))
print('A[1]',A[1],' ')
for i in range (2,19):
A[i]=T[i]+o*((T[i+1]+T[i-1]-2*T[i])*h-e*s*fp*((T[i]+273) **(4)))
print('a[i]',A[i],' ')
A[20]=T[20]+o*((2*T[19]+(t0+273)-3*T[20])*(h/2)-e*fp*s*((T[20]+273) **(4)))
print('a[20]',A[20],' ')
for i in range(1,20):
T[i]=A[i]
print (' ')