%% prog017 exponencial segunda versión

clearvars;
close all;
clc;

%en este caso usaremos un ciclo for para tener la sumatoria
%% datos inciales

x = input('x = ');  %valor de x
N = input('N = ');  %numero de elementos de la sumatoria
yverd = cos(x);

%% calculo principal con ciclo for para la sumatoria
yaprox = 1;

for n = 1:N
    yaprox = yaprox + ((-1)^n)* x^(2*n)/factorial(2*n);
end;

%% error

errr = abs(yverd - yaprox);
           
%% resultados           
fprintf('x = %e yverd = %e yaprox = %e error %e\n',x, yverd, yaprox, errr);