2D Barycentric Lagrange Interpolation

Interpolates a function on a rectangle.

This is a very fast, vectorized, 2D extension of barycentric lagrange interpolation for a product grid on a rectangle. See included file for demonstration of how to use the script.