>>328
どううまくいかんのだろうか。
最後のIntegrate[..]の行を

Func[x_] := NIntegrate[R /. r -> x , {sita, 0, 2*Pi}]
Plot[Func[x], {x, 0.001, 10}, PlotDivision -> 2]

に変えろ。
元のグラフに特異点があるので(数値積分の)結果の正確さは
まったく保証できないが。