Based on Python and Abaqus platform, a stochastic finite element algorithm program was developed. The program can automatically discretize the non-stationary random field and calculate the failure probability when the geometric parameters of the slope and the mean, correlation distance and variation coefficient of the soil shear strength parameters are provided. The algorithm is used to calculate the benckmark slope examples, and the results verify the reliability of the proposed method. Moreover, the program can also effectively solve the complex problems of multiple software interaction and subroutine compilation. The program is then redeveloped on Abaqus platform, which is a widely applicable finite element software. It was beneficial to the popularization of the program in practical engineering applications.