Circuit Theory/Thevenin-Norton/mupad3

numeric::linsolve([(15-va)/10 - (va-vb)/7 - (va-vb)/5 - va/3 = 0,
(va-vb)/7 - vb/2 + (va-vb)/5 = 0, 
iN = va/3 + (va-vb)/5],[va,vb,iN])