Ojo !... La matriz del sistema de ecuaciones siempre tiene que ser simétrica respecto de la diagonal principal. Ya sea de mallas como de nodos.
La solución a tu planteo seria:
Para el nodo 1 :
0.66 V1 - 0.33 V2 - 0.33 V3 = 12 ............(12= suma de las corrientes entrantes)
Para el nodo 2 :
-0.33 V1 + 1.53 V2 - 1 V3 = -4 ......................salientes
Para el nodo 3:
-0.33 V1 - 1 V2 + 1.47 V3 = -5 .......................salientes
La conductancias propias siempre son (+).
La conductancia propia del nodo 2 = 1/3 + 1/5 + 1/1 = 1.53 mhos .... vos pusiste= -0.88
La conductancia propia del nodo 3 = 1/7 + 1/3 + 1/1 = 1.47 mhos .... vos
Pusiste 0.52
El resultado final estaría dando;
V1= 27 volts........................V2= 9 volts ..................V3= 8.7 volts,
Por cualquier duda volvés a consultarnos.