Cdouble pj = poles[jp]; Cdouble hi = hi(pj,poles,zeros,gain); Cdouble hj = hr(pj,poles,zeros,gain); Cdouble hihj = hi.times(hj); double fb0,fb1,fb2,rb0,rb1,rb2,b0,b1,b2,a1,a2;
Cdouble pj = poles[jp]; Cdouble hi = hi(pj,poles,zeros,gain); Cdouble hj = hr(pj,poles,zeros,gain); Cdouble hihj = hi.times(hj); double fb0,fb1,fb2,rb0,rb1,rb2,b0,b1,b2,a1,a2;