home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c072 / 1.ddi / PRG3_3A.C < prev    next >
Encoding:
C/C++ Source or Header  |  1987-09-19  |  2.1 KB  |  73 lines
  1. /*Program 3_3a -- Integration of a user function 'func'
  2.    by Stephen R. Davis, 1987
  3.  
  4.   This function integrates a user function 'func' using the trapezoid
  5.   rule to evaluate a definite integral over the range 'a' to 'b' by
  6.   dividing it into 'steps' number of discrete intervals.
  7. */
  8.  
  9. #include <stdio.h>
  10.  
  11. /*prototype definitions --*/
  12. double integrate (double, double, unsigned);
  13. double func (double);
  14. double answer (double);
  15. int main (void);
  16.  
  17. /*Integrate - evaluate the definite integral of a user function
  18.               'func' over the interval 'a' to 'b' by using the
  19.               trapazoid rule on 'steps' number of subintervals*/
  20. double integrate (a, b, steps)
  21.     double a,b;
  22.     unsigned steps;
  23. {
  24.     double func ();
  25.     unsigned i;
  26.     double integral, spacing, x;
  27.  
  28.     spacing = (b - a) / steps;         /*divide the interval evenly*/
  29.     integral = (func(b) + func(a)) / 2.;
  30.     for (i = 1,x = a; i < steps; i++) { /*accumulate value over...*/
  31.          x += spacing;                  /*...the interval*/
  32.          integral += func(x);
  33.     }
  34.     return (integral * spacing);
  35. }
  36.  
  37. /*Func - a test function to integrate and its value.
  38.          (the integral of x**2 is equal to (x**3)/3).  Replace
  39.          func() and answer() with any desired function/integral
  40.          pair and recompile.*/
  41. double func (x)
  42.     double x;
  43. {
  44.     return (x*x);
  45. }
  46. double answer (x)
  47.     double x;
  48. {
  49.     return x*x*x/3.;
  50. }
  51.  
  52. /*Main - test routine integrate() w/ 'func'at various starting
  53.          value, stoping value and step size.  Compare with 'value'
  54.          for accuracy.*/
  55. main ()
  56. {
  57.     float a,b;
  58.     unsigned steps;
  59.  
  60.     printf ("Enter starting x, ending x and number of steps\n");
  61.     printf ("   (exit by entering starting x equal ending x)\n");
  62.     for (;;) {
  63.          printf (">");
  64.          scanf ("%f %f %d", &a, &b, &steps);
  65.      printf ("a = %f, b = %f, steps = %d\n", a, b, steps);
  66.          if (a == b)
  67.               break;
  68.          printf ("Integral is %f\n", integrate((double)a,
  69.                                                (double)b, steps));
  70.          printf ("(actual value is %f)\n", answer(b) - answer(a));
  71.     }
  72. }
  73.