ap=1
as=30
Os=2.35
n=buttord(1,Os,ap,as,'s')
%format long
epsilon=sqrt(10^(ap/10)-1)
[Hcit,Hjme]=butter(n,epsilon^(-1/n),'s') %èitatel a jmenovatel pøen. funkce
[z,p,H0]=butter(n,epsilon^(-1/n),'s')    %nuly, póly a násobná konstanta