Distributed generation including wind turbine (WT) and photovoltaic panel increased very fast in recent years around the world, challenging the conventional way of probabilistic load flow (PLF) calculation. Reliable and efficient PLF method is required to take into account such changing. This paper studies the maximum entropy probabilistic density function reconstruction method based on cumulant arithmetic of linearized load flow formulation, and then develops a maximum entropy based PLF (MEPLF) calculation algorithm for power system integrated with wind power generation (WPG). Comparing to traditional Gram–Charlier expansion based PLF (GC-PLF) calculation method, the proposed ME-PLF calculation algorithm can obtain more reliable and accurate probabilistic density functions (PDFs) of bus voltages and branch flows in various WT parameter scenarios. It can solve the limitation of GC-PLF calculation method that mistakenly gaining negative values in tail regions of PDFs. Linear dependence between active and reactive power injections of WPG can also be effectively considered by the modified cumulant calculation framework. Accuracy and efficiency of the proposed approach are validated with some test systems. Uncertainties yielded by the wind speed variations, WT locations, power factor fluctuations are considered.