In the formulation of bus impedance/admittance matrix, the load (which can be converted to constant impedance Zload = V^2/(P -jQ) for example) are often not accounted for.

If we modeled the loads as constant impedance as above, How can they be included into the diagonal elements of the Zbus?