Марковский процесс. Структура звезда, с

вершинами, где

. Получаем уравнение Колмогорова:

;

;

;

.
Решение:
Начальные условия

;
Используя нормировочное уравнение

.
Из получившихся 4 уравнении находим финальные вероятности

. С помощью матрицы интенсивностей переходов получается то же самое.
На графике функции, по

задано время

, по

- вероятность.
Нужно посмотреть как будет изменяться 4 графика, для 4 уравнении, со временем

. В маткад, с помощью оператора

- это метод рунге-кутты, нашли значения и построили графики.
Нужно найти время за которое графики доходят до установившееся время

(на графике: значение времени в котором графики пересекают "верх" и "ниж"). А установившееся время для этой структуры,

. Если бы это например была обычная функция типа

, мы вывели

. Или можно было посмотреть по значениям которые вышли после расчета в маткаде, но оно округлено и точное значение

определить трудно.
Результат вычислении в виде графика:

.
Как вычислить время за которое график достигает значения?