  For any $\beta>0$, we provide a tridiagonal matrix model and compute the
joint eigenvalue density of a random rank one non-Hermitian perturbation of
Gaussian and Laguerre $\beta$-ensembles of random matrices.
