This is play.c in view mode; [Download] [Up]
play(p,n_args)
float *p;
{
int jj;
float input[4];
double sig[500],errn,rms1,rms2,cc[40];
long i;
setnot(p[0],p[1],0);
for(i = 0; i<500; i++) {
getin(input,0);
sig[i] = (double)input[0];
}
alpol(sig,&errn,&rms1,&rms2,cc);
printf(" %e %e %e %e %e %e \n",errn,rms1,rms2,cc[0],cc[1],cc[10]);
for(jj=0; jj<20; jj++) printf(" coefs %f %d\n",cc[jj],jj);
jj = stable(cc);
printf(" jj = %d\n",jj);
/* if(jj < 0) */ fixit(cc);
for(jj=0; jj<20; jj++) printf(" coefs %e %d\n",cc[jj],jj);
printf(" %e %e %e %e %e %e \n",errn,rms1,rms2,cc[0],cc[1],cc[10]);
/* if(jj < 0) */ fixit(cc);
for(jj=0; jj<20; jj++) printf(" coefs %e %d\n",cc[jj],jj);
printf(" %e %e %e %e %e %e \n",errn,rms1,rms2,cc[0],cc[1],cc[10]);
}
These are the contents of the former NiCE NeXT User Group NeXTSTEP/OpenStep software archive, currently hosted by Netfuture.ch.