Разностные методы решения задач теплопроводности. Кузнецов Г.В - 69 стр.

UptoLike

69
Ниже приведен листинг программы для решения рассматриваемой
задачи
(на языке программирования Pascal)
uses crt;
const mf=102;
type
vector1=array[1..mf] of real;
vector2=array[1..mf,1..mf] of real;
var
{раздел описания переменных, которые мы будем использовать в
программе}
i, j, Nx, Ny : integer;
T : vector2;
alfa, beta : vector1;
ai, bi, ci, fi : real;
a, lamda, ro, c : real;
hx, hy, tau, t_end, time : real;
T0, L, H, Th, Tc : real;
f, g : text;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество пространственных узлов в пластине по оси
х, Nx');
Readln(Nx);
Writeln('Введите количество пространственных узлов в пластине по оси
y, Ny');
Readln(Ny);
Writeln('Введите окончание по времени, t_end');
Readln(t_end);
Writeln('Введите длину пластины, L');
Readln(L);
Writeln('Введите толщину пластины, H');
Readln(H);
Writeln('Введите коэффициент теплопроводности материала пластины,
lamda');
Readln(lamda);
Writeln('Введите плотность материала пластины, ro');
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите температуру на границе х = 0 области р
ешения, Th');