The invention discloses a polar coordinate Newton method load flow calculation method suitable for research purposes, and provides a Newton method load flow calculation method which is easy to modify and maintain for researchers performing further research on the basis of polar coordinate Newton method load flow calculation. For a power system with n nodes, node numbers are irrelevant to node types, a Jacobian matrix for load flow calculation is stored according to (2n) * (2n), elements of rows corresponding to delta Q and columns corresponding to delta U of PV nodes are 0, elements of rows and columns related to balance nodes are 0, but diagonal elements are not reset, and original calculation values are reserved. Therefore, a linear equation solving function with excellent performance provided by a programming language can be used for solving a correction equation for load flow calculation, the programming difficulty is reduced, the calculation speed is improved, and the reliability of the algorithm is ensured. Implementation examples show that the correction equation can be solved by using a linear equation solving function provided by a programming language, the load flow calculation speed is increased, and the programming difficulty is reduced.