The reliable and reproducible identification of gene interaction networks represents one of the grand challenges of both modern molecular biology and computational sciences. Approaches based on careful collection of literature data and network topological analysis, applied to unicellular organisms, have proven to offer results applicable to medical therapies. However, when little a priori knowledge is available, other approaches, not relying so strongly on previous literature, must be used. We propose here a novel algorithm (based on ordinary differential equations) able to infer the interactions occurring among genes, starting from gene expression steady state data.