Integral and nonnegativity preserving approximations of functions

In this paper we consider the problem of approximating a function by continuous piecewise linear functions that preserve the integral and nonnegativity of the original function.