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