BW=5;
[Af,Bf,Cf,Df]=ellip(8,0.1,90,BW*2*pi,'s');