delta = min([cspec(2:length(cspec))-cspec(1:length(cspec)-1); .01]);
delta = max(delta,.0001)
xs = 0:delta:max(cspec);

%{
%for j = 1:length(cspec)
%  y = cspec(j) ;
%  line([2,3],[y,y]) ;
%end
%
%axis([0 5 -1 cspec(length(cspec))+1])
%}

length(xs)

newys = [] ;
for j = 1:length(xs)
  newys(j) = sum(xs(j) > cspec) ;
end

plot(log(xs),newys,'-b') ;

%plot(xs,cspec+0.2*ones(size(cspec)),'.r') ;