void fpoly(float x, float p[], int np)
{
	int j;

	p[1]=1.0;
	for (j=2;j<=np;j++) p[j]=p[j-1]*x;
}
/* (C) Copr. 1986-92 Numerical Recipes Software $!6)D$|2$3. */