op=find_maximums(s,dz,no)

function op=find_maximums(s,dz,no)
% op=find_optimums(s,dz,no)
% s -signal
% dz - dead zone, +/-
% no - number of optimums to find
% op - found optimums


L=length(s);
lz=true(L,1); % leave zone, all free at ther begining
op=zeros(no,1);
for noc=1:no % optimums counter
    f=find(lz);
    [mx mxi]=max(s(f));
    op1=f(mxi);
    op(noc)=op1;
    i1=op1-dz;
    if i1<1
        i1=1;
    end
    i2=op1+dz;
    if i2>L
        i2=L;
    end
    lz(i1:i2)=false; % decrease leave zone near optimum
end