A new method to calculate multiple power flow solutions