#include "lpc.h"
fixit(coef)
double coef[];
{
        double r[50], th[50], b[50];
        int m;

        b[0] = (double) 1.0;
        for (m=0; m < NPOLE ;++m)  
                b[m + 1] = coef[m];
        solut(b, r, th);
        corect(r, th, coef);
        return;
}