function floating(b,p,min,max) F=[0]; for m=1:(b^p-1) for e=min:max F=[F,m*b^(e-p),-m*b^(e-p)]; end end F = sort(F)' Y=zeros(length(F),1); plot(F,Y,'.');