A simple and reliable algorithm for computing boundaries of power flow solutions due to system uncertainties