  We propose to use eigenvalue densities of unitary random matrix ensembles as
mass distributions in gravitational lensing. The corresponding lens equations
reduce to algebraic equations in the complex plane which can be treated
analytically. We prove that these models can be applied to describe lensing by
systems of edge-on galaxies. We illustrate our analysis with the Gaussian and
the quartic unitary matrix ensembles.
