At the beginning you have 3 variables unknown and 2 equations so you're able to solve it for S
and D
. There are contains 3 pairs of solutions of the form (f(V), g(V)) where f and g are expressions that contains V. Then you can iterate these pairs and use a substitution:
import sympy as sp
S, D, Vdu = sp.symbols('S D Vdu')
eq1 = sp.Eq(S, D)
eq2 = sp.Eq(Vdu, sp.pi**2*S*D**2/4)
solution = sp.solve((eq1, eq2), (S, D))
substituted_solution = [(s.subs({Vdu : 250e-6}), d.subs({Vdu : 250e-6})) for s, d in solution]
Thank you very much I was solved with two substitutions using .subs().subs(), but your method is more useful for me. I'm very grateful, thanks again :)