From: Craig Schroeder [schroeder@santel.net]
Sent: Thursday, December 23, 1999 01:59
To: JoshS@santsys.com
Subject: Definite Integration.

/* This program performs definite integration.  It integrates
a given function, f(x), on the interval [a,b].  N can be set as
a variable to control precision.  The algorithm in use is Simpson's
Rule, and it is VERY fast. */

#include 
#include 

double f(double x) {return cos(x)+sin(x);}  // function to integrate.

double I(double (*f)(double), double a, double b, unsigned long N)
{
  double w=(b-a)/N, s=0, k=a;
  for(unsigned long i=0; i