Gauss forward interpolation formula