9 double d_mod(x,y) doublereal *x, *y;
12 double drem(double, double);
17 double d_mod(doublereal *x, doublereal *y)
24 z = drem(xa = *x, ya);
34 if( (quotient = *x / *y) >= 0)
35 quotient = floor(quotient);
37 quotient = -floor(-quotient);
38 return(*x - (*y) * quotient );