%Figure 7.8
omega=logspace(-2,2,801);
theta=1;
delay_g=delay(theta,omega);
err = msub(1,delay_g);
max=[2 3.14;2 100; 2 Inf];
theta=[1 1; .02 1; 2 Inf];
clf
subplot(221);
vplot('liv,lm',err,max,'--',1,':',theta,':')
text(0.3,0.015,'1/DELAY')
ylabel('Magnitude'), xlabel('Frequency')
omega=logspace(-2,2,81);
tau=1;
g=nd2sys(1,[tau 1]);
g_g=frsp(g,omega);
err2 = msub(1,g_g);
tau=[1 1; .02 1; 2 Inf];
subplot(222); 
vplot('liv,lm',err2,1,':',tau,':')
axis([.01,100,.01,10])
text(0.3,0.015,'1/TAU')
xlabel('Frequency')
