function h = eigplot(poles, varargin)
% Plot eigenvalues

x = real(poles);
y = imag(poles);

h = plot(x, y, varargin{1:nargin-1});
