1 2 3 4 5 6 7 8 9 10 11 12
#include "f2c.h" #ifdef KR_headers double d_abs(x) doublereal *x; #else double d_abs(doublereal *x) #endif { if(*x >= 0) return(*x); return(- *x); }