num=[1 6.5 14 11.5 3];
den=[1 20 120 232 320];
u=mseries(100);
a=length(u);
[A,B]=getval(num,den);
[e,c,z,dif]=avlsm(u,A,B);
i=1:a;
figure(1)
plot(i,dif,'black')
figure(2)
for j=1:9
    subplot(3,3,j);
    plot(i,c(j,:),'red');
end