This MATLAB function returns the zeros z, poles p, and gain(s) k of the zero-pole- gain model sys.
www.mathworks.com