Se sabe que la ecuación de Schrödinger para los estados rotacionales de moléculas asimétricas es separable en coordenadas esferoconales e integrable en términos de funciones de Lamé. Sin embargo, la evaluación numérica de las últimas no se ha desarrollado eficientemente, limitando por lo tanto la aplicación práctica de tales soluciones. En este artículo, la evaluación matricial de los estados rotacionales se formula e implementa numéricamente para cualquier molécula asimétrica, usando la base familiar de armónicos esféricos. La matriz del Hamiltoniano - en un sistema de referencia fijo en la molécula y orientado a lo largo de los ejes principales - se construye en la base escogida y se muestra que se separa en bloques de ( l + 1) × ( l + 1) y l × l, para cada valor del número cuántico de momento angular l. La diagonalización de los bloques sucesivos conduce a valores precisos de las eigenenergías y los eigenvalores para todos los valores de los parámetros de asimetría. También se establece la conexión entre estos estados rotacionales y su representación en términos de funciones de Lamé, identificando al mismo tiempo una función generadora común para armónicos esféricos y armónicos esferoconales.The Schrödinger equation for the rotational states of asymmetric molecules is known to be separable in spheroconal coordinates and integrable in terms of Lamé functions. However, the numerical evaluation of the latter has not been developed efficiently, thereby limiting the practical application of such solutions. In this article, the matrix evaluation of the rotational states is formulated and implemented numerically for any asymmetric molecule, using the familiar bases of spherical harmonics. The matrix of the Hamiltonian - in a frame of reference fixed in the molecule and oriented along its principal axes - is constructed in the chosen basis and shown to separate into blocks of (l + 1) × (l + 1) AND l × l, for each value of the angular momentum quantum number l. The diagonalization of the successive blocks leads to accurate values of eigenenergies and eigenvectors for all values of the asymmetry parameters. The connection between these rotational states and their Lamé function representation is also established, identifying at the same time a common generating function for spherical harmonics and spheroconal harmonics.