Conventionally, polynomial filters are derived for evenly spaced points. Here, a derivation of polynomial filters for irregularly spaced points is provided and illustrated by example. The filter weights and variance reduction factors (VRFs) for both expanding memory polynomial (EMP) and fading-memory polynomial (FMP) filters are programmatically derived so that the expansion up to any degree can be generated. (Matlab was used for doing the symbolic weight derivations utilizing Symbolic Toolbox functions.) Order-switching and length-adaption are briefly considered. Outlier rejection and Cramer-Rao Lower Bound consistency are touched upon. In terms of performance, the VRF and its decay for the EMP filter is derived as a function of length (n) and the switch-over point is calculated where the VRFs of the EMP and FMP filters are equal. Empirical results verifying the derivation and implementation are reported.