Yeah I saw that right away, thats just a typo it's node1+3
Now after fixing the big mistake, which was just G9 should equal = G4+G5, now it works. So without my extra composite G's which aren't labelled the best, since I was just trying to save myself page space as I went
for node 2, eqn 1 would be
-(2*G3 + G4)*Va = -(G2 + 2*G3 + G4)*V1 + (G2 + 3*G3 + G4)*V2
and the supernode (1+3) would give eqn 2
iv - iu + (G4 + G5)*Va = (G1 + G2 + G4 + G5)*V1 - (G2 + G4)*V2
now finally it works, but I tried it 2-3 times and probably just had a silly mistake each time.
I need to be more careful, I was worried the method was wrong, and just working by fluke sometimes, so I thought maybe my assigned currents had to match the way I used them at node 2 or IDK...and in this case I did make them match