Changeset 2558 in ntrip for trunk/BNS/bnsutils.cpp
- Timestamp:
- Aug 10, 2010, 9:37:02 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNS/bnsutils.cpp
r1687 r2558 128 128 const ColumnVector& yi, // vector of the initial y-values 129 129 double dx, // the step size for the integration 130 ColumnVector (*der)(double x, const ColumnVector& y) 130 double* acc, // aditional acceleration 131 ColumnVector (*der)(double x, const ColumnVector& y, double* acc) 131 132 // A pointer to a function that computes the 132 133 // derivative of a function at a point (x,y) 133 134 ) { 134 135 135 ColumnVector k1 = der(xi , yi ) * dx;136 ColumnVector k2 = der(xi+dx/2.0, yi+k1/2.0 ) * dx;137 ColumnVector k3 = der(xi+dx/2.0, yi+k2/2.0 ) * dx;138 ColumnVector k4 = der(xi+dx , yi+k3 ) * dx;136 ColumnVector k1 = der(xi , yi , acc) * dx; 137 ColumnVector k2 = der(xi+dx/2.0, yi+k1/2.0, acc) * dx; 138 ColumnVector k3 = der(xi+dx/2.0, yi+k2/2.0, acc) * dx; 139 ColumnVector k4 = der(xi+dx , yi+k3 , acc) * dx; 139 140 140 141 ColumnVector yf = yi + k1/6.0 + k2/3.0 + k3/3.0 + k4/6.0;
Note:
See TracChangeset
for help on using the changeset viewer.