clc

close all

clear

syms s t

Fs = (2*s+1)/(s^3+7*s^2+10*s)

den = factor(s^3+7*s^2+10*s)

Ft = ilaplace(Fs, s, t);

Ft = vpa(Ft, 2)

ts = 0: 0.01: 10;

y = ones(1, 1001)*0.1;

Ft = subs(Ft, t, ts);

plot(ts, Ft, 'k', 'LineWidth', 2)

hold on;

plot(ts, y, 'r-.', 'LineWidth', 1)

legend('response', 'y=0.1');

grid on;