A general equation by which the Hessian matrix of an error function can be determined directly, has been derived. It was verified to be useful in optimization processes that include the Hessian matrix. A few design examples had shown that this method had accelerated the processes of finding the minimums. The advantage of this technique is the possibility of optimizing functions that composed of both the phases and magnitudes.