function [sig]=sigma(R)
sig=[];
for s=0:0.5:6
    sig=[sig sqrt(1/2/R/10^(s/10))];
end