%
clear
close
clc
%%
% f(x)=exp(x) x0=0
aa=-2;bb=2;
xx=linspace(aa,bb);
yy=zeros(1,length(xx));
a=(1)./factorial(0:5);
for i=1:length(a)
    figure(1)
    clf
    hold on
    yy=yy+a(i)*(xx-0).^(i-1);
    %
    hold on
    title(['f(x)=e^x   x_0=0   Polinomio de Taylor grado ',num2str(i-1)],'fontsize',16)
    axis([aa bb -1 3])
    plot(0,1,'k*')
    plot([aa bb],[0 0],'k-')
    box on
    plot(xx,exp(xx),'linewidth',2)
    pause
    plot(xx,yy,'r-','linewidth',2)
    pause
end
clear
close all
clc
%%
% f(x)=sin(x) x0=pi
aa=0;bb=7;
xx=linspace(aa,bb);
yy=zeros(1,length(xx));
a=[0 -1 0 1 0 -1];
a=a./factorial(0:5);
for i=1:length(a)
    yy=yy+a(i)*(xx-pi).^(i-1);
    %
    figure(2)
    clf
    hold on
    axis([aa bb -1 1])
    title(['f(x)=sen(x)   x_0=\pi   Polinomio de Taylor grado ',num2str(i-1)],'fontsize',16)
    plot(pi,0,'k*')
    plot([aa bb],[0 0],'k-')
    box on
    plot(xx,sin(xx),'linewidth',2)
    pause
    plot(xx,yy,'r-','linewidth',2)
    %
    pause
end