%% prog016 exponencial primera versión

clearvars;
close all;
clc;

%% datos inciales

x = input('x = ');
yverd = cos(x);

%% primera aproximación

yaprox = 1 - x^2/factorial(2) ...
           + x^4/factorial(4) ...
           - x^6/factorial(6) ...
           + x^8/factorial(8) ...
           - x^10/factorial(10) ...
           + x^12/factorial(12) ...
           - x^14/factorial(14) ...
           + x^16/factorial(16) ...
           - x^18/factorial(18) ...
           + x^20/factorial(20);

%% error

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