def P_Phi_p(fv, **cosmo):
P_Phi_p = power_primordial(fv, **cosmo)
loglog(cc.k, P_Phi_p, next(linecycler), linewidth = 2, label = '$k^{-2}$-factor: %s' %fv)
legend(loc = 'best', prop = {'size':9})
axvline(x = 0.17 * 10**-2, color='r', ls = '--') # vertical line at k_eq
axvline(x = 0.14 * 10**-1, color='b', ls = '-.') # vertical line at k_eq
xlim((cc.k_min, cc.k_max))
xlabel('k')
ylabel('$P^p_{Phi}$(k)')
def P_Phi(model, fv, **cosmo):
P_Phi_p = power_primordial(fv, **cosmo); g_Delta_m = var.grwoth_factor(model,**cosmo)[0]
g_Phi = var.grwoth_factor(model,**cosmo)[1]
P_Delta = (g_Delta_m * g_Phi)**2 * P_Phi_p
loglog(cc.k, P_Delta, next(linecycler), linewidth = 2, label = '$k^{-2}$-factor: %s' %fv)
legend(loc = 'best', prop = {'size':9})
axvline(x = 0.17 * 10**-2, color='r', ls = '--') # vertical line at k_eq
axvline(x = 0.14 * 10**-1, color='b', ls = '-.') # vertical line at k_eq
xlim((cc.k_min, cc.k_max))
xlabel('k')
ylabel('$P_m$(k, a = 1)')
