On the numerical computation of Diophantine rotation numbers of analytic circle maps