Exact Integration for singular Zienkiewicz and Guzman-Neilan Finite Elements with Implementation

We develop a recursive integration formula for a class of rational polynomials in 2D. Based on this, we present implementations of finite elements that have rational basis functions. Specifically, we provide simple Matlab implementations of the singular Zienkiewicz and the lowest-order Guzman-Neilan finite element in 2D.