#función que calcula e^-x de 3 formas diferentes
#consideramos el verdadero valor el de la función del lenguaje

error_e_x = function(x,E) {
  S1=0 #acumulador de la forma de division
  S2=0 #acumulador de la forma alternante
  tabla=c(0,exp(-x),0,S2,exp(-x)-0,exp(-x)-S2)
  n=0
  while (abs(exp(x)-S1) >E || abs(exp(-x)-S2) >E) {	    
    S1 = S1 + (x^n)/factorial(n)
    S2 = S2 + ((-1)^n * x^n)/factorial(n)
    n  = n + 1
    tabla=rbind(tabla,c(n,exp(-x),1/S1,S2,exp(-x)-1/S1,exp(-x)-S2))    
  }#while
  tabla
}

#para probar la función
error_e_x(1,0.0001)