%prog019ExponencialV03

clearvars;
close all;
clc;

%en este caso usaremos un ciclo for para tener la sumatoria y sin factorial
%sin el elevado a (^)

%% datos inciales

x = input('x = ');  %valor de x
N = input('N = ');  %numero de elementos de la sumatoria
yverd = cos(x);

%% calculo principal con ciclo for para la sumatoria

yaprox = 1;                 %valor incial aproximado
errr = abs(yverd - yaprox); %error inicial

fprintf('yverd   = %f --- yaprox = %f\n',yverd,yaprox);
fprintf('errr    = %f\n',errr);

%para calcular sumatoria
for n = 1:N%--------------------------------------------->inicio ciclo for externo
    
    %para calcular factorial
    fakto = 1;
    for k = 1:(2*n)%------------->inicio for interno
        fakto = fakto*k;
    end;%-------------------->fin for interno

    %para calcular el elevado
    elevado = 1;
    for k = 1:(2*n)%----------->inicio for interno
        elevado = elevado*x;
    end;%-------------------->fin for interno
     
    yaprox = yaprox + (-1)^n*elevado/fakto;
    
    %error
    errr = abs(yverd - yaprox);
    
    %escribe resultados paso a paso
    disp(' ')
    fprintf('n       = %f\n ',n);
    fprintf('elevado = %f\n ',elevado);
    fprintf('fakto   = %f\n',fakto);
    fprintf('yverd   = %f --- yaprox = %f\n',yverd,yaprox);
    fprintf('errr    = %f\n',errr);
    
    
end;%--------------------------------------------------->fin del for externo
  
%% error

errr = abs(yverd - yaprox);
           
%% resultados           
    disp(' ')
    fprintf('n       = %d\n ',n);
    fprintf('elevado = %f\n ',elevado);
    fprintf('fakto   = %f\n',fakto);
    fprintf('yverd   = %f --- yaprox = %f\n',yverd,yaprox);
    fprintf('errr    = %f\n',errr);