Uno de los problemas que aborda la Mecánica Celeste es la determinación de las figuras de equilibrio de los cuerpos celestes. Para investigar su solución mediante métodos directos, se precisa evaluar el potencial generado por su autogravitación, el generado por su fuerza centrí­fuga y el generado por la fuerza de atracción entre los cuerpos. Los métodos clásicos de Finlay y Kopal que afrontan estos problemas, para determinar el potencial autogravitatorio en las configuraciones de equilibrio, emplean desarrollos en serie de los potenciales interior y exterior del potencial autogravitatorio. Estos métodos incurren en el error de suponer la convergencia en capas donde resulta cuestionable dicha convergencia para estos desarrollos en serie. En este trabajo se han elaborado unos algoritmos que contemplan toda la casuí­stica y que permiten una manipulación eficiente del producto de polinomios de Legendre, del producto de funciones asociadas de Legendre y del producto de armónicos esféricos como combinacióon lineal de ellos mismos, respectivamente. Se han obtenido, para primer y segundo orden en las amplitudes, los desarrollos correctos para los potenciales interior y exterior del potencial autogravitatorio para configuraciones de equilibrio aisladas, y , en primer orden de amplitudes, los mismos potenciales para los sistemas binarios próximos. Se ha elaborado un método analí­tico, en primer orden respecto de las amplitudes, para la determinación del potencial de marea en sistemas binarios próximos en el cual se manifiesta la forma de la componente secundaria del sistema.