How to Plot Inverse Laplace with Heaviside?

Hello. I want to plot an inverse laplace transform using simple MATLAB functions however, I get an error message when I try to. Does this have to do with the heaviside functions? The equation is provided below:

clc
syms s t

Q = 10/(s*(s^2+100*s+1600))-10*exp(-3*s)/(s*(s^2+100*s+1600));
pretty(Q)

q = ilaplace(Q);
pretty(q)

plot(t,q)

Answer 1

[Star Strider]

Use the fplot function —

clc
syms s t

Q = 10/(s*(s^2+100*s+1600))-10*exp(-3*s)/(s*(s^2+100*s+1600));
pretty(Q)

          10                 exp(-3 s) 10
--------------------- - ---------------------
    2                       2
s (s  + 100 s + 1600)   s (s  + 100 s + 1600)

q = ilaplace(Q)

q = 

pretty(q)

exp(-80 t)   exp(-20 t)                    / exp(240 - 80 t)   exp(60 - 20 t)     1  \       1
---------- - ---------- - heaviside(t - 3) | --------------- - -------------- + ---- | 10 + ---
    480          120                       \       4800             1200        1600 /      160

figure
fplot(q, [0 5])
grid
xlabel('t')
ylabel('q(t)')
axis('padded')

.