190   p->b0 = p->b1 = p->b2 = p->a1 = p->a2 = 0;  in start()
220       p->b2 =  (1 - cos(w0))/2;  in start()
229       p->b2 =  (1 + cos(w0))/2;  in start()
238       p->b2 =  -sin(w0)/2;  in start()
247       p->b2 =  -alpha;  in start()
256       p->b2 =   1;  in start()
265       p->b2 =   1 + alpha;  in start()
276       p->b2 =   1 - alpha*A;  in start()
287       p->b2 =    A*( (A+1) - (A-1)*cos(w0) - 2*sqrt(A)*alpha );  in start()
300       p->b2 =    A*( (A+1) + (A-1)*cos(w0) - 2*sqrt(A)*alpha );  in start()
337       p->b2 = 1 + sin(w0);  in start()
380         double b_re = p->b0 + p->b1 * cos(-y) + p->b2 * cos(-2 * y);  in start()
382         double b_im = p->b1 * sin(-y) + p->b2 * sin(-2 * y);  in start()
385         p->b0 *= g; p->b1 *= g; p->b2 *= g;  in start()
387       mult = (p->b0 + p->b1 + p->b2) / (p->a0 + p->a1 + p->a2);  in start()