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