Glass piecewise linear ODE models are frequently used for simulation of neural and gene regulatory networks Efficient computational tools for automatic synthesis of such models are highly desirable However  the existing algorithms for the identification of desired models are limited to four-dimensional networks  and rely on numerical solutions of eigenvalue problems