I am looking to find the proper way to include PV generation at load buses of IEEE 14 bus system to analyse its voltage stability. My main focus is on line losses, loadability margin, reactive power margin, bus participation factor and I would like to see how PV penetration impacts these factors. I have developed MATLAB m code for the base system and want to move to PV penetration.