Как оценить символическую переменную в SymbolicC ++?

Я пробовал пару библиотек компьютерной алгебры для C ++, чтобы использовать с курсом векторного исчисления, который я беру. У меня проблемы с нелинейными уравнениями в GiNaC и в SymbolicC ++ это действительно сработало.

Вот простой пример, но проблема в том, что я не могу понять, как вычислить число и, возможно, привести его к двойному или плавающему:

#include <iostream>
#include "symbolicc++.h"
using namespace std;
int main(void)
{
Symbolic x("x"), y("y");

Equation e1 = (x^2) + (y^2) == 13;
Equation e2 = (x^2) - y == 7;

Equations eqs = {e1, e2};
list<Symbolic> symbs = {x, y};
list<Equations> sols = solve(eqs, symbs);

Symbolic x_sol, y_sol;
int i = 1;
for( auto iter1 = sols.begin(); iter1 != sols.end(); iter1++)
{
x_sol = x.subst((*(*iter1).begin()));
y_sol = y.subst((*(--(*iter1).end())));
cout << "p" << i << " = {" << x_sol << ", " << y_sol << "};" << endl;
i++;
}
return 0;
}

С этим выводом я могу скопировать и вставить его в ginsh и он оценивает просто отлично, но остается в развернутом виде в SymbolicC ++.

Точный вывод, который я получаю, выглядит следующим образом:

p1 = {1/2*(-2*(25)^(1/2)+26)^(1/2), -1/2*(25)^(1/2)-1/2};
p2 = {1/2*(2*(25)^(1/2)+26)^(1/2), 1/2*(25)^(1/2)-1/2};
p3 = {-1/2*(-2*(25)^(1/2)+26)^(1/2), -1/2*(25)^(1/2)-1/2};
p4 = {-1/2*(2*(25)^(1/2)+26)^(1/2), 1/2*(25)^(1/2)-1/2};

Как я могу оценить подобные выражения и привести их к doubles?