Actualizo para separar el procedimiento buscar SolucionNR en varios para poder controlar NR en cada iteracion, se agrega ademas funciones en matrices esparsas para restar un coeficiente, (resta_my) esto para poder modificar la matriz MY del problema entre iteraciones
