We present a method of solving functional equations of the type where f, F: P → P are unknown functions acting on an integral domain P and parameteres are given. We prove that under some assumptions on the parameters involved, all solutions to such kind of equations are polynomials. We use this method to solve some concrete equations of this type. For example, the equation (1) for f, F: ℝ → ℝ is solved without any regularity assumptions. It is worth noting that (1) stems from a well-known quadrature rule used in numerical analysis.