On the higher rank numerical range of the shift operator

For any n-by-n complex matrix T and any $1\leqslant k\leqslant n$, let $\Lambda_{k}(T)$ the set of all $\lambda\in \C$ such that $PTP=\lambda P$ for some rank-k orthogonal projection $P$ be its higher rank-k numerical range. It is shown that if $\bbS$ is the n-dimensional shift on ${\C}^{n}$ then its rank-k numerical range is the circular disc centred in zero and with radius $\cos\dfrac{k\pi}{n+1}$ if $1