home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c072 / 1.ddi / PRG3_3B.C < prev    next >
Encoding:
C/C++ Source or Header  |  1987-09-19  |  2.1 KB  |  72 lines
  1. /*Program 3_3b -- Integration of a user function using Trapezoid Rule
  2.   by Stephen R. Davis
  3.  
  4.   This program demostrates how by allowing a function to be passed
  5.   to another function, C allows some very general purpose routines
  6.   to be designed.  This is a more general routine than that of 3_3a
  7.   since it makes no restrictions on the name of the user function.
  8.   This is the type of function which could be included in a C library
  9.   for later use.
  10. */
  11.  
  12. #include <stdio.h>
  13.  
  14. /*prototype definitions --*/
  15. double integrate (double (*)(double), double, double, unsigned);
  16. double func (double), answer (double);
  17. int main (void);
  18.  
  19. /*Integrate - evaluate the definite integral of the function pointed at
  20.               by 'fn' by applying the trapezoid rule to the function
  21.               over the range 'a' to 'b' by dividing it into 'steps'
  22.               number of intervals.*/
  23. double integrate (fn, a, b, steps)
  24.     double (*fn)(),a,b;
  25.     unsigned steps;
  26. {
  27.     unsigned i;
  28.     double integral,spacing,x;
  29.  
  30.     spacing = (b - a) / steps;
  31.     integral = ((*fn)(b) + (*fn)(a)) / 2;
  32.     for (i = 1,x = a; i < steps; i++) {
  33.          x += spacing;
  34.          integral += (*fn)(x);
  35.     }
  36.     return (integral * spacing);
  37. }
  38.  
  39. /*Experiment - a sample user program to integrate.  (It has renamed
  40.                from 'func' to avoid any confusion with the pointer
  41.                to a function 'fn' in 'Integrate')*/
  42. double experiment (x)
  43.     double x;
  44. {
  45.     return (x*x);
  46. }
  47. double answer (x)
  48.     double x;
  49. {
  50.     return x*x*x/3.;
  51. }
  52.  
  53. /*Main - same as before*/
  54. main ()
  55. {
  56.     float a,b;
  57.     unsigned steps;
  58.  
  59.     printf ("Enter starting x, ending x and number of steps\n");
  60.     printf ("    (exit by entering starting x equal ending x)\n");
  61.     for (;;) {
  62.          printf (">");
  63.          scanf ("%f %f %d", &a, &b, &steps);
  64.      printf ("a = %f, b = %f, steps = %d\n", a, b, steps);
  65.          if (a == b)
  66.               break;
  67.          printf ("Integral is %f\n", integrate(experiment,
  68.                           (double)a, (double)b, steps));
  69.          printf ("(actual value is %f)\n", answer (b) - answer (a));
  70.     }
  71. }
  72.