load pbam.txt;
x=pbam;

[z,p,k]=butter(10,.1);
[num,den]=zp2tf(z,p,k);
w=.1:.1:3;
%[m,p]=dbode(num,den,w);
%loglog(w,m);


%y=dlsim(num,den,x);
%plot(y(1:100));

xu=[];
for i=1:500
  xu=[xu; x(i);zeros(5,1)];
end

y=dlsim(num,den,xu);
plot(y(1:500));
