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

#define EPSILON 1e-10

int main(void)
{
  int n;
  double rad, ckrad, sinrad = 0.0;
  double k = 1.0, l, radrad;

  printf("input rad =");
  scanf("%lf", &rad);

  ckrad = l = rad;
  radrad = rad * rad;
  
  for (n = 0; n < 5 && fabs(ckrad) >= EPSILON * sinrad; n++) {
    if ((n % 2) == 0)
      sinrad += ckrad;
    else
      sinrad -= ckrad;
    k *= (2 * k + 1) * (2 * (k + 1) + 1);
    l *= radrad;
    ckrad = l / k;
  }

  
  printf("sin(rad)        = %.15lf\n", sin(rad));
  printf("taylor sin(rad) = %.15lf\n", sinrad);


  return 0;
}