                                                                                  СтацРешRND
****ОПРЕДЕЛЕНИЕ СТАЦИОНАРНЫХ РЕШЕНИЙ МОДЕЛЕЙ ГГС МЕТОДОМ ГОМОТОПИИ
    с использованием датчика случайных чисел для задания  компонент
    x0(i), i=1,2,...,n, стартового решения в методе гомотопии.

Пусть требуется найти стационарные решения модели, представленной автономной системой вида:

i=1,2,...,n,

dx(i)/dt = ALFA/[1 + BETA*z(i)] - x(i),                                  (1)

где z(i) - функция, представляющая регуляторные связи (первого рода), так что стационарное
решение определяется из системы нелинейных уравнений:

f(i) = ALFA/[1 + BETA*z(i)] - x(i) = 0,  i=1,2,...,n.                    (2)

Рассмотрим систему нелинейных уравнений с параметром q вида:

F(i,q) = ALFA/[1 + BETA*r(i)*z(i)] - x(i) = 0,  i=1,2,...,n,             (3)                 
                                                                       
r(i) = q + (1-q)*[(ALFA/x0(i) - 1)/BETA/z0(i), 0<=q<=1,

где z0(i) - значения z(i), вычисленные при x(i) = x0(i). При этом система (3) при q = 0
имеет точное решение x(i) = x0(i), i=1,2,...,n, а при q = 1  определяет решение системы 
(2), которое, таким образом, находится методом продолжения по параметру q. Правые части
(1) задаются с помощью матрицы связей.

Значения компонент x0(i), i=1,2,...,n, принадлежащие гиперкубу с ребром ALFA, задаются с 
помощью датчика случайных чисел. Вычисление всех стационарных решений системы (2) требует
достаточно большого числа "бросаний".


