xinit('fstyp12.1.ps')
    hd=[0*ones(1,15) ones(1,10) 0*ones(1,39)];//desired samples
    hst1=fsfirlin(hd,1);//filter with no sample in the transition
    hd(15)=.5;hd(26)=.5;//samples in the transition bands
    hst2=fsfirlin(hd,1);//corresponding filter
    pas=1/prod(size(hst1))*.5;
    fg=0:pas:.5;//normalized frequencies grid
    n=prod(size(hst1))
    plot(fg(1:n),hst1);
    plot2d(fg(1:n)',hst2',[-1],"000");

