In this paper, we consider the problem of eigenvalue computation in a class of dynamical system equations arising from modelling turbulent dispersion. An alternative numerical approach is suggested for computing eigenvalues. The dynamical system equations are considered as differential-algebraic (DA) equations. The eigenvalues of the DA equations are then calculated as roots of the boundary condition at the end point. Numerical examples are included to illustrate the efficiency of the suggested approach.