Решение инженерно-технических задач в среде Mathcad
    Математический пакет Mathcad [23, 42, 44] позволяет специалистам в конкретной научно-технической области быстро освоить работу на компьютере и реализовать на них математические модели, не вдаваясь в тонкости программирования. Mathcad - это зарегистрированная торговая марка фирмы MathSoft, Inc. (http://www.mathsoft.com), представляющая собой сокращение английской фразы: Mathematical Computer Aided Design - математическое проектирование с помощью компьютера.
    Ниже приведены три примера решения типовых задач теплоэнергетики [68], которые дают представления об особенностях работы в среде Mathcad.
    Задача 1. Расчет термодинамических свойств воды и водяного пара. Ниже представлено формирование функции пользователя vw(p), предназначенной для определения удельного объема кипящей воды на линии насыщения (vw) от давления (p), и функции hss(T, p), предназначенной для определения удельной энтальпии перегретого пара (hss) от температуры (T) и давления (p).
    Расчет удельного объема кипящей воды на линии насыщения сводится к интерполяции (здесь задействована встроенная Mathcad-функция interp) кубическими сплайнами (cspline) табличных данных, хранящихся в двух векторах: P - табличные значения давления и V - табличные значения удельного объема. Элементы векторов - величины размерные: давление измеряется в мегапаскалях (106Pa), а удельный объем - в литрах на килограмм (m3/1000kg). В векторах P и V (они транспонируются из векторов-строк в вектора-столбцы) показаны только крайние элементы (далее будут показаны экраны дисплея, отображающие работу в среде Mathcad):
    Интерполяция сплайном для определения удельного объема кипящей воды ведется в диапазоне давлений от 611 Па до 22,088 МПа. Функция error отслеживает значение аргумента p, прерывает выполнение вычислительной процедуры и выдает пользовательское сообщение об ошибке ("Давление низкое" или "Давление высокое"), если значение аргумента не попадает в оговоренный диапазон (эту работу выполняет оператор if).
    Созданную функцию с одним аргументом можно вызвать в традиционном виде (функция), а также в виде префиксного и постфиксного операторов.
    При этом в расчетах возможны различные размерности давления, объема и массы (здесь kgf - килограмм силы, galUK - британский галлон, lb - фунт и lbf - фунт силы и ft - фут). Далее представлено решение обратной задачи через встроенную Mathcad-функцию root (поиск корня уравнения) - определение давления при заданном значении удельного объема кипящей воды.
    Расчет удельной энтальпии перегретого пара ведется также через кубический сплайн, но уже двухмерной интерполяцией. Табличные значения энтропии (килоджоули на килограмм) заносятся в матрицу M, "боковик" которой (без первого элемента) значения температуры (градусы Цельсия), а "шапка" - давления (мегапаскали). Далее эти значения с помощью встроенной функции submatrix разносятся по матрицам TP (у нее два столбца и 18 строк) и M (квадратная матрица с порядком 18). Так достигается компромисс между требованиями наглядного представления табличных данных и спецификой двумерной сплан-интерполяции в среде Mathcad.
    К пользовательским сообщениям об ошибках ввода неверных значений аргументов (температура и давление пара) добавлено новое - "Ошибка: измените p и/или T". Ошибка возникает в том случае, если рассчитанное значение энтальпии окажется меньше наименьшего табличного значения (1930,8 кДж/кг), что может случиться на границах матрицы H, где вместо отсутствующих термодинамических табличных данных записано: -1.
    Пользовательскую функцию hss(T, p) с двумя аргументами можно вызывать в традиционной форме (функция), а также инфиксным и древовидным оператором.
    Аналогичным образом построены и другие функции термодинамических свойств воды и водяного пара: удельная энтальпия кипящей воды на линии насыщения hw(p), удельная энтальпия сухого пара на линии насыщения hs(p), удельная энтропия кипящей воды на линии насыщения sw(p), удельная энтропия сухого пара на линии насыщения ss(p), удельная энтропия перегретого пара sss(T, p) и др. ............