1.前言
今天在推导公式时，发现这破方程组人力根本计算不出结果，我就想可不可以使用编程解决，我发现MATLAB的solve()函数、和Python的sympy库都可以解决这个问题，在本篇文章中，使用Python sympy库来解决这个问题。
2. 问题
一个圆经过点（sqrt(3)(2r+t)/4，(2r+t)/4），（r，0），
另一个圆经过点（sqrt(3)(2r+t)/4，(2r+t)/4），（sqrt(3)(2r-t)/4，-(2*r-t)/4），
分别求两圆心坐标。
3.程序
3.1 导入库
# 导入sympy库
from sympy import *
3.2 定义字符
# 定义字符
a, b, r, t = symbols('a b r t')
3.3 确定方程
#将所有项移至左端，右端为零
eq1 = (sqrt(3)*(2*r+t)/4 - a)**2 + ((2*r+t)/4 - b)**2 - r**2
eq2 = (r - a)**2 + b**2 - r**2
eq3 = (sqrt(3)*(2*r-t)/4 - a)**2 + ((2*r-t)/4 + b)**2 - r**2
3.4 求解函数
# 求解函数，[eq1, eq2]为函数，[a, b]为未知数
c = solve([eq1, eq2], [a, b])
d = solve([eq1, eq3], [a, b])
3.5 输出函数
pretty为手写形式函数，更易读
# 输出结果，pretty为手写形式函数，更易读
print(c, '\n'*2, pretty(d))
