function s = range(T,W, p, pm)
  v=[];
  for i=1:length(pm)
   v = [v ((pm(i)-1)/(p*W))*(3e8/2)];
  end
  s=v;