%% prog20230421 esponencial version 00

clearvars;
clear all;
close all;
clc;

%% lectura de datos
x = input('numero al que se calculara la funcion exponencial x = ');
N = input('numero de componentes de la sumatoria             N = ');

yverd = exp(x);   %valor calculado con la funcion intrinseca de octave
yaprx = 0;        %valor incial de la funcion aproximada por sumatoria

%% calculo principal
yaprx = 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) + x^11/factorial(11) + x^12/factorial(12) + ...
            x^13/factorial(13) + x^14/factorial(14) + x^15/factorial(15);

erro = abs(yverd - yaprx); %calulamos el error

%% escribimos los resultados
fprintf('yverd = %e\n',yverd);
fprintf('yaprx = %e\n',yaprx);
fprintf('error = %e\n',erro);