numeric::linsolve([(15-va)/17 - va/8 - va/2],[va]) leftCurrent = (15-1.2903)/17; rightCurrent= 1.2903/8; Vth = 7*leftCurrent + 5*rightCurrent