def Pk_bbks(self, k, z): """ Matter Power Spectra Pk in units if h^{-3}Mpc^{3} as a function of k in units of [h Mpc^{-1}] and z; Transfer function is taken to be BBKS Ref: Bardeen et. al., Astrophys. J., 304, 15 (1986) """ return self.A0bbks() * k**self.ns * tf.Tbbks( k, self.Omega_m_z(0), self.Ob0, self.h)**2. * self.D_plus_z(z)**2.
def integrand_bbks(self, k): return k**(self.ns + 2.) * tf.Tbbks(k, self.Omega_m_z(0), self.Ob0, self.h)**2. * (self.Wf(k))**2.