Numerical integration based on spline approximation