Integrated modified least squares estimation and (fixed-b) inference for systems of cointegrating multivariate polynomial regressions