A method has been proposed for the combination of P-wave first-motion directions and S-wave polarization data for the numerical determination of earthquake focal mechanism. The method takes into account the influence of nearness of stations with inconsistent P-wave polarity observations, with respect to the assumed nodal planes. The mechanism solutions for six earthquakes selected from different geographic locations and depth ranges have been determined. Equal area projections of the nodal planes together with the P-wave first-motion and S-wave polarization data are presented for each earthquake. The quality of resolution of nodal plane determination on the basis of P-wave data, S-wave polarization, and the combination of P and S-wave data according to the present method, is discussed.