f=[8,2,1],
g=[1.0000,3.3000,9.9000,2.7000],
z=roots(f),
p=roots(g),
k=f(1)/g(1),
[r,p,h]=residue(f,g);
t=0:0.1:10;
x=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(2)*t);
plot(t,x);