code: def normalization(self):
          
          norm_r = self.normalization_range()        
          return quad(self, norm_r[0]*self._j_ref, norm_r[1]*self._j_ref)[0]
docs: The normalization 
          i.e., the intergral of the function over the normalization_range