Polynomial operations: Numerical performance in matrix diophantine equation