En este articulo se presenta un algoritmo para la solucion numerica de sistemas de ecuaciones no lineales. Para este proposito el sistema de ecuaciones se convierte en una funcion de valor real para luego ser minimizada, en el sentido global, en un dominio inicial dado (una caja en n) usando analisis de intervalos. El algoritmo disenado tiene la capacidad de determinar la existencia o no de soluciones al sistema de ecuaciones en una caja dada. Las soluciones del sistema de ecuaciones, si existen dentro de la caja dada, son expresadas mediante encerramientos por subcajas cuyo tamano es menor que la exactitud establecida. No hay restriccion acerca de la relacion entre el numero de ecuaciones y el numero de incognitas del sistema. Se realiza ademas un analisis de la convergencia del algoritmo y se muestran los resultados de su aplicacion para algunos problemas de prueba.