A parasitic reconfigurable antenna array is a low-power approach for beamforming using passive tunable elements. Prior work on reconfigurable antennas in communication theory is based on ideal radiation pattern abstractions. It does not address the problem of physical realizability. Beamforming with parasitic elements is inherently difficult because mutual coupling creates non-linearity in the beamforming gain objective. We develop a multi-port circuit-theoretic model of the hybrid array with parasitic elements and antennas with active RF chain validated through electromagnetic simulations with a dipole array. We then derive the beamforming weight of the parasitic element using the theoretical beam pattern expression for the case of a single active antenna and multiple parasitic elements. We show that the parasitic beamforming is challenging because the weights are subject to coupled magnitude and phase constraints. We simplify the beamforming optimization problem using a shift-of-origin transformation to the typical unit-modulus beamforming weight. With this transformation, we derive a closed-form solution for the reconfigurable parasitic reactance. We generalize this solution to the multi-active multi-parasitic hybrid array operating in a multi-path channel. Our proposed hybrid architecture with parasitic elements outperforms conventional architectures in terms of energy efficiency.