An easy-to-use variant of the finite difference method (FDM), the polynomial finite difference method (PFDM), for the numerical solution of ordinary differential equations (ODEs) and differential-algebraic equations (DAEs), is presented. Compared to the traditional implementation of the FDM, the PFDM approach has two major advantages: straightforward implementation, and easily adjustable accura...