The paper, based on the technique evaluation-interpolation, proposes a new numerical algorithm for the computation of the determinant of a bivariate polynomial matrix. The numerical algorithm uses the Newton bivariate polynomial interpolation with equidistant points. The algorithm is applied to the programming environment of MATLAB and is compared to the corresponding build-in function of MATLAB det() which uses symbolic operations.