k = 1;
z = [];
p = [0,1,roots([1 1 0 1])'];
sys = zpk(z,p,k)