remainder function

double remainder(double x, double y)

Implementation

double remainder(double x, double y) {
  double n = rint(x / y);
  return x - n * y;
}