/* demean, detrend, de t^2, ... de t^n */


float *farray1();
float Legpolyn();

denthpow(data,ndata,n)
float *data;
int ndata,n;
{
int i,j;
float a,*p,x,x0,slop;
  
  p=farray1(0,ndata-1);

  x0=-1.0+1.0/(float)ndata;
  slop=2.0/(float)ndata;
  for(i=0;i<=n;i++)
  { a=0.0;
    for(j=0;j<ndata;j++)
    { x=x0+slop*(float)j;
      p[j]=Legpolyn(i,x);
      a+=data[j]*p[j];
    }
    a*=slop;
    for(j=0;j<ndata;j++) data[j]-=a*p[j];
  }
  free_farray1(p,0);
  
}