added setsid man page.

[unix-history] / lib / libF77 / pow_di.c

#include "f2c.h"

#ifdef KR_headers
double pow_di(ap, bp) doublereal *ap; integer *bp;
#else
double pow_di(doublereal *ap, integer *bp)
#endif
{
double pow, x;
integer n;

pow = 1;
x = *ap;
n = *bp;

if(n != 0)
       {
       if(n < 0)
               {
               n = -n;
               x = 1/x;
               }
       for( ; ; )
               {
               if(n & 01)
                       pow *= x;
               if(n >>= 1)
                       x *= x;
               else
                       break;
               }
       }
return(pow);
}