Можно по другому расписать
G11*u10-G12*u20-G13*u30=0
G33*u30-G32*u20-G31*u10=i(t)
входящий в узел ток берут с +, выходящий - с минусом.
я так тоже пробовал расписывать, у меня не сходятся токи с контурным методом,
может я не так расписываю проводимости?
G11=1/z14+1/z12
G12=1/z12
G13=1/z13
G33=1/z23+1/zн
G32=1/z23
G31=1/z13