%% prog016 exponencial primera versión

clearvars;
close all;
clc;

%% datos inciales

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

%% primera aproximación

yaprox = 1 + x^1/factorial(1) ...
           + x^2/factorial(2) ...
           + x^3/factorial(3) ...
           + x^4/factorial(4) ...
           + x^5/factorial(5) ...
           + x^6/factorial(6) ...
           + x^7/factorial(7) ...
           + x^8/factorial(8) ...
           + x^9/factorial(9) ...
           + x^10/factorial(10);

%% error

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