Согласно разработанному алгоритму, на первом этапе мы задаеми описываем все необходимые для расчетов переменные.
m - число шагов по времени
n - число шагов по координате
E - модуль упругости материала
ro - плотность материала
L - длина бруска
T - %время моделирования
dx=L/n - шаг по расстоянию
dt=T/m - шаг по времени
cur=dt/dx*(E/ro)^0.5 - число Курранта
Далее задаем начальные и граничные условия для расчета модели.Задаем циклический процесс для расчета. Для организации вычислений с заданным числом повторяющихся циклов используем конструкцию имеющую следующий вид:
forvar=Выражение. Инструкция..... Инструкция end
Результатом решения является матрица значений продольного перемещения на каждом шаге по времени и по расстоянию (рис. 2.2).
Рисунок 2.2. – фрагмент матрицы решений