Курсовая работа: Исследование неявного метода Эйлера для линейной системы ОДУ с постоянным и переменным шагом

Жесткая матрица А:

function A=a();

A=[-50 50;

50 -50.1];

Нежесткая матрица А:

function B=b();

B=[5/2; 0];

Матрица U:

function U=u(t);

U=[1];

Переменные

- Выходные переменные:

tout – выходной вектор времени;

yout – выходное значение функций (соответствует вектору времени);

eout – выходное значение ошибок определения функций (высчитано по методу Рунге).

- Внутренние переменные:

A , B , U – соответствующие матрицы (вектора) из формулы Х = АХ+В U ( t ) .

t – текущее время;

y – текущие значения функций;

ym1 – предыдущие значения функций;

y p 1 – следующие значения функций;

y t –значения функций, высчитанное с половинчатым шагом (для определения ошибки по методу Рунге);

h – текущий шаг;

h1 – предыдущий шаг;

n – вектор с единицами с размерностями как y;

I – единичная матрица;

eam – текущая ошибка, вычисляемая для определения допустимого шага;

- Входы:

funA, funB, funU – внешние функции вычисления A, B и U;