Numerical Solution of Integro-Differential Equations with Local Polynomial Regression