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
);