Abstract: The split coefficient matrix method is one of finite difference methods, which performs discretization on the convective term according to eigenvalues of the equations. That method can describe the real information propagation direction of flow field. The system equations were integrated based on the split coefficient matrix method in the paper. The generated block triple diagonal algebra matrix was solved with the method of forward elimination and backward substitution. The developed code was assessed by several separated effect tests. These tests were carried out through convergent-divergent nozzles and two-phase flow compression test. Good agreements between the computation and analysis results prove the validity and applicability of the method. The split coefficient matrix method was successfully applied to solve the two-phase flow equations. The results provide a new point of view for the thermal-hydraulic code development.