'check for bona fide density correction
currently only checks that density integrates to 1
`       non-negativity - NotImplementedYet'

def _verify(self):

    intdomain = self.limits
    self._corfactor = (1.0 / integrate.quad(self.evaluate, *intdomain)[0])
    return self._corfactor
