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

UptoLike

53
Ниже приведен листинг программы для решения рассматриваемой
задачи (на языке программирования Pascal)
uses crt;
const mf=500;
type
vector=array[1..mf] of real;
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, N, N1, N2 : integer;
T, alfa, beta : vector;
ai, bi, ci, fi : real;
a1, lamda1, ro1, c1 : real;
a2, lamda2, ro2, c2 : real;
h, tau, t_end, time : real;
T0, Tl, Tr, L : real;
f, g : text;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество промежутков в первой части пластины,
N1');
Readln(N1);
Writeln('Введите количество промежутков во
второй част
и пластины,
N2');
Readln(N2);
Writeln('Введите окончание по времени, t_end');
Readln(t_end);
Writeln('Введите толщину пластины, L');
Readln(L);
Writeln('Введите коэффициент теплопроводности материала первой
части пластины, lamda1');
Readln(lamda1);
Writeln('Введите коэффициент теплопроводности материала второй
части пластины, lamda2');
Readln(lamda2);
Writeln('Введите плотность материала первой части пластины, ro1');
Readln(ro1);
Writeln('Введите плотность материала второй части пластины, ro2');
Readln(ro2);
Writeln('Введите теплоемкость материала первой части пластины, c1');