%%
%% This is file `bodeplot.sty',
%% generated with the docstrip utility.
%%
%% The original source files were:
%%
%% bodeplot.dtx  (with options: `package')
%% This is a generated file.
%% Copyright (C) 2021-2024 by Rushikesh Kamalapurkar
%% This file may be distributed and/or modified under the
%% conditions of the LaTeX Project Public License, either
%% version 1.3c of this license or (at your option) any later
%% version. The latest version of this license is in:
%% http://www.latex-project.org/lppl.txt
%% and version 1.3c or later is part of all distributions of
%% LaTeX version 2006/05/20 or later.
\NeedsTeXFormat{LaTeX2e}[2018/04/01]
\DeclareRelease{}{2024-02-06}{bodeplot-2024-02-06.sty}
\DeclareCurrentRelease{}{2024/04/18}
\ProvidesPackage{bodeplot}[2024/04/18 v1.2 Generate Bode, Nichols, and Nyquist plots]
\RequirePackage{pdftexcmds}
\RequirePackage{ifplatform}
\RequirePackage{xparse}
\RequirePackage{pgfplots}
  \pgfplotsset{compat=1.18}
  \usepgfplotslibrary{groupplots}

\newif\if@pgfarg\@pgfargfalse
\DeclareOption{pgf}{
  \@pgfargtrue
}
\newif\if@declutterarg\@declutterargfalse
\DeclareOption{declutter}{
  \@declutterargtrue
}
\newif\if@radarg\@radargfalse
\DeclareOption{rad}{
  \@radargtrue
}
\newif\if@hzarg\@hzargfalse
\DeclareOption{Hz}{
  \@hzargtrue
}
\ProcessOptions\relax
\newcommand{\n@mod}[2]{(#1)-((round((#1)/(#2)))*(#2))}
\newcommand{\n@mod@p}[2]{(#1)-((floor((#1)/(#2)))*(#2))}
\newcommand{\n@mod@n}[2]{(#1)-((floor((#1)/(#2))+1)*(#2))}
\if@pgfarg
  \newcommand{\n@pow}[2]{(#1)^(#2)}
\else
  \newcommand{\n@pow}[2]{(#1)**(#2)}
  \newcounter{gnuplot@id}
  \setcounter{gnuplot@id}{0}
  \if@declutterarg
    \edef\bodeplot@prefix{gnuplot/\jobname}
  \else
    \edef\bodeplot@prefix{\jobname}
  \fi
  \tikzset{
    gnuplot@prefix/.style={
      id=\arabic{gnuplot@id},
      prefix=\bodeplot@prefix
    }
  }
  \ifwindows\else
    \if@declutterarg
      \immediate\write18{mkdir -p gnuplot}
    \fi
  \fi
\fi
\newif\if@babel\@babelfalse
\AtBeginDocument{%
  \@ifpackageloaded{babel}{%
    \@babeltrue
    \let\shorthand@list\@empty
    \def\do#1{%
      \begingroup
        \lccode`\~=`#1\relax
        \lowercase{\ifbabelshorthand~{\g@addto@macro\shorthand@list{~}}{}}
      \endgroup
    }
    \dospecials
  }{}
}
\pgfplotsset{
  bode@style/.style = {
    label style={font=\footnotesize},
    tick label style={font=\footnotesize},
    grid=both,
    major grid style={color=gray!80},
    minor grid style={color=gray!20},
    x label style={at={(ticklabel cs:0.5)},anchor=near ticklabel},
    y label style={at={(ticklabel cs:0.5)},anchor=near ticklabel},
    scale only axis,
    samples=200,
    width=5cm,
    log basis x=10
  }
}
\pgfplotsset{freq@filter/.style = {}}
\def\freq@scale{1}
\pgfplotsset{freq@label/.style = {xlabel = {Frequency (rad/s)}}}
\pgfplotsset{ph@x@label/.style = {xlabel={Phase (deg)}}}
\pgfplotsset{ph@y@label/.style = {ylabel={Phase (deg)}}}
\def\ph@scale{180/pi}
\if@radarg
  \pgfplotsset{ph@y@label/.style = {ylabel={Phase (rad)}}}
  \pgfplotsset{ph@x@label/.style = {xlabel={Phase (rad)}}}
  \def\ph@scale{1}
\fi
\if@hzarg
  \def\freq@scale{2*pi}
  \pgfplotsset{freq@label/.style = {xlabel = {Frequency (Hz)}}}
  \if@pgfarg
    \pgfplotsset{freq@filter/.style = {x filter/.expression={x-log10(2*pi)}}}
  \fi
\fi
\tikzset{
  phase unit/.initial={deg},
  phase unit/.default={deg},
  phase unit/.is choice,
  phase unit/deg/.code={
    \renewcommand{\ph@scale}{180/pi}
    \pgfplotsset{ph@x@label/.style = {xlabel={Phase (deg)}}}
    \pgfplotsset{ph@y@label/.style = {ylabel={Phase (deg)}}}
  },
  phase unit/rad/.code={
    \renewcommand{\ph@scale}{1}
    \pgfplotsset{ph@y@label/.style = {ylabel={Phase (rad)}}}
    \pgfplotsset{ph@x@label/.style = {xlabel={Phase (rad)}}}
  },
  frequency unit/.initial={rad},
  frequency unit/.default={rad},
  frequency unit/.is choice,
  frequency unit/Hz/.code={
    \renewcommand{\freq@scale}{2*pi}
    \pgfplotsset{freq@label/.style = {xlabel = {Frequency (Hz)}}}
    \if@pgfarg
      \pgfplotsset{freq@filter/.style = {x filter/.expression={x-log10(2*pi)}}}
    \fi
  },
  frequency unit/rad/.code={
    \renewcommand{\freq@scale}{1}
    \pgfplotsset{freq@label/.style = {xlabel = {Frequency (rad/s)}}}
  }
}
\def\get@interval@start#1:#2\@nil{#1}
\def\get@interval@end#1:#2\@nil{#2}
\newcommand*{\MagK}[2]{(20*log10(abs(#1)))}
\newcommand*{\MagKAsymp}{\MagK}
\newcommand*{\MagKLin}{\MagK}
\newcommand*{\PhK}[2]{((#1<0?-pi:0)*\ph@scale)}
\newcommand*{\PhKAsymp}{\PhK}
\newcommand*{\PhKLin}{\PhK}
\newcommand*{\MagDel}[2]{0}
\newcommand*{\PhDel}[2]{(-#1*t*\ph@scale)}
\newcommand*{\MagPole}[2]
  {(-20*log10(sqrt(\n@pow{#1}{2} + \n@pow{t - (#2)}{2})))}
\newcommand*{\MagPoleLin}[2]{(t < sqrt(\n@pow{#1}{2} + \n@pow{#2}{2}) ?
  -20*log10(sqrt(\n@pow{#1}{2} + \n@pow{#2}{2})) :
  -20*log10(t)
  )}
\newcommand*{\MagPoleAsymp}{\MagPoleLin}
\newcommand*{\PhPole}[2]{((#1 > 0 ? (#2 > 0 ?
  (\n@mod@p{-atan2((t - (#2)),-(#1))}{2*pi}) :
  (-atan2((t - (#2)),-(#1)))) :
  (-atan2((t - (#2)),-(#1))))*\ph@scale)}
\newcommand*{\PhPoleLin}[2]{
  ((abs(#1)+abs(#2) == 0 ? -pi/2 :
  (t < (sqrt(\n@pow{#1}{2} + \n@pow{#2}{2}) /
    (\n@pow{10}{sqrt(\n@pow{#1}{2}/(\n@pow{#1}{2} + \n@pow{#2}{2}))})) ?
  (-atan2(-(#2),-(#1))) :
  (t >= (sqrt(\n@pow{#1}{2} + \n@pow{#2}{2}) *
    (\n@pow{10}{sqrt(\n@pow{#1}{2}/(\n@pow{#1}{2} + \n@pow{#2}{2}))})) ?
  (#2>0?(#1>0?3*pi/2:-pi/2):-pi/2) :
  (-atan2(-(#2),-(#1)) + (log10(t/(sqrt(\n@pow{#1}{2} + \n@pow{#2}{2}) /
    (\n@pow{10}{sqrt(\n@pow{#1}{2}/(\n@pow{#1}{2} +
    \n@pow{#2}{2}))}))))*((#2>0?(#1>0?3*pi/2:-pi/2):-pi/2) + atan2(-(#2),-(#1)))/
    (log10(\n@pow{10}{sqrt((4*\n@pow{#1}{2})/
    (\n@pow{#1}{2} + \n@pow{#2}{2}))}))))))*\ph@scale)}
\newcommand*{\PhPoleAsymp}[2]{((t < (sqrt(\n@pow{#1}{2} + \n@pow{#2}{2})) ?
  (-atan2(-(#2),-(#1))) :
  (#2>0?(#1>0?3*pi/2:-pi/2):-pi/2))*\ph@scale)}
\newcommand*{\MagZero}{0-\MagPole}
\newcommand*{\MagZeroLin}{0-\MagPoleLin}
\newcommand*{\MagZeroAsymp}{0-\MagPoleAsymp}
\newcommand*{\PhZero}{0-\PhPole}
\newcommand*{\PhZeroLin}{0-\PhPoleLin}
\newcommand*{\PhZeroAsymp}{0-\PhPoleAsymp}
\newcommand*{\MagCSPoles}[2]{(-20*log10(sqrt(\n@pow{\n@pow{#2}{2}
    - \n@pow{t}{2}}{2} + \n@pow{2*#1*#2*t}{2})))}
\newcommand*{\MagCSPolesLin}[2]{(t < #2 ? -40*log10(#2) : - 40*log10(t))}
\newcommand*{\MagCSPolesAsymp}{\MagCSPolesLin}
\newcommand*{\PhCSPoles}[2]{((-atan2((2*(#1)*(#2)*t),(\n@pow{#2}{2}
  - \n@pow{t}{2})))*\ph@scale)}
\newcommand*{\PhCSPolesLin}[2]{((t < (#2 / (\n@pow{10}{abs(#1)})) ?
  0 :
  (t >= (#2 * (\n@pow{10}{abs(#1)})) ?
  (#1>0 ? -pi : pi) :
  (#1>0 ? (-pi*(log10(t*(\n@pow{10}{#1})/#2))/(2*#1)) :
    (pi*(log10(t*(\n@pow{10}{abs(#1)})/#2))/(2*abs(#1))))))*\ph@scale)}
\newcommand*{\PhCSPolesAsymp}[2]{((#1>0?(t<#2?0:-pi):(t<#2?0:pi))*\ph@scale)}
\newcommand*{\MagCSZeros}{0-\MagCSPoles}
\newcommand*{\MagCSZerosLin}{0-\MagCSPolesLin}
\newcommand*{\MagCSZerosAsymp}{0-\MagCSPolesAsymp}
\newcommand*{\PhCSZeros}{0-\PhCSPoles}
\newcommand*{\PhCSZerosLin}{0-\PhCSPolesLin}
\newcommand*{\PhCSZerosAsymp}{0-\PhCSPolesAsymp}
\newcommand*{\MagCSPolesPeak}[3][]{
  \draw[#1,->] (axis cs:{#3},{-40*log10(#3)}) --
  (axis cs:{#3},{-40*log10(#3)-20*log10(2*abs(#2))})
}
\newcommand*{\MagCSZerosPeak}[3][]{
  \draw[#1,->] (axis cs:{#3},{40*log10(#3)}) --
  (axis cs:{#3},{40*log10(#3)+20*log10(2*abs(#2))})
}
\newcommand*{\MagSOPoles}[2]{
  (-20*log10(sqrt(\n@pow{#2 - \n@pow{t}{2}}{2} + \n@pow{#1*t}{2})))}
\newcommand*{\MagSOPolesLin}[2]{
  (t < sqrt(abs(#2)) ? -20*log10(abs(#2)) : - 40*log10(t))}
\newcommand*{\MagSOPolesAsymp}{\MagSOPolesLin}
\newcommand*{\PhSOPoles}[2]{((-atan2((#1)*t,((#2) - \n@pow{t}{2})))*\ph@scale)}
\newcommand*{\PhSOPolesLin}[2]{((#2>0 ?
  \PhCSPolesLin{(#1/(2*sqrt(#2)))}{(sqrt(#2))} :
  (#1>0 ? -pi : pi)))}
\newcommand*{\PhSOPolesAsymp}[2]{((#2>0 ?
  \PhCSPolesAsymp{(#1/(2*sqrt(#2)))}{(sqrt(#2))} :
  (#1>0 ? -pi : pi)))}
\newcommand*{\MagSOZeros}{0-\MagSOPoles}
\newcommand*{\MagSOZerosLin}{0-\MagSOPolesLin}
\newcommand*{\MagSOZerosAsymp}{0-\MagSOPolesAsymp}
\newcommand*{\PhSOZeros}{0-\PhSOPoles}
\newcommand*{\PhSOZerosLin}{0-\PhSOPolesLin}
\newcommand*{\PhSOZerosAsymp}{0-\PhSOPolesAsymp}
\newcommand*{\MagSOPolesPeak}[3][]{
  \draw[#1,->] (axis cs:{sqrt(abs(#3))},{-20*log10(abs(#3))}) --
  (axis cs:{sqrt(abs(#3))},{-20*log10(abs(#3)) -
    20*log10(abs(#2/sqrt(abs(#3))))});
}
\newcommand*{\MagSOZerosPeak}[3][]{
  \draw[#1,->] (axis cs:{sqrt(abs(#3))},{20*log10(abs(#3))}) --
  (axis cs:{sqrt(abs(#3))},{20*log10(abs(#3)) +
    20*log10(abs(#2/sqrt(abs(#3))))});
}
\newcommand{\BodeZPK}[4][approx/true]{
  \parse@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
  \build@ZPK@plot{\func@mag}{\func@ph}{\opt@approx}{#2}
    \edef\temp@cmd{\noexpand\begin{groupplot}[
      bode@style,
      xmin=#3,
      xmax=#4,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=2.5cm,
      xmode=log,
      group style = {group size = 1 by 2,vertical sep=0.25cm},
      \opt@group
    ]}
    \temp@cmd
      \edef\temp@mag@cmd{\noexpand\nextgroupplot [ylabel={Gain (dB)}, xmajorticks=false, \optmag@axes]
      \noexpand\addplot [freq@filter, variable=t, thick, \optmag@plot]}
      \edef\temp@ph@cmd{\noexpand\nextgroupplot [ph@y@label, freq@label, \optph@axes]
      \noexpand\addplot [freq@filter, variable=t, thick, trig format plots=rad, \optph@plot]}
      \if@pgfarg
        \temp@mag@cmd {\func@mag};
        \optmag@commands
        \temp@ph@cmd {\func@ph};
        \optph@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@mag@cmd gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set dummy t;
          set logscale x 10;
          set xrange [#3*\freq@scale:#4*\freq@scale];
          set samples \pgfkeysvalueof{/pgfplots/samples};
          plot \func@mag;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($1/(\freq@scale)):($2);
        };
        \optmag@commands
        \stepcounter{gnuplot@id}
        \temp@ph@cmd gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set dummy t;
          set logscale x 10;
          set xrange [#3*\freq@scale:#4*\freq@scale];
          set samples \pgfkeysvalueof{/pgfplots/samples};
          plot \func@ph;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($1/(\freq@scale)):($2);
        };
        \optph@commands
      \fi
    \end{groupplot}
  \end{tikzpicture}
}
\AtBeginDocument{%
  \if@babel
  \let\Orig@BodeZPK\BodeZPK
  \renewcommand{\BodeZPK}{%
    \expandafter\shorthandoff\expandafter{\shorthand@list}
    \BodeZPK@Shorthandoff
  }
  \newcommand{\BodeZPK@Shorthandoff}[4][]{%
    \Orig@BodeZPK[#1]{#2}{#3}{#4}
    \expandafter\shorthandon\expandafter{\shorthand@list}
  }
  \fi
}
\newcommand{\BodeTF}[4][]{
  \parse@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
  \build@TF@plot{\func@mag}{\func@ph}{#2}
    \edef\temp@cmd{\noexpand\begin{groupplot}[
      bode@style,
      xmin=#3,
      xmax=#4,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=2.5cm,
      xmode=log,
      group style = {group size = 1 by 2,vertical sep=0.25cm},
      \opt@group
    ]}
    \temp@cmd
      \edef\temp@mag@cmd{\noexpand\nextgroupplot [ylabel={Gain (dB)}, xmajorticks=false, \optmag@axes]
      \noexpand\addplot [freq@filter, variable=t, thick, \optmag@plot]}
      \edef\temp@ph@cmd{\noexpand\nextgroupplot [ph@y@label, freq@label, \optph@axes]
      \noexpand\addplot [freq@filter, variable=t, thick, trig format plots=rad, \optph@plot]}
      \if@pgfarg
        \temp@mag@cmd {\func@mag};
        \optmag@commands
        \temp@ph@cmd {\n@mod{\func@ph}{2*pi*\ph@scale}};
        \optph@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@mag@cmd gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set dummy t;
          set logscale x 10;
          set xrange [#3*\freq@scale:#4*\freq@scale];
          set samples \pgfkeysvalueof{/pgfplots/samples};
          plot \func@mag;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($1/(\freq@scale)):($2);
        };
        \optmag@commands
        \stepcounter{gnuplot@id}
        \temp@ph@cmd gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set dummy t;
          set logscale x 10;
          set trange [#3*\freq@scale:#4*\freq@scale];
          set samples \pgfkeysvalueof{/pgfplots/samples};
          plot '+' using (t) : ((\func@ph)/(\ph@scale)) smooth unwrap;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($1/(\freq@scale)):($2*\ph@scale);
        };
        \optph@commands
      \fi
    \end{groupplot}
  \end{tikzpicture}
}
\AtBeginDocument{
  \if@babel
  \let\Orig@BodeTF\BodeTF
  \renewcommand{\BodeTF}{%
    \expandafter\shorthandoff\expandafter{\shorthand@list}
    \BodeTF@Shorthandoff
  }
  \newcommand{\BodeTF@Shorthandoff}[4][]{%
    \Orig@BodeTF[#1]{#2}{#3}{#4}
    \expandafter\shorthandon\expandafter{\shorthand@list}
  }
  \fi
}
\newcommand{\addBodeZPKPlots}[3][true/{}]{
  \foreach \approx/\opt in {#1} {
    \gdef\plot@macro{}
    \gdef\temp@macro{}
    \ifnum\pdf@strcmp{#2}{phase}=0
      \build@ZPK@plot{\temp@macro}{\plot@macro}{\approx}{#3}
    \else
      \build@ZPK@plot{\plot@macro}{\temp@macro}{\approx}{#3}
    \fi
    \if@pgfarg
      \edef\temp@cmd{\noexpand\addplot [freq@filter, domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, thick, trig format plots=rad, \opt]}
      \temp@cmd {\plot@macro};
    \else
      \stepcounter{gnuplot@id}
      \edef\temp@cmd{\noexpand\addplot [variable=t, thick, \opt]}
      \temp@cmd gnuplot [raw gnuplot, gnuplot@prefix]
      { set table $meta;
        set dummy t;
        set logscale x 10;
        set xrange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
        set samples \pgfkeysvalueof{/pgfplots/samples};
        plot \plot@macro;
        set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
        plot "$meta" using ($1/(\freq@scale)):($2);
      };
    \fi
  }
}
\newcommand{\addBodeTFPlot}[3][thick]{
  \gdef\plot@macro{}
  \gdef\temp@macro{}
  \ifnum\pdf@strcmp{#2}{phase}=0
    \build@TF@plot{\temp@macro}{\plot@macro}{#3}
  \else
    \build@TF@plot{\plot@macro}{\temp@macro}{#3}
  \fi
  \if@pgfarg
    \ifnum\pdf@strcmp{#2}{phase}=0
      \edef\temp@cmd{\noexpand\addplot [freq@filter, domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
      \temp@cmd {\n@mod{\plot@macro}{2*pi}};
    \else
      \edef\temp@cmd{\noexpand\addplot [freq@filter, domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, #1]}
      \temp@cmd {\plot@macro};
    \fi
  \else
    \stepcounter{gnuplot@id}
    \ifnum\pdf@strcmp{#2}{phase}=0
      \addplot [variable=t, #1] gnuplot [raw gnuplot, gnuplot@prefix]
      { set table $meta;
        set dummy t;
        set logscale x 10;
        set trange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
        set samples \pgfkeysvalueof{/pgfplots/samples};
        plot '+' using (t) : ((\plot@macro)/(\ph@scale)) smooth unwrap;
        set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
        plot "$meta" using ($1/(\freq@scale)):($2*\ph@scale);
      };
    \else
      \addplot [variable=t, #1] gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set dummy t;
          set logscale x 10;
          set xrange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
          set samples \pgfkeysvalueof{/pgfplots/samples};
          plot \plot@macro;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($1/(\freq@scale)):($2);
        };
    \fi
  \fi
}
\newcommand{\addBodeComponentPlot}[2][thick]{
  \if@pgfarg
    \edef\temp@cmd{\noexpand\addplot [freq@filter, domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
    \temp@cmd {#2};
  \else
    \stepcounter{gnuplot@id}
    \addplot [variable=t, #1] gnuplot [raw gnuplot, gnuplot@prefix]
    { set table $meta;
      set dummy t;
      set logscale x 10;
      set xrange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
      set samples \pgfkeysvalueof{/pgfplots/samples};
      plot #2;
      set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
      plot "$meta" using ($1/(\freq@scale)):($2);
    };
  \fi
}
\AtBeginDocument{%
  \if@babel
    \AddToHook{env/BodePhPlot/begin}{\expandafter\shorthandoff\expandafter{\shorthand@list}}
    \AddToHook{env/BodePhPlot/end}{\expandafter\shorthandon\expandafter{\shorthand@list}}
  \fi
}
\NewDocumentEnvironment{BodePhPlot}{O{}mm+b}{
  \parse@env@opt{#1}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \edef\temp@cmd{\noexpand\begin{semilogxaxis}[
      ph@y@label,
      freq@label,
      bode@style,
      xmin={#2},
      xmax={#3},
      domain=#2:#3,
      height=2.5cm,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      #4
    \end{semilogxaxis}
  \end{tikzpicture}
}{}
\AtBeginDocument{%
  \if@babel
    \AddToHook{env/BodeMagPlot/begin}{\expandafter\shorthandoff\expandafter{\shorthand@list}}
    \AddToHook{env/BodeMagPlot/end}{\expandafter\shorthandon\expandafter{\shorthand@list}}
  \fi
}
\NewDocumentEnvironment{BodeMagPlot}{O{}mm+b}{
  \parse@env@opt{#1}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \edef\temp@cmd{\noexpand\begin{semilogxaxis}[
      bode@style,
      freq@label,
      xmin={#2},
      xmax={#3},
      domain=#2:#3,
      height=2.5cm,
      ylabel={Gain (dB)},
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      #4
    \end{semilogxaxis}
  \end{tikzpicture}
}{}
\AtBeginDocument{%
  \if@babel
    \AddToHook{env/BodePlot/begin}{\expandafter\shorthandoff\expandafter{\shorthand@list}}
    \AddToHook{env/BodePlot/end}{\expandafter\shorthandon\expandafter{\shorthand@list}}
  \fi
}
\NewDocumentEnvironment{BodePlot}{O{}mm+b}{
  \parse@env@opt{#1}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \edef\temp@cmd{\noexpand\begin{semilogxaxis}[
      bode@style,
      freq@label,
      xmin={#2},
      xmax={#3},
      domain=#2:#3,
      height=2.5cm,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      #4
    \end{semilogxaxis}
  \end{tikzpicture}
}{}
\newcommand*{\add@feature}[3]{
  \ifcat$\detokenize\expandafter{#1}$
    \xdef#1{\unexpanded\expandafter{#1 0+#2}}
  \else
    \xdef#1{\unexpanded\expandafter{#1+#2}}
  \fi
  \foreach \y [count=\n] in #3 {
    \xdef#1{\unexpanded\expandafter{#1}{\y}}
    \xdef\Last@LoopValue{\n}
  }
  \ifnum\Last@LoopValue=1
    \xdef#1{\unexpanded\expandafter{#1}{0}}
  \fi
}
\newcommand{\build@ZPK@plot}[4]{
  \foreach \feature/\values in {#4} {
    \ifnum\pdf@strcmp{\feature}{z}=0
      \foreach \z in \values {
        \ifnum\pdf@strcmp{#3}{linear}=0
          \add@feature{#2}{\PhZeroLin}{\z}
          \add@feature{#1}{\MagZeroLin}{\z}
        \else
          \ifnum\pdf@strcmp{#3}{asymptotic}=0
            \add@feature{#2}{\PhZeroAsymp}{\z}
            \add@feature{#1}{\MagZeroAsymp}{\z}
          \else
            \add@feature{#2}{\PhZero}{\z}
            \add@feature{#1}{\MagZero}{\z}
          \fi
        \fi
      }
    \fi
    \ifnum\pdf@strcmp{\feature}{p}=0
      \foreach \p in \values {
        \ifnum\pdf@strcmp{#3}{linear}=0
          \add@feature{#2}{\PhPoleLin}{\p}
          \add@feature{#1}{\MagPoleLin}{\p}
        \else
          \ifnum\pdf@strcmp{#3}{asymptotic}=0
            \add@feature{#2}{\PhPoleAsymp}{\p}
            \add@feature{#1}{\MagPoleAsymp}{\p}
          \else
            \add@feature{#2}{\PhPole}{\p}
            \add@feature{#1}{\MagPole}{\p}
          \fi
        \fi
      }
    \fi
    \ifnum\pdf@strcmp{\feature}{k}=0
      \ifnum\pdf@strcmp{#3}{linear}=0
        \add@feature{#2}{\PhKLin}{\values}
        \add@feature{#1}{\MagKLin}{\values}
      \else
        \ifnum\pdf@strcmp{#3}{asymptotic}=0
          \add@feature{#2}{\PhKAsymp}{\values}
          \add@feature{#1}{\MagKAsymp}{\values}
        \else
          \add@feature{#2}{\PhK}{\values}
          \add@feature{#1}{\MagK}{\values}
        \fi
      \fi
    \fi
    \ifnum\pdf@strcmp{\feature}{d}=0
      \ifnum\pdf@strcmp{#3}{linear}=0
        \PackageError {bodeplot} {Linear approximation for pure delays is not
        supported.} {Plot the true Bode plot using `true' instead of `linear'.}
      \else
        \ifnum\pdf@strcmp{#3}{asymptotic}=0
          \PackageError {bodeplot} {Asymptotic approximation for pure delays is not
          supported.} {Plot the true Bode plot using `true' instead of `asymptotic'.}
        \else
          \ifdim\values pt < 0pt
            \PackageError {bodeplot} {Delay needs to be a positive number.}
          \fi
          \add@feature{#2}{\PhDel}{\values}
          \add@feature{#1}{\MagDel}{\values}
        \fi
      \fi
    \fi
  }
}
\newcommand{\build@TF@plot}[3]{
  \gdef\num@real{0}
  \gdef\num@im{0}
  \gdef\den@real{0}
  \gdef\den@im{0}
  \gdef\loop@delay{0}
  \foreach \feature/\values in {#3} {
    \ifnum\pdf@strcmp{\feature}{num}=0
      \foreach \numcoeff [count=\numpow] in \values {
        \xdef\num@degree{\numpow}
      }
      \foreach \numcoeff [count=\numpow] in \values {
        \pgfmathtruncatemacro{\currentdegree}{\num@degree-\numpow}
        \ifnum\currentdegree = 0
          \xdef\num@real{\num@real+\numcoeff}
        \else
          \ifodd\currentdegree
            \xdef\num@im{\num@im+(\numcoeff*(\n@pow{-1}{(\currentdegree-1)/2})*%
              (\n@pow{t}{\currentdegree}))}
          \else
            \xdef\num@real{\num@real+(\numcoeff*(\n@pow{-1}{(\currentdegree)/2})*%
              (\n@pow{t}{\currentdegree}))}
          \fi
        \fi
      }
    \fi
    \ifnum\pdf@strcmp{\feature}{den}=0
      \foreach \dencoeff [count=\denpow] in \values {
        \xdef\den@degree{\denpow}
      }
      \foreach \dencoeff [count=\denpow] in \values {
        \pgfmathtruncatemacro{\currentdegree}{\den@degree-\denpow}
        \ifnum\currentdegree = 0
          \xdef\den@real{\den@real+\dencoeff}
        \else
          \ifodd\currentdegree
            \xdef\den@im{\den@im+(\dencoeff*(\n@pow{-1}{(\currentdegree-1)/2})*%
              (\n@pow{t}{\currentdegree}))}
          \else
            \xdef\den@real{\den@real+(\dencoeff*(\n@pow{-1}{(\currentdegree)/2})*%
              (\n@pow{t}{\currentdegree}))}
          \fi
        \fi
      }
    \fi
    \ifnum\pdf@strcmp{\feature}{d}=0
      \xdef\loop@delay{\values}
    \fi
  }
  \xdef#2{((atan2((\num@im),(\num@real))-atan2((\den@im),%
    (\den@real))-\loop@delay*t)*(\ph@scale))}
  \xdef#1{(20*log10(sqrt((\n@pow{\num@real}{2})+(\n@pow{\num@im}{2})))-%
    20*log10(sqrt((\n@pow{\den@real}{2})+(\n@pow{\den@im}{2}))))}
}
\newcommand{\parse@opt}[1]{
  \gdef\optmag@axes{}
  \gdef\optph@axes{}
  \gdef\optph@plot{}
  \gdef\optmag@plot{}
  \gdef\opt@group{}
  \gdef\opt@approx{}
  \gdef\optph@commands{}
  \gdef\optmag@commands{}
  \gdef\opt@tikz{}
  \foreach \obj/\typ/\opt in {#1} {
    \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{plot}=0
      \ifnum\pdf@strcmp{\unexpanded\expandafter{\typ}}{mag}=0
        \xdef\optmag@plot{\unexpanded\expandafter{\opt}}
      \else
        \ifnum\pdf@strcmp{\unexpanded\expandafter{\typ}}{ph}=0
          \xdef\optph@plot{\unexpanded\expandafter{\opt}}
        \else
          \xdef\optmag@plot{\unexpanded\expandafter{\opt}}
          \xdef\optph@plot{\unexpanded\expandafter{\opt}}
        \fi
      \fi
    \else
      \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{axes}=0
        \ifnum\pdf@strcmp{\unexpanded\expandafter{\typ}}{mag}=0
          \xdef\optmag@axes{\unexpanded\expandafter{\opt}}
        \else
          \ifnum\pdf@strcmp{\unexpanded\expandafter{\typ}}{ph}=0
            \xdef\optph@axes{\unexpanded\expandafter{\opt}}
          \else
            \xdef\optmag@axes{\unexpanded\expandafter{\opt}}
            \xdef\optph@axes{\unexpanded\expandafter{\opt}}
          \fi
        \fi
      \else
        \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{group}=0
          \xdef\opt@group{\unexpanded\expandafter{\opt}}
        \else
          \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{approx}=0
            \xdef\opt@approx{\unexpanded\expandafter{\opt}}
          \else
            \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{commands}=0
              \ifnum\pdf@strcmp{\unexpanded\expandafter{\typ}}{ph}=0
                \xdef\optph@commands{\unexpanded\expandafter{\opt}}
              \else
                \xdef\optmag@commands{\unexpanded\expandafter{\opt}}
              \fi
            \else
              \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{tikz}=0
                \xdef\opt@tikz{\unexpanded\expandafter{\opt}}
              \else
                \xdef\optmag@plot{\unexpanded\expandafter{\optmag@plot},
                  \unexpanded\expandafter{\obj}}
                \xdef\optph@plot{\unexpanded\expandafter{\optph@plot},
                  \unexpanded\expandafter{\obj}}
              \fi
            \fi
          \fi
        \fi
      \fi
    \fi
  }
}
\newcommand{\parse@env@opt}[1]{
  \gdef\opt@axes{}
  \gdef\opt@tikz{}
  \foreach \obj/\opt in {#1} {
    \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{axes}=0
      \xdef\opt@axes{\unexpanded\expandafter{\opt}}
    \else
      \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{tikz}=0
        \xdef\opt@tikz{\unexpanded\expandafter{\opt}}
      \else
        \xdef\opt@axes{\unexpanded\expandafter{\opt@axes},
          \unexpanded\expandafter{\obj}}
      \fi
    \fi
  }
}
\newcommand{\NyquistZPK}[4][]{
  \parse@N@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
  \build@ZPK@plot{\func@mag}{\func@ph}{}{#2}
    \edef\temp@cmd{\noexpand\begin{axis}[
      bode@style,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=5cm,
      xlabel={$\Re$},
      ylabel={$\Im$},
      samples=500,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      \addplot [only marks,mark=+,thick,red] (-1 , 0);
      \edef\temp@cmd{\noexpand\addplot [variable=t, thick, trig format plots=rad, \unexpanded\expandafter{\opt@plot}]}
      \if@pgfarg
        \temp@cmd ( {\n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale))},
          {\n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))} );
        \opt@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@cmd gnuplot [parametric, gnuplot@prefix] {
          \n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale)),
          \n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))
        };
        \opt@commands
      \fi
    \end{axis}
  \end{tikzpicture}
}
\AtBeginDocument{%
  \if@babel
  \let\Orig@NyquistZPK\NyquistZPK
  \renewcommand{\NyquistZPK}{%
    \expandafter\shorthandoff\expandafter{\shorthand@list}
    \NyquistZPK@Shorthandoff
  }
  \newcommand{\NyquistZPK@Shorthandoff}[4][]{%
    \Orig@NyquistZPK[#1]{#2}{#3}{#4}
    \expandafter\shorthandon\expandafter{\shorthand@list}
  }
  \fi
}
\newcommand{\NyquistTF}[4][]{
  \parse@N@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \build@TF@plot{\func@mag}{\func@ph}{#2}
    \edef\temp@cmd{\noexpand\begin{axis}[
      bode@style,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=5cm,
      xlabel={$\Re$},
      ylabel={$\Im$},
      samples=500,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      \addplot [only marks, mark=+, thick, red] (-1 , 0);
      \edef\temp@cmd{\noexpand\addplot [variable=t, thick, trig format plots=rad, \unexpanded\expandafter{\opt@plot}]}
      \if@pgfarg
        \temp@cmd ( {\n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale))},
          {\n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))} );
        \opt@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@cmd gnuplot [parametric, gnuplot@prefix] {
          \n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale)),
          \n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))
        };
        \opt@commands
      \fi
    \end{axis}
  \end{tikzpicture}
}
\AtBeginDocument{%
  \if@babel
  \let\Orig@NyquistTF\NyquistTF
  \renewcommand{\NyquistTF}{%
    \expandafter\shorthandoff\expandafter{\shorthand@list}
    \NyquistTF@Shorthandoff
  }
  \newcommand{\NyquistTF@Shorthandoff}[4][]{%
    \Orig@NyquistTF[#1]{#2}{#3}{#4}
    \expandafter\shorthandon\expandafter{\shorthand@list}
  }
  \fi
}
\newcommand{\addNyquistZPKPlot}[2][]{
  \gdef\func@mag{}
  \gdef\func@ph{}
  \build@ZPK@plot{\func@mag}{\func@ph}{}{#2}
  \if@pgfarg
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
    \temp@cmd ( {\n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale))},
      {\n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))} );
  \else
    \stepcounter{gnuplot@id}
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, #1]}
    \temp@cmd gnuplot [parametric, gnuplot@prefix] {
      \n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale)),
      \n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))
    };
  \fi
}
\newcommand{\addNyquistTFPlot}[2][]{
  \gdef\func@mag{}
  \gdef\func@ph{}
  \build@TF@plot{\func@mag}{\func@ph}{#2}
  \if@pgfarg
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
    \temp@cmd ( {\n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale))},
      {\n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))} );
  \else
    \stepcounter{gnuplot@id}
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, #1]}
    \temp@cmd gnuplot [parametric, gnuplot@prefix]{
      \n@pow{10}{((\func@mag)/20)}*cos((\func@ph)/(\ph@scale)),
      \n@pow{10}{((\func@mag)/20)}*sin((\func@ph)/(\ph@scale))
    };
  \fi
}
\AtBeginDocument{%
  \if@babel
    \AddToHook{env/NyquistPlot/begin}{\expandafter\shorthandoff\expandafter{\shorthand@list}}
    \AddToHook{env/NyquistPlot/end}{\expandafter\shorthandon\expandafter{\shorthand@list}}
  \fi
}
\NewDocumentEnvironment{NyquistPlot}{O{}mm+b}{
  \parse@env@opt{#1}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \edef\temp@cmd{\noexpand\begin{axis}[
      bode@style,
      height=5cm,
      domain=#2:#3,
      xlabel={$\Re$},
      ylabel={$\Im$},
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      \addplot [only marks,mark=+,thick,red] (-1 , 0);
      #4
    \end{axis}
  \end{tikzpicture}
}{}
\newcommand{\parse@N@opt}[1]{
  \gdef\opt@axes{}
  \gdef\opt@plot{}
  \gdef\opt@commands{}
  \gdef\opt@tikz{}
  \foreach \obj/\opt in {#1} {
    \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{axes}=0
      \xdef\opt@axes{\unexpanded\expandafter{\opt}}
    \else
      \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{plot}=0
        \xdef\opt@plot{\unexpanded\expandafter{\opt}}
      \else
        \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{commands}=0
          \xdef\opt@commands{\unexpanded\expandafter{\opt}}
        \else
          \ifnum\pdf@strcmp{\unexpanded\expandafter{\obj}}{tikz}=0
            \xdef\opt@tikz{\unexpanded\expandafter{\opt}}
          \else
            \xdef\opt@plot{\unexpanded\expandafter{\opt@plot},
              \unexpanded\expandafter{\obj}}
          \fi
        \fi
      \fi
    \fi
  }
}
\newcommand{\NicholsZPK}[4][]{
  \parse@N@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
  \build@ZPK@plot{\func@mag}{\func@ph}{}{#2}
    \edef\temp@cmd{\noexpand\begin{axis}[
      ph@x@label,
      bode@style,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=5cm,
      ylabel={Gain (dB)},
      samples=500,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      \edef\temp@cmd{\noexpand\addplot [variable=t, thick, trig format plots=rad, \opt@plot]}
      \if@pgfarg
        \temp@cmd ( {\func@ph} , {\func@mag} );
        \opt@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@cmd gnuplot [raw gnuplot, gnuplot@prefix]
        { set table $meta;
          set logscale x 10;
          set dummy t;
          set samples \pgfkeysvalueof{/pgfplots/samples};
          set trange [#3*\freq@scale:#4*\freq@scale];
          plot '+' using (\func@mag) : ((\func@ph)/(\ph@scale));
          unset logscale x;
          set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
          plot "$meta" using ($2*\ph@scale):($1);
        };
        \opt@commands
      \fi
    \end{axis}
  \end{tikzpicture}
}
\AtBeginDocument{%
  \if@babel
  \let\Orig@NicholsZPK\NicholsZPK
  \renewcommand{\NicholsZPK}{%
    \expandafter\shorthandoff\expandafter{\shorthand@list}
    \NicholsZPK@Shorthandoff
  }
  \newcommand{\NicholsZPK@Shorthandoff}[4][]{%
    \Orig@NicholsZPK[#1]{#2}{#3}{#4}
    \expandafter\shorthandon\expandafter{\shorthand@list}
  }
  \fi
}
\newcommand{\NicholsTF}[4][]{
  \parse@N@opt{#1}
  \gdef\func@mag{}
  \gdef\func@ph{}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \build@TF@plot{\func@mag}{\func@ph}{#2}
    \edef\temp@cmd{\noexpand\begin{axis}[
      ph@x@label,
      bode@style,
      domain=#3*\freq@scale:#4*\freq@scale,
      height=5cm,
      ylabel={Gain (dB)},
      samples=500,
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      \edef\temp@cmd{\noexpand\addplot [variable=t, thick, trig format plots=rad, \opt@plot]}
      \if@pgfarg
        \temp@cmd ( {\n@mod{\func@ph}{2*pi*\ph@scale}} , {\func@mag} );
        \opt@commands
      \else
        \stepcounter{gnuplot@id}
        \temp@cmd gnuplot [raw gnuplot, gnuplot@prefix]
          { set table $meta1;
            set logscale x 10;
            set dummy t;
            set samples \pgfkeysvalueof{/pgfplots/samples};
            set trange [#3*\freq@scale:#4*\freq@scale];
            plot '+' using (\func@mag) : ((\func@ph)/(\ph@scale));
            unset logscale x;
            set table $meta2;
            plot "$meta1" using ($1):($2) smooth unwrap;
            set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
            plot "$meta2" using ($2*\ph@scale):($1);
          };
        \opt@commands
      \fi
    \end{axis}
  \end{tikzpicture}
}
\AtBeginDocument{
  \if@babel
    \let\Orig@NicholsTF\NicholsTF
    \renewcommand{\NicholsTF}{%
      \expandafter\shorthandoff\expandafter{\shorthand@list}
      \NicholsTF@Shorthandoff
    }
    \newcommand{\NicholsTF@Shorthandoff}[4][]{%
      \Orig@NicholsTF[#1]{#2}{#3}{#4}
      \expandafter\shorthandon\expandafter{\shorthand@list}
    }
    \AddToHook{env/NicholsChart/begin}{\expandafter\shorthandoff\expandafter{\shorthand@list}}
    \AddToHook{env/NicholsChart/end}{\expandafter\shorthandon\expandafter{\shorthand@list}}
  \fi
}
\NewDocumentEnvironment{NicholsChart}{O{}mm+b}{
  \parse@env@opt{#1}
  \edef\temp@cmd{\noexpand\begin{tikzpicture} [\unexpanded\expandafter{\opt@tikz}]}
  \temp@cmd
    \edef\temp@cmd{\noexpand\begin{axis}[
      ph@x@label,
      bode@style,
      domain=#2:#3,
      height=5cm,
      ylabel={Gain (dB)},
      \unexpanded\expandafter{\opt@axes}
    ]}
    \temp@cmd
      #4
    \end{axis}
  \end{tikzpicture}
}{}
\newcommand{\addNicholsZPKChart}[2][]{
  \gdef\func@mag{}
  \gdef\func@ph{}
  \build@ZPK@plot{\func@mag}{\func@ph}{}{#2}
  \if@pgfarg
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
    \temp@cmd ( {\func@ph} , {\func@mag} );
  \else
    \stepcounter{gnuplot@id}
    \addplot [#1] gnuplot [raw gnuplot, gnuplot@prefix]
    { set table $meta;
      set logscale x 10;
      set dummy t;
      set samples \pgfkeysvalueof{/pgfplots/samples};
      set trange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
      plot '+' using (\func@mag) : ((\func@ph)/(\ph@scale));
      unset logscale x;
      set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
      plot "$meta" using ($2*\ph@scale):($1);
    };
  \fi
}
\newcommand{\addNicholsTFChart}[2][]{
  \gdef\func@mag{}
  \gdef\func@ph{}
  \build@TF@plot{\func@mag}{\func@ph}{#2}
  \if@pgfarg
    \edef\temp@cmd{\noexpand\addplot [domain=\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale, variable=t, trig format plots=rad, #1]}
    \temp@cmd ( {\n@mod{\func@ph}{2*pi*\ph@scale}} , {\func@mag} );
  \else
    \stepcounter{gnuplot@id}
    \addplot [#1] gnuplot [raw gnuplot, gnuplot@prefix]
    { set table $meta1;
      set logscale x 10;
      set dummy t;
      set samples \pgfkeysvalueof{/pgfplots/samples};
      set trange [\freq@scale*\pgfkeysvalueof{/pgfplots/domain}*\freq@scale];
      plot '+' using (\func@mag) : ((\func@ph)/(\ph@scale));
      unset logscale x;
      set table $meta2;
      plot "$meta1" using ($1):($2) smooth unwrap;
      set table "\bodeplot@prefix\arabic{gnuplot@id}.table";
      plot "$meta2" using ($2*\ph@scale):($1);
    };
  \fi
}
\endinput
%%
%% End of file `bodeplot.sty'.