Numerical ranges of KMS matrices