Sympyを用いた連立方程式解出力の形式について

Question

前提・実現したいこと

Pythonを使って高次元連立方程式を解き，結果をcsvやtxt形式で出力しようと思っています．

発生している問題・エラーメッセージ

下のプログラムを実行したところ，結果は以下のようになりました．

[{a: -1, b: 2}, {a: 1, b: 2}, {a: -2sqrt(5)/3, b: -5/3}, {a: 2sqrt(5)/3, b: -5/3}]

結果に平方根（sqrt)が含まれており扱いづらく感じています．
できればルート内も計算し，少数値として結果を取得したいのですが，上手なやり方がわからず困っております．

どなたかよい方法を教えていただけませんか．

該当のソースコード

import sympy as sp

a = sp.Symbol('a')
b = sp.Symbol('b')

#数式は適当なものを当て込んでいます．
expr1 = a**2 + b**2 -5
expr2 = 3*a**2 + b-5

result = sp.solve([expr1, expr2])

print(result)

Answer 1

数値解がほしいということでしょうか？
その場合、sp.N() で解析解を数値解に変換できます。

import sympy as sp

a = sp.Symbol("a")
b = sp.Symbol("b")

# 数式は適当なものを当て込んでいます．
expr1 = a ** 2 + b ** 2 - 5
expr2 = 3 * a ** 2 + b - 5

result = sp.solve([expr1, expr2])

result = [{k: sp.N(v) for k, v in sol.items()} for sol in result]
print(result)
# [{a: -1.00000000000000, b: 2.00000000000000},
#  {a: 1.00000000000000, b: 2.00000000000000},
#  {a: -1.49071198499986, b: -1.66666666666667},
#  {a: 1.49071198499986, b: -1.66666666666667}]

Comments

[monkey.0435]

早速のご回答ありがとうございます！
望みの形式で出力することができました！