def METHOD_NAME(k, ct_range=np.pi/2):
    """Expectation value of cos_theta for a given k"""
    ct = np.cos(np.linspace(0, ct_range, 500))
    pb = _prob_boltz(ct, k)
    return (ct*pb).sum()/pb.sum()