Numerical integration based on optimal nodal splines