#include "f2c.h" #undef abs #include double d_sqrt (doublereal * x) { return (sqrt (*x)); }