Polynomial Genetic Programming for Response Surface Modeling