allpole filter…
This runs an allpole filter. Good luck.
From the source code:
float allpole(float x, int *jcount, int npoles, float *past, float *c)