The new algorithm of solving harmonic balance equations which used in electronic CAD systems is presented. The new algorithm is based on implementation to harmonic balance equations the ideas of model order reduction methods. This algorithm allows significantly reduce the size of memory for storing of model equations and reduce of computational costs.