Abstract We introduce a new method to derive the orbital parameters of spectroscopic binary stars by a non-linear least squares analysis of residuals. Using measured radial velocity data of four double lined spectroscopic binary systems RZ Cas, CC Cas, V1130 Tau and Y Cygni, we find corresponding orbital and spectroscopic elements. Our numerical results are in good agreement with those obtained by others using more traditional methods.