// Updated(3-8-07)
// 6.5

xi = 0.1*rand(1,10000,'normal');  
a = 1; b = [];
d = [1 1 -0.5];
ar = armac(a,b,d,1,1,1);
v = arsimul(ar,xi);
z = [v xi];

// Plot noise, plant output and ACF
subplot(2,1,1), plot(v(1:500))
label('',4,'','v',4)
subplot(2,1,2), plot(xi(1:500))
label('',4,'n','xi',4)
getf plotacf.sci;
xset('window',1)
plotacf(v,1,11,1);