Circuit Theory/Thevenin-Norton/mupad

numeric::linsolve([
(15-va)/10 - (va-vb)/7 - (va-vc)/12,
(va-vb)/7 - vb/2 - (vb-vc)/3,
(va-vc)/12 + (vb-vc)/5 - vc/3],[va,vb,vc])

va = 5.392644577; vb = 1.167273696; vc = 1.107311006;
(va-vc)/12