C65.m

num=[0 0 0 1];
den=conv(conv([1 0],[1 2]),[1 4]);
rlocus(num,den),grid