Object-oriented programming is a solution for problems in the development, maintenance, and update of large software such as power system analysis software. However, many applications in the electrical industry critically depend on the computational efficiency of the implementation. In order to be flexible and reduce the computation time, this paper represents an efficient method for constructing a Jacobian matrix and for factorizing it, and designs the class hierarchy for power system. This method is applied to four different power systems for load flow calculation. The performance of the object-oriented program developed in C++ is assessed in computer simulation.