toms886_test, a MATLAB code which calls toms886(), which implements an interpolation procedure based on "Padua points", defined in the square [-1,+1]^2, whose interpolating power is especially good. It is possible to map these points to the general rectangle, ellipse or triangle to do interpolation on these regions as well.
The computer code and data files made available on this web page are distributed under the MIT license
toms886, a MATLAB code which defines the Padua points for interpolation in a 2D region, including the rectangle, triangle, and ellipse, by Marco Caliari, Stefano de Marchi, Marco Vianello. This is a version of ACM TOMS algorithm 886.