var('vi1 ik il r1 r2 r3')
eq1=Eq(0+vi1-r1*(1e-3-il)-r2*(1e-3-ik),0)
eq2=Eq(0-r2*(ik-1e-3)-3-r3*ik,0)
eq3=Eq(vi1-2*(-r2*(ik-1e-3))+3-r1*(il-1e-3)
ans=solve([eq1,eq2,eq3,vi1)],[vi1,ik,il])
ans[vi1].subs({r1:1e3,r2:2e3,r3:1e3})
5.00000000000000
ans[ik].subs({r1:1e3,r2:2e3,r3:1e3})
-0.000333333333333333
ans[il].subs({r1:1e3,r2:2e3,r3:1e3})
-0.00133333333333333
(-r2*(ac[ik]-1e-3)).subs({r1:1e3,r2:2e3,r3:1e3})
2.66666666666667