We use the method introduced by Karami & Mohebi (2007), and Karami & Teimoorinia (2007) which enable us to derive the orbital parameters of the spectroscopic binary stars by the nonlinear least squares of observed vs. curve fitting (o–c). Using the measured experimental data for radial velocities of the four double-lined spectroscopic binary systems PV Pup,HD141929, EE Cet and V921 Her,we find both the orbital and the combined spectroscopic elements of these systems. Our numerical results are in good agreement with those obtained using the method of Lehmann-Filh´es.