We implement in Matlab a Gauss-like cubature formula on bivariate
domains whose boundary is a piecewise smooth Jordan curve (curvilinear
polygons). The key tools are Green’s integral formula, together with
the recent software package Chebfun to approximate the boundary curve
close to machine precision by piecewise Chebyshev interpolation.
Several tests are presented, including some comparisons of this new
routine ChebfunGauss with the recent SplineGauss that approximates
the boundary by splines.