En este trabajo analizamos variaciones y aplicaciones de un método numérico autoconsistente para estimar el estado fundamental de sistemas de espín 1/2 en dos dimensiones. El mencionado método se basa en el uso de la transformación de Jordan-Wigner sobre una red de espines bidimensional junto con una aproximación de campo medio multiparamétrica que permite resolver el sistema fermiónico resultante. En una primer etapa, investigamos redes cuadradas con condiciones de contorno abiertas y analizamos la incidencia de estas condiciones de borde en el estado fundamental del sistema. Luego, una vez establecida la viabilidad del método con condiciones abiertas, lo aplicamos en una segunda etapa, sobre una red bidimensional frustrada, obtenida al incorporar enlaces diagonales sobre la red cuadrada original. En esta última red observamos al variar el grado de frustración, una transición entre la fase de Neel y una fase con magnetización nula en el eje preferencial de cuantización.