Abstract: Based on the rule of FEM software command of ANSYS, and the function library of MATLAB and its real-time visualization, we accomplished an irregular boundary FEM model with according to the tectonic boundary of North China. Furthermore, based on the uniqueness of FEM node coordinate and its number indication, we made use of MATLAB to produce APDL code to put the irregular load on the boundary nodes and to put the irregular dislocation vectors on inner source. We also provided an example on calculation of Coulomb stress. In the post-processing of result plotting, with the help of function of ellipsoid in MATLAB, we transferred the node coordinate in Cartesian into the geographical coordinates. Thereafter we could display the numerical simulation result in numerous general geographical software regardless of subtle image or deep comparison.