% f(x)=sin(x), x0=0, f^(n)(x0)=0,1,0,-1,0,...
x0 = 0 ;
x  = linspace(-pi, pi, 1000) ;
df = [0 1 0 -1 0 1 0] ;
y  = serie_Taylor(df, x, x0) ;
figure(1) ;
hold on ;
plot(x,sin(x)) ;
plot(x,y) ;
legend('sin(x)','Taylor (7 termini)') ;

% f(x)=exp(x), x0=0, f^(n)(x0)=1
x0 = 0 ;
x  = linspace(-pi, pi, 1000) ;
df = [1 1 1 1 1] ;
y  = serie_Taylor(df, x, x0) ;
figure(2) ;
hold on ;
plot(x,exp(x)) ;
plot(x,y) ;
legend('exp(x)','Taylor (5 termini)') ;