This work presents dynamic stability analysis of general high speed milling process by higher order semi-analytical method. The first order full-discretization method (1 st FDM) and second order full-discretization method (2 nd FDM) are presented. These methods based on the direct integration scheme. The governing mathematical model applied is the delay-differential equation with single time periodic delay taking regenerative effect into account. The stability lobes diagrams are presented for single degree of freedom mechanical model and two degree of freedom mechanical model. Up-milling and down-milling stability chart are presented for various radial immersion ratio in order to compare accuracy of up-and down-milling. It is demonstrated that for full-immersion down-milling and up-milling stability lobes diagram are the same approximatively but in other cases downmilling is more accuracy than up-milling both for single degree of freedom and two degree of freedom.the computational time of calculation of eigenvalue is also variable for different computational parameter. The rate of convergence for half immersion and low immersion is presented for variable parameters.