def compute_zT(rho, seebeck, kappa, temperature):
    """Returns the thermoelectric figure of merit.
    Note: All inputs should have SI units.
    """
    return seebeck ** 2 * temperature / (rho * kappa)