Based on the regressive character of chaotic motion in nonlinear dynamic systems, a numerical regression algorithm is developed, which can be used to research the dark-lines passing through chaotic regions in bifurcation plots. The dark-lines of the parabolic mapping are obtained by using the numerical regression algorithm, and compared with those that are accurately acquired through dark-line equations. Thus the validity of this algorithm is proved. Furthermore,for the Brussel oscillation system and the piecewise linear dynamic system of a gear pair, the dark-lines are researched by using the regression algorithm. By researching the dark-lines in the bifurcation plots of nonlinear dynamic systems,the periodic windows embedded in chaotic regions can be ascertained by tangential points of dark-lines, and the tuing points of chaotic attractors can be also obtained by intersected points. The results show that this algorithm is helpful to analyse dynamic behaviour of systems and control chaotic motion.