[oneday2010]

syms x y n
x=-20:20;
y=x.^n/(sym('n!'))
n=10;
yy=double(subs(y));
plot(x,yy);