We present the analysis of O-C diagrams for six detached eclipsing binary systems with 3rd and/or 4th body solutions. These solutions were obtained by the fitting of O-C diagrams using genetic algorithms and Monte-Carlo simulation. O-C diagrams were formed using times of minima collected from publications and our observations.