Python解算多项式

from sympy import *

#定义符号变量
a,b,c,a1,a2,a3,b1,b2,R,x = symbols("a,b,c,a1,a2,a3,b1,b2,R,x")

#公式
a = ((-(4*R**2*(cos(b1))**2*(sin(b1)*(pow((x**2)*sin(a3)**2-R**2,2))+cos(b1)*R)**2)/x**2*sin(a2)**2*sin(a3)**2+R**2+(4*R*cos(b1)*(sin(b1)*pow((x**2*sin(a3)**2-R**2),2))+cos(b1)*R)/sin(a2)+sin(a3))/x**2)

#解算
result = simplify(a)
pprint(result, use_unicode=True) #打印标准格式

#写入到文件
with open('test.txt','a') as file0:
    print(result,file=file0)