Todas las fundaciones de los edificios en zonas de peligro sísmico deben analizarse de acuerdo con la carga vertical resultante excéntrica. El propósito de este trabajo es desarrollar una aplicación computacional para calcular la distribución biaxial de presiones o tensiones en la base de zapatas rectangulares rígidas por el método de Czerniak (1963). Según el lugar donde se aplica la carga vertical resultante en zapatas, esta puede localizarse por una parte dentro del núcleo central y, por otra, fuera del mismo con cuatro posibles ubicaciones referidas a los cuadrantes respecto del centroide, formado así cinco zonas de análisis. Para establecer las presiones sobre la superficie de contacto debe definirse la ubicación del eje neutro del sistema. Dado que la posición del eje neutral no se conoce inicialmente se propone una aplicación computacional que realice un proceso de aproximaciones sucesivas hasta que la convergencia de resultados establezca el porcentaje de área en contacto con la superficie de suelo. La ejecución del procedimiento de cálculo se muestra a través de ejemplos.