home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / DETOAD.C < prev    next >
Encoding:
Text File  |  1984-08-07  |  1022 b   |  33 lines
  1.    detoad(iter,h,x,y,dy,ddy,iflag)
  2.  
  3.       /* subroutine to solve third-order ordinary differential  */
  4.       /* equations y''' = f(x,y,y',y'') by Adams' method.       */
  5.  
  6.       int iflag,iter;
  7.       float dy[],ddy[],h,x[],y[];
  8.  
  9.     {
  10.       int i,lda;
  11.       float u1,u2,u3,u4;
  12.       double ytrpri();
  13.  
  14.       lda = 4;
  15.       detork(lda,h,x,y,dy,ddy,iflag);
  16.       u1 = ytrpri(x[0],y[0],dy[0],ddy[0],iflag);
  17.       u2 = ytrpri(x[1],y[1],dy[1],ddy[1],iflag);
  18.       u3 = ytrpri(x[2],y[2],dy[2],ddy[2],iflag);
  19.       u4 = ytrpri(x[3],y[3],dy[3],ddy[3],iflag);
  20.  
  21.       for(i = 3; i <= iter-2; i++)
  22.       {
  23.        ddy[i+1] = ddy[i] + h * (55.*u4-59.*u3 + 37.*u2-9.*u1)/24.;
  24.        dy[i+1] = dy[i]+h*(55.*ddy[i]-59.*ddy[i-1]+37.*ddy[i-2]-9.*
  25.                                                     ddy[i-3])/24.;
  26.        y[i+1] = y[i]+h*(55.*dy[i]-59.*dy[i-1]+37.*dy[i-2]-9.*dy[i-3])/24.;
  27.        u1 = u2;
  28.        u2 = u3;
  29.        u3 = u4;
  30.        u4 = ytrpri(x[i+1],y[i+1],dy[i+1],ddy[i+1],iflag);
  31.       }
  32.     }
  33.