code: def normal_log_q(self,z):
                       
          means, scale = self.get_means_and_scales()
          return ss.norm.logpdf(z,loc=means,scale=scale)
docs: The unnormalized log posterior components for mean-field normal family (the quantity we want to approximate)
          RAO-BLACKWELLIZED!