Legendre polynomials using GiNaC

legendre.cc

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

using namespace std;
using namespace GiNaC;

ex LegendrePoly(const symbol &x, int n)
{
    ex PKer = pow(x * x - 1, n);
    return normal(1 / (pow(2, n) * factorial(n)) * diff(PKer, x, n));
}

int main()
{
    const int n = 10;

    symbol x("x");

    for (int i = 0; i < n; ++i)
        cout << "P_" << i << "(x) == " << LegendrePoly(x, i) << endl;

    return 0;
}