home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / ADMBAS.C next >
Encoding:
Text File  |  1984-07-17  |  641 b   |  34 lines
  1.    admbas(n,h,f,b)
  2.  
  3.       /*purpose...integral of given function by Adams-Bashforth.*/
  4.  
  5.      int n;
  6.      float h,f[],*b;
  7.  
  8.     {
  9.  
  10.      int i,n3;
  11.      float sum;
  12.  
  13.       sum = 0.;
  14.  
  15.       if(n == 6) goto ne6;
  16.       if(n < 6) goto nl6;
  17.  
  18.        n3 = n - 4;
  19.        for(i = 3; i <= n3 ; i++)
  20.         sum += f[i];
  21.        sum *= 24.;
  22.  
  23. ne6:  sum -= 31. * f[2] + 33. * f[n-3];
  24.       goto nne6;
  25. nl6:  if(n < 5) goto nl5;
  26.       sum -= 22. * f[2];
  27. nne6: sum += 28. * f[1] - 4. * f[n-2];
  28.       goto nge5;
  29. nl5:  sum = sum + 37. * f[1] - 59. * f[2];
  30. nge5: sum = sum - 9. * f[0] + 55. * f[n-1];
  31.       *b = (sum*h)/24.;
  32.     }
  33.  
  34.