Hyper2d: A numerical code for hyperinterpolation on rectangles

Hyperinterpolation at Morrow–Patterson–Xu cubature points for the product Chebyshev measure provides a simple and powerful polynomial approximation method on rectangles. Here, we present an accurate and efficient Matlab/Octave implementation of the hyperinterpolation formula, accompanied by several numerical tests.