Hello,

I am trying to replicate the figure. 7 in the article arXiv:1103.2683 which is in the attachment. However there are some problems in my mathematica code. First there are divergences near \rho==2 and the plot mismatches the one in the aricle arXiv:1103.2683 near \rho==0. Second if we look the plot carefully we find that there are some strange points which seems that the alogrithm isn't stable.

ListPlot[TabrAreaR, AxesOrigin -> {0, 0},  PlotRange -> {{0, 2.5}, {-2, 2}}]

The article and my mathematica notebook are in the attachment. I don't know how to deal with these problems. I suspect that we should resolve the two warnings in NIntegrate and NDSolve. Thank you very much!

Best regards, Xu