Matlab - Euler

Saludos

estoy tratando de demostrar la cte de euler exp^jθ =-1
usando taylor sabemos que esa exponencial es igual a
Ʃ=(x^n)/n! de cero hasta inf.
hay es donde aparece mi duda, como puedo hacer para que la funcion tome el valor del factorial?
ya he creado la funcion [factorial (n)]
syms x n
f=((x^n)/factorial(n)
symsum (f,n,0,inf)

despues de correr el programa marca un error que no se reconoce la funcion factorial (n)
alguien puede orientarme en el caso....
gracias

aux=1;
for n=2:x
aux=aux*n;
end

Esta es la forma rápida, aux será el factorial de x(mayor de 2, aunque cambiando l dos por un 1 tambien te vale, no olvides inicializar aux), la forma elegante es la recursividad.