function h=ideal_lp(M,wp);

k=M-1
s1 = k/2;
s2 = [0:1:k];
s3 = s2 - s1 + eps;
h = sin(wp*s3) ./ (pi*s3);

