home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / INTEGS.C < prev    next >
Encoding:
Text File  |  1984-06-28  |  550 b   |  24 lines
  1.    integs(m,y,xlow,xhi,ans)
  2.  
  3.       /*this function evaluates the definite integral of a */
  4.       /*function by Simpson's one-third rule.              */
  5.  
  6.       int m;
  7.       float y[],xlow,xhi,*ans;
  8.  
  9.     {
  10.       int xm,mm4,i;
  11.  
  12.       xm = m;
  13.       mm4 = m - 4;
  14.       *ans = 0.0;
  15.       if (m <= 1) return;
  16.       *ans = y[0] + 4.0*y[m - 2] + y[m-1];
  17.       if(m != 3)
  18.       {
  19.        for(i = 1; i <= mm4; i += 2)
  20.         *ans = *ans + 4.0*y[i] + y[i + 1] + y[i + 1];
  21.       }
  22.       *ans = *ans * (xhi - xlow)/(3.0*(xm - 1.0));
  23.     }
  24.