Resolver esta ecuación en wolfram alpha,
… y su equivalente usando sympy
from sympy import symbols, var
v1,v2,v3,vo,r1,r2,r3,r4,rf=symbols('v1,v2,v3,vo,r1,r2,r3,r4,rf')
from sympy.solvers import solve
from sympy import init_printing
init_printing(use_unicode=False, wrap_line=False, no_global=True)
solve(v1/r1 + v2/r2 - ((v3*(r4/(r3+r4)))-vo)/rf, vo)
r1*r2*r4*v3 - r1*r3*rf*v2 - r1*r4*rf*v2 - r2*r3*rf*v1 - r2*r4*rf*v1
[-------------------------------------------------------------------]
r1*r2*(r3 + r4)