Facebook
From Gamboge Moth, 7 Years ago, written in C.
Embed
Download Paste or View Raw
Hits: 262
  1. double rk4( double x0, double y0, double h, double (*fun)(double, double))
  2. {
  3.         double y1;
  4.         double k1,k2,k3,k4;
  5.         k1 = h*fun(x0     ,y0      );
  6.         k2 = h*fun(x0+h/2.,y0+k1/2.);
  7.         k3 = h*fun(x0+h/2.,y0+k2/2.);
  8.         k4 = h*fun(x0+h   ,y0+k3   );
  9.         y1 = y0 + (k1+2.*k2+2.*k3+k4)/6.;
  10.         return y1;
  11. }