The centralized model of electricity generation cannot meet the growing demand of users at all times due to the low availability of fossil fuels, their high cost, and the environmental pollution they cause. Most distributed generation from renewable sources is integrated into radial distribution systems due to its availability, low operating cost, and absence of pollutants. Integrating photovoltaic distributed generation in a random location and size can increase losses, costs, and large fluctuations in the voltage profile. So, studying the optimal inclusion of photovoltaic distributed generation in distribution systems is essential to provide consumers with reliable and safe electricity. This paper presents the Beta Hill Climbing algorithm to solve the problem of locating and sizing photovoltaic distributed generation in radial distribution systems. Annual active power loss reduction, voltage deviation, and voltage stability index were used as objective functions. To validate the algorithm, the standard IEEE-33 circuit was used in the MATLAB environment and its results were compared with those obtained by genetic algorithms and particle swarm optimization.