def p(self,mu): return self.Nv*fermi_dirac_one_half((self.Ev-mu)/kT)
def p(self, mu): return self.Nv * fermi_dirac_one_half((self.Ev - mu) / kT)
def n(self,mu): return self.Nc*fermi_dirac_one_half((mu-self.Ec)/kT)
def n(self, mu): return self.Nc * fermi_dirac_one_half((mu - self.Ec) / kT)