Expansion over a system of shifts of pyramidal functions

A method for approximation of functions of two variables by a linear combination of nonnegative piecewise linear functions with a compact support is presented. The crucial idea of this method consists in an “integral” calculation method for the coefficients. The accuracy of the approximation in the spaces of continuous and integrable functions is proved to have the same order as the best approximation by piecewise linear functions.