1 /* 2 * "@(#)r_mod.c 1.1" 3 */ 4 5 double r_mod(x,y) 6 float *x, *y; 7 { 8 double floor(), quotient = *x / *y; 9 if (quotient >= 0.0) 10 quotient = floor(quotient); 11 else 12 quotient = -floor(-quotient); 13 return(*x - (*y) * quotient ); 14 } 15