a = zeros(1,50);

for k=1:24
    a(2*k-1) =(-1)^(k-1)/factorial(k-1);
end

for k=1:50
        fprintf('%10.6f ', a(k));
end

x =  -1:0.1:1;

plot(x, exp(-x.^2), 'b-');
hold on
plot(x, Evalua_Polinomio(x, a), 'r+');



