function lm = CS4640_local_max_1D(h)
%

N = length(h);
lm = zeros(1,N);

for p = 2:N-1
    if h(p)>=max([h(p-1:p+1)])
        lm(p) = h(p);
    end
end
