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

UptoLike

118
Ниже приведен листинг программы для решения рассматриваемой
задачи (на языке программирования Pascal)
uses crt;
const mf=500;
R=8.31;
eps=1e-5;
Patm=101325;
type
vector=array[1..mf] of real;
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, N : integer;
T, Tn, alfa, beta : vector;
ai, bi, ci, fi, d, d1 : real;
lamda, ro, c, h, tau : real;
q, A, qev, P0, M : real;
T0, L, t_end, time : real;
f, g : text;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество пространственных узлов, N');
Readln(N);
Writeln('Введите окончание по времени, t_end');
Readln(t_end);
Writeln('Введите толщину пластины, L');
Readln(L);
Writeln('Введите коэффициент теплопроводн
ости материала пластины,
lamda');
Readln(lamda);
Writeln('Введите плотность материала пластины, ro');
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите коэффициент аккомодации, A');
Readln(A);
Writeln('Введите плотность теплового потока на границе, q');
Readln(q);
Writeln('Введите тепловой эффект фазового перехода, qev');
Readln(qev);