Fully modified least squares estimation and inference for systems of cointegrating polynomial regressions