Gauss interpolation formula