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