#include "f2c.h" #undef abs #include double r_nint (real * x) { return ((*x) >= 0 ? floor (*x + .5) : -floor (.5 - *x)); }