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