#include <iostream>
#include <ginac/ginac.h>

int main()
{
	using namespace std;
	using namespace GiNaC;

	// Replace numeric(1,2) with 0.5 to give a solution
	ex y1 = numeric(1) + numeric(1,2) * I;
	ex y2 = numeric(1);
	ex y3 = numeric(3, 5) + numeric(-4, 5)*I;

	cout << "Calculating..." << endl;
	ex res = Li( lst(1,1,2), lst(y1,y2,y3) );
	cout << res.evalf() << endl;

	return 0;
}