#include "f2c.h"

#ifdef KR_headers

double pow();

double pow_dd(ap, bp) doublereal *ap, *bp;

#else

#undef abs

#include "math.h"

double pow_dd(doublereal *ap, doublereal *bp)

#endif

{

return(pow(*ap, *bp) );

}