an>
        <span class="comment">% Plot the step response</span>
        <span class="comment">%</span>

        T = linspace(0,2*delay,1000);
        <span class="keyword">for</span> inp = 1 : 6,
            figure(6*xi-5+inp);
            GQ = G + sparse(inp,inp,delay,n,n);
            A = -inv(C)*GQ;
            B = -A*ones(n,1);
            Y = simple_step(A,B,T(2),length(T));
            hold <span class="string">off</span>; plot(T,Y,<span class="string">'-'</span>);  hold <span class="string">on</spa