dobinskis formula