Approximation by gridding
Using "arildno"'s method
%