/* 
 * arquivo: exp2epsilon.c
 * ----------------------
 * Dado x e epsilon (reais) e n, um inteiro, exp calcula exp(x) 
 * atraves da soma dos termos da serie 
 *
 * 1 + x + x^2/2! + x^3/3! + ...
 *
 * ate o valor absoluto da parcela se tornar <= epsilon.
 * Paramos quando encontramos um tal termo.
 */

#include <stdio.h>
#include <math.h>

int main()
{
  int i;
  double x, epsilon, s = 1.0, t = 1.0;

  printf("Forneca o valor de epsilon: ");
  scanf("%lf", &epsilon);
  printf("Forneca o valor de x: ");
  scanf("%lf", &x);

  i = 1;
  while (t > epsilon || t < -epsilon) {
    t *= x / i;
    s += t;
    i++;
  }

  printf("exp(%g) = %.20g [nossa serie]\n", x, s);
  printf("exp(%g) = %.20g [biblioteca math]\n", x, exp(x));
  return 0;
}