switch(this.integMethod){
// Trapezium rule
case 0: ncoef = Complex.oneDarray(2);
ncoef[0].reset(this.deltaT*this.cConst,0.0D);
ncoef[1].reset(this.deltaT*this.cConst,0.0D);
super.zNumer=new ComplexPoly(1);
super.zNumer.resetPoly(ncoef);
super.zNumerDeg=1;
dcoef = Complex.oneDarray(2);
dcoef[0].reset(this.bConst*this.deltaT - 2*this.aConst,0.0D);
dcoef[1].reset(this.bConst*this.deltaT + 2*this.aConst,0.0D);
super.zDenom=new ComplexPoly(1);
super.zDenom.resetPoly(dcoef);
super.zDenomDeg=1;
super.zZeros = Complex.oneDarray(1);
super.zZeros[0].reset(-1.0D, 0.0D);
super.zPoles = Complex.oneDarray(1);
super.zPoles[0].reset((2.0D*this.aConst-super.deltaT*this.bConst)/(2.0D*this.aConst+super.deltaT*this.bConst), 0.0D);
break;
// Backward rectangulr rule
case 1: ncoef = Complex.oneDarray(2);
ncoef[0].reset(0.0D,0.0D);
ncoef[1].reset(this.cConst*this.deltaT,0.0D);
super.zNumer=new ComplexPoly(1);
super.zNumer.resetPoly(ncoef);
super.zNumerDeg=1;
dcoef = Complex.oneDarray(2);
dcoef[0].reset(this.bConst*this.deltaT + this.aConst,0.0D);
dcoef[1].reset(this.aConst,0.0D);
super.zDenom=new ComplexPoly(1);
super.zDenom.resetPoly(dcoef);
super.zDenomDeg=1;
super.zZeros = Complex.oneDarray(1);
super.zZeros[0].reset(0.0D, 0.0D);
super.zPoles = Complex.oneDarray(1);
super.zPoles[0].reset(this.aConst/(super.deltaT*this.bConst+this.aConst), 0.0D);
break;
// Foreward rectangular rule
case 2: ncoef = Complex.oneDarray(1);
ncoef[0].reset(this.cConst*this.deltaT,0.0D);
super.zNumer=new ComplexPoly(0);
super.zNumer.resetPoly(ncoef);
super.zNumerDeg=0;
dcoef = Complex.oneDarray(2);
dcoef[0].reset(-this.aConst,0.0D);
dcoef[1].reset(this.bConst*this.deltaT - this.aConst,0.0D);
super.zDenom=new ComplexPoly(1);
super.zDenom.resetPoly(dcoef);
super.zDenomDeg=1;
super.zPoles = Complex.oneDarray(1);
super.zPoles[0].reset(this.aConst/(super.deltaT*this.bConst-this.aConst), 0.0D);
break;