Optimizing the milling process requires a priori knowledge of many process variables. However, the ability to include both milling stability and accuracy information is limited because current methods do not provide simultaneous milling stability and accuracy predictions. The method described within this paper, called Temporal Finite Element Analysis (TFEA), provides an approach for simultaneous prediction of milling stability and surface location error. This paper details the application of this approach to a multiple mode system in two orthogonal directions. The TFEA method forms an approximate analytical solution by dividing the time in the cut into a finite number of elements. The approximate solution is then matched with the exact solution for free vibration to obtain a discrete linear map. The formulated dynamic map is then used to determine stability, steady-state surface location error, and to reconstruct the time series for a stable cutting process. Solution convergence is evaluated by simply increasing the number of elements and through comparisons with numerical integration. Analytical predictions are compared to several different milling experiments. An interesting period two behavior, which was originally believed to be a flip bifurcation, was observed during experiment. However, evidence is presented to show this behavior can be attributed to runout in the cutter teeth.