math/test-math.c

   1 #include <ansidecl.h>
   2 #include <stdio.h>
   3 #include <math.h>
   4 #include <stdlib.h>
   5
   6 int
   7 DEFUN_VOID(main)
   8 {
   9   CONST char str[] = "123.456";
  10   double x,h,li,lr,a,lrr;
  11
  12   x = atof (str);
  13
  14   printf ("%g %g\n", x, pow (10.0, 3.0));
  15
  16   x = sinh(2.0);
  17
  18   printf("sinh(2.0) = %g\n", x);
  19
  20   x = sinh(3.0);
  21
  22   printf("sinh(3.0) = %g\n", x);
  23
  24   h = hypot(2.0,3.0);
  25
  26   printf("h=%g\n", h);
  27
  28   a = atan2(3.0, 2.0);
  29
  30   printf("atan2(3,2) = %g\n", a);
  31
  32   lr = pow(h,4.0);
  33
  34   printf("pow(%g,4.0) = %g\n", h, lr);
  35
  36   lrr = lr;
  37
  38   li = 4.0 * a;
  39
  40   lr = lr / exp(a*5.0);
  41
  42   printf("%g / exp(%g * 5) = %g\n", lrr, a, lr);
  43
  44   lrr = li;
  45
  46   li += 5.0 * log(h);
  47
  48   printf("%g + 5*log(%g) = %g\n", lrr, h, li);
  49
  50   printf("cos(%g) = %g,  sin(%g) = %g\n", li, cos(li), li, sin(li));
  51
  52   x = drem(10.3435,6.2831852);
  53
  54   printf("drem(10.3435,6.2831852) = %g\n", x);
  55
  56   x = drem(-10.3435,6.2831852);
  57
  58   printf("drem(-10.3435,6.2831852) = %g\n", x);
  59
  60   x = drem(-10.3435,-6.2831852);
  61
  62   printf("drem(-10.3435,-6.2831852) = %g\n", x);
  63
  64   x = drem(10.3435,-6.2831852);
  65
  66   printf("drem(10.3435,-6.2831852) = %g\n", x);
  67
  68
  69   printf("x%8.6gx\n", .5);
  70   printf("x%-8.6gx\n", .5);
  71   printf("x%6.6gx\n", .5);
  72
  73   {
  74     double x = atof ("-1e-17-");
  75     printf ("%g %c= %g %s!\n",
  76             x,
  77             x == -1e-17 ? '=' : '!',
  78             -1e-17,
  79             x == -1e-17 ? "Worked" : "Failed");
  80   }
  81
  82   return 0;
  83 }