This is solut.c in view mode; [Download] [Up]
#include "lpc.h" solut(b, r, th) double b[], r[], th[]; { int j, kerr; double coe[128], rootr[128], rooti[128]; for (j=0; j < NPP1 ;++j) coe[j] = b[NPP1-(j+1)]; /* coe[j] = b[j]; */ kerr = 1; dproot(coe, rootr, rooti, kerr); for (j=0; j < NPOLE ;++j) { r[j]=sqrt(pow(rootr[j],(double)2.0)+pow(rooti[j],(double)2.0)); th[j]=atan2(rooti[j], rootr[j]); printf(" r and th %e %e\n",r[j],th[j]); } return; }
These are the contents of the former NiCE NeXT User Group NeXTSTEP/OpenStep software archive, currently hosted by Netfuture.ch.