%%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

    elevado = 1;
   for i=1:n %ciclo for para calcular la potencia necesaria -----> inicio ciclo for interno
     elevado = elevado * x;
   end; %------------------> fin ciclo for interno

  terminoActual = elevado/factorial(n); %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