%%prog20230419 funcion exponencial
clearvars;
clear all;
close all;
clc;

% el valor de la funcion exp(x) puede ser calculado de la siguiente forma
% y = 1 + x^1/factorial(1) + x^2/factorial(2) + x^3/factorial(3) + ... + x^N/factorial(N)

x = input('x = ');
N = input('N = ');
yverdadero = exp(x);

sumatoria = 1;
for n=1:N % --------------------> inicio ciclo for externo

  facto = 1;
  for i=1:n %--------------------> inicio ciclo for interno
    facto = facto * i;
  end;%---------------------------> fin ciclo for interno

  terminoActual = x^n/facto; %realizamos el cálculo del termino actual
  sumatoria = sumatoria + terminoActual; %realizamos a la variable acumulativa
  error = abs(yverdadero - sumatoria); %acá calculamos el error entre el valor calculado y el valor verdadero
                                       %ocupando la función exp(x) incluida en octave
  fprintf('x = %f, n = %f, expCalculado = %f expVerdadero = %f,error = %e\n',x,n,sumatoria,yverdadero,error)

  %pause(0.5);
end; %--------------> fin ciclo for externo