Recent, problems on the voltage-instability have been paid attention in power system and methods to find the limit of voltage-stability, concerned with these problems, were developed. However, these methods are short of precision on the limit of voltage-instability. Here, using the second-order load flow, constraint equation(d Pi/d Vi=0) and its patial differentiations are precisely formulated. Also, since the taylor series expansion of power flow equations terminates at the second-order terms, partial differentiations of constraint equation, that is Hessian, are constant. Then, Hessian matrix are calculated once during iteration process.