num = [0 0 10 10];
den = [1 4 0 0];
sys = tf(num, den);