ПОИСК Статьи Рисунки Таблицы Метод Эйлера для систем дифференциальных уравнений из "Компьютеры Применение в химии" Е — фермент 8 — субстрат Е8 — фермент-субстратный комплекс Р — продукт. [c.230] что итерационные формулы для одного дифференциального уравнения и системы дифференциальных уравнений аналогичны. Однако вторые слагаемые в формулах для системы уравнений if и/2) зависят как от л , так и отз, и2 (впрочем, могут и не зависеть). [c.231] В подпрограмме 10000, вычисляюшей правые части дифференциальных уравнений, используются индексированные переменные D(I) и У(1). В этой подпрограмме для данной конкретной задачи вычисляются две правые части двух дифференциальных уравнений и полученные значения присваиваются переменным D(l) и D(2) соответственно. [c.233] В основной программе после описания массивов переменной N присваивается значение 2 и тем самым задается размерность системы дифференциальных уравнений. Вместо одного начального условия теперь вводится N начальных условий, по одному для каждого дифференциального уравнения системы. Это осуществляется в строках 300—600. Еще одно отличие в основной программе касается вывода данных, поскольку для системы уравнений надо выводить не одно, а N значений (строки 1400—1440). [c.233] Чтобы решить систему дифференциальных уравнений, отличающуюся от системы, решенной в приведенном вместе с программой примере, надо в строке 200 задать другое значение N и в подпрограмме 10000 элементам массива D( ) присвоить арифметические выражения, соответствующие правым частям дифференциальных уравнений. [c.233] Задание 130. В программе СИСТ-ЭЙЛЕР , так же как и в программе ДУ-ЭЙЛЕР1 , не предусмотрено прекращения итерационной процедуры по критерию сходимости рещения. Дополните программу СИСТ-ЭЙЛЕР проверкой результата на устойчивость и измените вывод данных так, чтобы промежуточные результаты выводились на экран только по желанию пользователя. [c.234] Величина Н пропорциональна энтальпии реакции D — коэффициент теплопередачи через стенку реактора. [c.235] Составьте программу для рещения соответствующей системы дифференциальных уравнений и найдите зависимость с от времени при заданных значениях параметров 7 , АГр, E ,H,D и а при начальном условии (/ = 0) с = Сц, Г = Т . (Обратите внимание, что в случае сильноэкзотермической реакции и малого отвода тепла при вьшолнении программы могут встретиться трудности, связанные с приближенным характером численного метода, например отрицательные концентрации в промежуточных результатах, неустойчивость решения, неправдоподобно большие значения температуры и многое другое.) При достаточно малом шаге интегрирования, т. е. если отрезок, ограниченный начальным и конечным значениями аргумента, разделить на очень большое число частичных отрезков, эти трудности устраняются и ЭВМ выдает вполне разумные результаты. [c.235] Вернуться к основной статье