uses MtxExpr, Math387, MtxVec, SignalUtils, MtxVecTee, MtxVecEdit, IirFilters,
       LinearSystems;

  procedure TForm1.Button1Click(Sender: TObject);
  var z,p, num,den, FreqFr,Response: Vector;
      Order: integer;
      k,Bw,Wc: double;
      WcArray: TDoubleArray; //modified 3dB frequency
  begin
      SetLength(WcArray,2);
      Order := ChebyshevIIOrder([1,3,6,9],0.2,50,ftBandpass,WcArray,True);
      ChebyshevIIAnalog(Order,50,z,p,k);  //design analog protype
      Wc := Sqrt(WcArray[0]*WcArray[1]); //modified 3dB frequency
      Bw := WcArray[1]-WcArray[0];
      LowpassToBandpass(z,p,k,Wc,Bw);  //frequency transformation in s-domain
      ZeroPoleToTransferFun(num,den,z,p,k);
      FreqFr.Length := 1000;         //Define the frequency grid (logarithmic)
      LogRamp(FreqFr,-1,1); //between 0.1 (=10^(-1)) and 10 (=10^1) rad/sec
      FrequencyResponseS(num,den,FreqFr,Response); //Laplace
      DrawIt(Response); //Y axis linear, X axis logarithmic;
  end;

  #include "MtxExpr.hpp"
  #include "MtxVecEdit.hpp"
  #include "MtxVecTee.hpp"
  #include "SignalUtils.hpp"
  #include "IirFilters.hpp"
  #include "LinearSystems.hpp"

  void __fastcall TForm1::BitBtn1Click(TObject *Sender)
  {
      sVector z,p, num,den, FreqFr, Response, WcArray;
      int Order;
      double k,Wc,Bw;

      WcArray.Size(2);
      Order = ChebyshevIIOrder(OPENARRAY(double,(1,3,6,9)),0.2,50,ftBandpass,WcArray.PValues1D(0),WcArray.Length-1,true);
      ChebyshevIIAnalog(Order,50,z,p,k);  //design analog protype
      Wc = Sqrt(WcArray[0] * WcArray[1]); //modified 3dB frequency
      Bw = WcArray[1] - WcArray[0];
      LowpassToBandpass(z,p,k,Wc,Bw);  //frequency transformation in s-domain
      ZeroPoleToTransferFun(num,den,z,p,k);
      FreqFr.Length = 1000;         //Define the frequency grid (logarithmic)
      LogRamp(FreqFr,-1,1); //between 0.1 (=10^(-1)) and 10 (=10^1) rad/sec
      FrequencyResponseS(num,den,FreqFr,Response); //Laplace
      DrawIt(Response); //Y axis linear, X axis logarithmic;
  }