Gregory–Newton interpolation formula