Instead of matching the function at the nodes, find the best interpolant minimizing the mean square difference between the approximation and the exact function. Still use finite elements, but linear approximations within elements.