a=[21:40]; b=[25,32,28,38]; [cc,ii]=ismember(b,a); ii=sort(ii);
if sum(cc)~=length(b), error('sth is error !'); end
for j=1:length(ii)-1, pp=ii(j)+1:ii(j+1); a(pp)=a(pp)-j; end
if ii(end)~=length(a), j=j+1; pp=ii(j)+1:length(a); a(pp)=a(pp)-j; end
a(ii)=[]; a