5.Каждому элементу вектора характеризующего температуру стенки канала в любой момент времени присваивается постоянное значение температуры Число циклов присвоения значений вектору увеличивают на один, так-так его элементов на один должно быть больше чем число интервалов разбиений.

6.Для вычисления матрицы определяющей температуру цилиндра по радиусу в любой момент времени используем два вложенных цикла For. Во внутреннем цикле предусмотрено изменение радиуса цилиндра, и вычисление температурного поля в заданный момент времени.

7.При переходе к внешнему циклу отсчёт по времени увеличивается на единицу. Значение производной температуры по радиусу в любой момент времени равно нулю и поэтому, чтобы учесть ещё одно краевое условие при переходе от внешнего цикла к внутреннему значение последней температуры копируется два раза.

8.После получения матрицы температур надо построить график. Чтобы координатные оси были проградуированные удобно для использования в матрице температур переставляют столбцы. Осуществляется это с использованием двух вложенных циклов.

9.Далее следует вывод графика и градуировка его осей.