fp=2000;fs=10000;
fm=44100;fny=fm/2;
Ap=0.5;As=50;
ap=1-10^(-(Ap-0.4)/20);
as=10^(-(As+10)/20);

[N,Fo,Mo,W]=remezord([fp fs],[1 0],[ap as],44100);
B=remez(N,Fo,Mo,W);

[H,F]=freqz(B,1,500,fm);
plot(F,abs(H));hold on;zoom;grid;
%plot([fp,fs],-[Ap As],'*r');hold off;
fprintf(1,'N = %d\n',N);