#include <ncbi.h>
#include <gishlib.h>

/* 
integer power function

Originally submission by John Spouge, 6/25/90
Added to shared library by WRG
*/
double
fct_powi(x, n)
	double	x;				/* argument */ 
	register int	n;		/* power */
{
	register int	i;
	double	y;

	y = 1.;
	for (i = ABS(n); i > 0; i /= 2) {
		if (i&1)
			y *= x;
		x *= x;
	}
	if (n >= 0)
		return y;
	return 1./y;
}