Уравнение теплопроводности (перевести код)

uses crt; const mf=500; type
 
vector=array[1..mf] of real;
 
var {раздел описания переменных, которые мы будем использовать в программе}
 
i, j, N
 
: integer;
 
T, alfa, beta
 
: vector;
 
ai, bi, ci, fi
 
: real;
 
lamda, ro, c, h, tau
 
: real;
 
T0, Th, R, t_end, time
 
: real;
 
f, g
 
: text;
 
begin
 
 
 
clrscr;
 
 
 
{с клавиатуры вводим все необходимые входные параметры}
 
Writeln('Введите количество пространственных узлов, N'); Readln(N);
 
Writeln('Введите окончание по времени, t_end'); Readln(t_end);
 
Writeln('Введите радиус цилиндра, R');
 
Readln(R);
 
Writeln('Введите коэффициент теплопроводности материала пластины, lamda');
 
Readln(lamda);
 
Writeln('Введите плотность материала пластины, ro');
 
Readln(ro);
 
Writeln('Введите теплоемкость материала пластины, c'); Readln(c);
 
Writeln('Введите начальную температуру, T0');
 
Readln(T0);
 
Writeln('Введите температуру на границе r=R, Th');
 
Readln(Th);
 
{определяем расчетный шаг сетки по пространственной координате} h:=R/(N-1);
 
{определяем расчетный шаг сетки по времени} tau:=t_end/100.0;
 
{определяем поле температуры в начальный момент времени} for i:= 1 to N do
 
28
 
T[i]:=T0;
 
{проводим интегрирование нестационарного уравнения теплопроводности}
 
time:=0;
 
while time<t_end do {используем цикл с предусловием} begin
 
{увеличиваем переменную времени на шаг τ} time:=time+tau;
 
{определяем начальные прогоночные коэффициенты на основе левого граничного условия}
 
alfa[1]:=1.0;
 
beta[1]:=0.0;
 
{цикл с параметром для определения прогоночных коэффициентов по формуле (8)}
 
for i:= 2 to N-1do begin
 
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с трехдиагональной матрицей}
 
ai:=0.5*lamda*(2*i-1)/(sqr(h)*(i-1));ci:=0.5*lamda*(2*i-3)/(sqr(h)*(i-1));bi:=ai+ci+ro*c/tau;fi:=-ro*c*T[i]/tau;
 
{alfa[i], beta[i] – прогоночные коэффициенты} alfa[i]:=ai/(bi-ci*alfa[i-1]); beta[i]:=(ci*beta[i-1]-fi)/(bi-ci*alfa[i-1]);
 
end;
 
{определяем значение температуры на правой границе}
 
T[N]:=Th;
 
{используя соотношение (7) определяем неизвестное поле температуры}
 
for i:= N-1downto 1 do T[i]:=alfa[i]*T[i+1]+beta[i];
 
end;
 
{выводим результат в файл}
 
Assign(f,'res.txt');
 
Rewrite(f);
 
Writeln(f,'Радиус цилиндра R = ',R:6:4);
 
Writeln(f,'Число узлов по координате N = ',N);
 
Writeln(f,'Коэффициент теплопроводности материала пластины lamda =
 
',lamda:6:4);
 
Writeln(f,'Плотность материала пластины ro = ',ro:6:4);
 
29
 
 
Writeln(f,'Теплоемкость материала пластины с = ',c:6:4);
 
Writeln(f,'Начальная температура T0 = ',T0:6:4);
 
Writeln(f,'Температура на границе r = R, Th = ',Th:6:4); Writeln(f,'Результат получен с шагом по координате h = ',h:6:4); Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4); Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4); close(f);
 
Assign(g,'tempr.txt');
 
Rewrite(g);
 
for i:=1 to N do
 
writeln(g,' ',h*(i-1):6:3,'',T[i]:8:5); close(g);
 
end.