%% prog20230510 exponencial while 01

clearvars;
clear all;
close all;
clc;

%este prpograma calcula el valor de la funcion exponencial exp(x) usando una sumatoria
%y al mismo tiempo determina el numero de elementos de la sumatoria N que son
%necesarios para calcular el valor de la funcion con un error maximo de una
%millonesima 1e-6


%% lectura de datos
x   = input('numero al que se calculara la funcion exponencial x = ');
tol = input('error maximo que se tolera                        tol = ');

yverd = exp(x);
yaprx = 1;
N = 1;

%calculamos el error de forma previa
erro = abs(yverd - yaprx);

%mientras el error sea mas grande que la tolerancia
%el ciclo while se repetira
%cuando el error sea menor que la tolerancia se deja de repetir
while (erro >= tol) && (N <= 100)

      yaprx = 1;

      for n = 1:N
          yaprx = yaprx + x^n/factorial(n);
      end;

      erro = abs(yverd - yaprx);
      fprintf('N = %d --> yverd = %e --> yaprx = %e --> erro = %e\n',N,yverd,yaprx,erro);
      N = N + 1;
end;

disp(' ');
disp('resultados finales');
fprintf('N = %d --> yverd = %e --> yaprx = %e --> erro = %e\n',N-1,yverd,yaprx,erro);