Try
For[r = 1, r <= 32, r++, For[j = 1, j <= Length[XandY], j = j + 0.1, If[Sqrt[(y[j] - 4.8)^2 + (x[j] - 2.6)^2] <= r, plotCount[[r]] = plotCount[[r]] + 1]]]
Thanks, I just had a brainwave though! Basically just used <= r/10
Cheers for the response!
