poly_roots(lpc_in, re, im, conv, 10, 20, 7);
poly_roots(lpc_in, re, im, conv, 10, 30, 7);
void lpc_flat(float gamma, float *lpc_in, float *lpc_out, int order){   int i;   float re[10], im[10], conv[10];   float alpha;   alpha = 1/(4-4*gamma);   poly_roots(lpc_in, re, im, conv, 10, 20, 7);   for (i=0;i<order;i++)   {      float fact,tmp;      float radius = sqrt(re[i]*re[i]+im[i]*im[i]);      tmp=1-radius;      if (tmp>2-2*gamma)         fact = tmp;      else         fact = alpha*tmp*tmp-gamma+1;      fact = (1-fact)/(radius+.001);      re[i]*=fact;      im[i]*=fact;   }   poly(re, im, lpc_out, order);}