Poisson–Charlier function