Instead of matching the function at the block-corners, find the best interpolant minimizing the mean square difference between the approximation and the exact function. Still use finite elements, but bilinear approximations.