四元二次方程组的求解

# 4元2次方程组的计算。
# 应用在惠斯通电桥测电阻后求每个独立电阻的阻值。
import sympy


a = 1.75  # R1//(R2+R3+R4)
b = 4.00  # R2//(R1+R3+R4)
c = 1.75  # R3//(R1+R2+R4)
d = 3.00  # R4//(R1+R2+R3)
x1, x2, x3, x4, x_sum = sympy.symbols("x1, x2,x3,x4,x_sum",positive=True)
eq0 = x1 + x2 + x3 + x4 - x_sum
eq1 = x1 * (x_sum - x1) - a * x_sum
eq2 = x2 * (x_sum - x2) - b * x_sum
eq3 = x3 * (x_sum - x3) - c * x_sum
eq4 = x4 * (x_sum - x4) - d * x_sum


s = sympy.solve([eq0, eq1, eq2, eq3, eq4], [x1, x2, x3, x4, x_sum], dict=True)
print(s)