clear
clc
close all

%%% Creamos un vector x

x = -2*pi:0.1:2*pi;

y = cos(x);

fig = figure();
set(fig, 'color', 'white')
plot(x,y, 'LineWidth', 2)
grid on
xlabel('x')
ylabel('y')

%%% Utilizamos el vector x para hacer una estimación del coseno
%%% utilizando la expancion en serie de Taylor

N = 1;

yest = 0;

for k=0:N
    yest = yest + (-1)^k*(x.^(2*k))./factorial(2*k);
end;

hold on

plot(x, yest, 'r-', 'LineWidth', 2);

legend('Funcion original', 'Funcion Aproximada Taylor')