//poles and zeroes s=%s sys=syslin('c',1/((2*s+1)*(1+5*s)*(s))) plzr(sys)