A numerical method is outlined for computing roots of determinants of lambda matrices. Convergence of the method is quadratic as long as the derivative of the determinant does not vanish at the root. When the derivative is zero, the method may still converge in special cases. Three examples of mechanics problems giving rise to lambda matrices are included.