code: def k_bn(self, n, Re):
          
          bn = self.b_n(n)
          k = bn*Re**(-1./n)
          return k, bn
docs: returns normalisation of the sersic profile such that Re is the half light radius given n_sersic slope