for (int i1=0,i2=0,jp=0; i1<_n1; ++jp) { Cdouble pj = poles[jp]; Cdouble hi = hi(pj,poles,zeros,gain); Cdouble hj = hr(pj,poles,zeros,gain); Cdouble hihj = hi.times(hj);
for (int i1=0,i2=0,jp=0; i1<_n1; ++jp) { Cdouble pj = poles[jp]; Cdouble hi = hi(pj,poles,zeros,gain); Cdouble hj = hr(pj,poles,zeros,gain); Cdouble hihj = hi.times(hj);