simpson integration in expectation() and normalise