function y = inhour(lambda, Gamma, beta, rho)
% Find the roots of the inhour equation
% function y = inhour(lambda, Gamma, beta, rho)
% (c) 2007 Alain Hebert, Ecole Polytechnique de Montreal
  [nn, dd]=residue(lambda.*beta,-lambda,[]) ;
  y=sort(roots(conv([-Gamma, rho-sum(beta)],dd)+[0 0 nn])) ;