Matlab - funcion en Matlab u Octave que devuelve los coeficientes- del polinomio

function legendreP = legendre(n)

syms x

    if n == 0

        legendreP = 1;

    elseif n == 1

        legendreP = 1-x;

    else

        legendreP = expand(((2*n-1-x)*legendre(n-1)- (n-1)*legendre(n-2))/n);

    end

>> legendreP = legendre(3)

legendreP =

- x^3/6 + (3*x^2)/2 - 3*x + 1