A method of formalizing the analysis of asymptotic properties of solutions to systems of differential equations with distributed time-delays and Boolean-type nonlinearities is offered. Such objects arise in many applications, but of most importance are systems coming from gene regulatory networks (GRN). The dynamics of GRN are governed by sigmoid-type nonlinearities which are close to the step functions. This is due to the fact that genes are only activated if certain concentrations are close to the respective threshold values. The delay effects arise from the time required to complete transcription, translation and diffusion to the place of action of a protein.
We describe an algorithm of localizing stationary points in the presence of delays as well as stability analysis around such points. The basic technical tool consists in replacing step functions with the so-called "logoid functions", combined with a special modification of the well-known "linear chain trick", and investigating the smooth systems thus obtained.
A significant part of this framework is based on asymptotic analysis of singularly perturbed matrices, where we apply Mathematica to be able to derive exact stability criteria. This work is a brief review of the results presented in Ref. 1.