def delta_bias(model, f_NL, **cosmo):
delta_bias = var.galaxy_bias(model, f_NL, **cosmo )[1]
semilogx(cc.k, delta_bias, next(linecycler), linewidth = 2, label = model + "-$f_{NL}$_%s" %f_NL)
legend(loc = 'best')
xlim((0.3 * 10**-3, 0.3* 10**-2))
xlabel('k')
ylabel('$\Delta b_g$(k, a = 1)')
def bias(model, f_NL, **cosmo):
bias = var.galaxy_bias(model, f_NL, **cosmo)[0]
loglog(cc.k, bias, next(linecycler), linewidth = 2, label = model+ "-$f_{NL}$_%s" %f_NL)
legend(loc = 'best')
xlim((0.3 * 10**-3, 0.3* 10**-2))
xlabel('k')
ylabel('$b_g$(k, a = 1)')
