home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / DEFOAD.C < prev    next >
Encoding:
Text File  |  1984-08-07  |  629 b   |  27 lines
  1.    defoad(iter,h,x,y,dy,iflag)
  2.  
  3.       /* subroutine to solve first-order ordinary differential */
  4.       /* equations y'=f(x,y) by Adam's method.                 */
  5.  
  6.       int iflag,iter;
  7.       float dy[],h,x[],y[];
  8.     {
  9.       int i,j;
  10.       float xn,yn;
  11.       double yprime();
  12.  
  13.       xn = x[0];
  14.       yn = y[0];
  15.       dy[0] = yprime(xn,yn,iflag);
  16.       j = 4;
  17.       defork(j,h,x,y,dy,iflag);
  18.  
  19.       for(i = 3; i <= iter-2; i++)
  20.       {
  21.        y[i+1]=y[i]+h*(55.*dy[i]-59.*dy[i-1]+37.*dy[i-2]-9.*dy[i-3])/24.;
  22.        xn = x[i+1];
  23.        yn = y[i+1];
  24.        dy[i+1] = yprime(xn,yn,iflag);
  25.       }
  26.     }
  27.