порядку дозволеного Щодо похідної
В
(13)
Система (13) у більш компактному вігляді запісується у векторній ФОРМІ
(14)
Тут - вектор стовпець невідоміх функцій,
- Вектор функція правих частин.
До векторному діференціального рівняння (14) можна застосуваті ВСІ методи розглянуті Вище в цьом розділі (Завдяк лінійній структурі всех Розглянуто методів). При цьом у формулах наведених Вище УСІ Величини векторні крім змінної х і Кроку h. p align="justify"> Розглянемо задачу Коші для системи двох ОДУ Першого порядку, де рівняння запісані в Розгорнутим вігляді
(15)
формули методу Рунге-Кутта 4-го порядку точності для Вирішення (15) Такі:
В В
(16)
Контроль правільності Вибори Кроку h у разі Використання методу Рунге-Кутта четвертого порядку точності для системи (15) может буті організованій за помощью обчислення на шкірному кроці h параметрів
(17)
Если Величини порядку кількох сотих одініці, то розрахунок Триває з тим же кроком, ЯКЩО больше однієї десятої, то крок слід Зменшити, ЯКЩО ж менше однієї сотої, то крок можна збільшити
2. Розробка програми для ЕОМ
В В В В В В В В
В
В
В
Висновки
У розрахунково-графічному завданні булу досліджена гомогенна ізотермічна Реакція, склад ее математична модель та знайдені кінетічні залежності компонентів Реакції.
Для Виконання поставленої задачі БУВ Обраний з обгрунтуванням метод чісельної реалізації математичної МОДЕЛІ. ЦІМ методом Було Звертаючись - метод Рунге-Кутті. p align="justify"> Для програмної реалізації методу булу Використана система Mathcad у якій на Основі раніше обраної Функції булу Створена програма для Вирішення поставленої задачі.
Список літератури
1. Кравцов А.В., Новиков А.А., Коваль П.І. Методи аналізу хіміко технологічних процесів. - Томськ: вид-во. ТПУ, 1994.-76с.
. Електронний курс з MathCAD
. Чисельні методи. Навчальний посібник
. Про чисельних методах MathCAD 12 керівництво
. Урмаєв А.С. Основі моделювання на АВМ. - М.: Головна ре...