Modelling three dimensional gene regulatory networks

We consider the three-dimensional gene regulatory network (GRN in short). This model consists of ordinary differential equations of a special kind, where the nonlinearity is represented by a sigmoidal function and the linear part is present also. The evolution of GRN is described by the solution vector X(t), depending on time. We describe the changes that system undergoes if the entries of the regulatory matrix are perturbed in some way.