In this paper we consider computing determinants of polynomial matrices symbolically. Determinant computation of matrices with polynomial entries in a small number of variables is of particular interest since it commonly appears in solving engineering design problems. A parallel algorithm based on multivariate Newton polynomial interpolation with "cut-surface" (total degree bound) is presented and its efficiency is demonstrated by showing computational results for some practical examples from control system design.