def reduced_field(eta, a, mu0, E, kB, T):
    return eta * a**2 * mu0 * E / (kB * T)