#include "f2c.h" #define log10e 0.43429448190325182765 #undef abs #include double d_lg10 (doublereal * x) { return (log10e * log (*x)); }