In this paper, a probabilistic load flow analysis is proposed in order to deal with probabilistic problems related to the power system. Due to increasing trend of penetration of renewable energy sources in power system brought two factors: One is uncertainty, and another one is dependence. Uncertainty and dependence factor increase risk associated with power system operation and planning. In this proposed model these two factors is considered. Gaussian Copula theory is proposed to establish the probability distribution of correlated input random variables. Three sampling methods are used with Monte Carlo simulation as simple random sampling, Box-Muller sampling, and Latin hypercube sampling in order to evaluate the accuracy of the proposed method. The main advantages of this model are as: It can establish any type of correlation between input random variable with the help of Copula theory, it is free from the restrictions of Pearson coefficient of correlation, it is unconstrained by the marginal distribution of input random variables, and uncertainty is established with photovoltaic generation this is the main source of uncertainty. Additional, in order to evaluate the accuracy and efficiency of the proposed model a real load and photovoltaic generation data is adopted. For accuracy evaluation purpose two comparative test system is adopted as modified IEEE 14 and IEEE 118-bus test system.