seqfreq

function [g] = sf(a,b)
c=[];
for i=1:length(b);
   clear g;
   g(1:b(i)) = [a(i)];
   c=[c g];
end
x=sum(b);
while (((nnz(diff(c))+nnz(diff(fliplr(c)))) < (2*(x-1))))
   r=randperm(x);
   g(1:x)=c(r(1:x));
   c=g;
end