float breson(x,a,out,nvals)
float *x,*a,*out;
{
	int i;
	for(i=0; i<nvals; i++) {
		out[i] = *a * x[i] + *(a+1) * *(a+3) - *(a+2) * *(a+4);
		*(a+4) = *(a+3);
		*(a+3) = out[i];
	}
}