Abstract: In order to solve the problems existing in subsidence prediction for irregular working surface, through the research on the triangulation method in irregular working face, the upper and lower bounds of each integral region for double integrals are derived. Considering the different shapes of the triangle and the conversion between different coordinate systems, a triangulation algorithm for irregular working face is designed. A subsidence prediction program is developed based on MATLAB software and the triangulation algorithm. The program is not only applicable to the surface subsidence prediction of gently inclined working surface under different shapes, but also used to process the results directly without needing to use the third party auxiliary software, which can realize the real time visualization of the subsidence prediction data.