In this paper, a probabilistic method is proposed to analyze the steady-state operating conditions of an active electrical distribution system with Wind (WD) and Photovoltaic (PV) generation plants. This method takes into account the uncertainties of power load demands and power production from renewable generation systems and combines Monte Carlo simulation techniques and multi-linearized power flow equations. The power flow equations include models of wind turbine and PV generation units and multi-linearization is accomplished by applying a criterion based on the total active power of the system. The method properly extends a probabilistic method proposed in the relevant literature for traditional passive electrical distribution systems to the field of an active electrical distribution system with WD and PV generation units. Numerical applications are presented and discussed with reference to a 17-bus test distribution system characterized by WD and PV systems connected at different busbars. The results obtained with the proposed algorithm are compared with the results obtained using a Monte Carlo simulation algorithm that included non-linear power flow equations. (C) 2014 Elsevier Ltd. All rights reserved.