#include "f2c.h" #undef abs #include double pow_dd (doublereal * ap, doublereal * bp) { return (pow (*ap, *bp)); }