Abstract. Machining accuracy of a milled surface is influenced by process dynamics. Surface location error (SLE) in milling determines final dimensional accuracy of the finished surface. Therefore, it is critical to predict, control, and minimize SLE. In traditional methods, the effects of uncertain factors are usually ignored during prediction of SLE, and this would tend to generate estimation errors. In order to solve this problem, this paper presents methods for probabilistic analysis of SLE in milling. A dynamic model for milling process is built to determine relationship between SLE and cutting parameters using full-discretization method (FDM). Monte-Carlo simulation (MCS) method and artificial neural network (ANN) based MCS method are proposed for predicting reliability of the milling process. Finally, a numerical example is used to evaluate the accuracy and efficiency of the proposed method.