1340/CH2/EX2.3/2_3.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

clc;
s = %s;
syms t;
num = 32;
den = s*(s^2+12*s+32);
y = syslin('c',num/den);
disp(y);

[A] = pfss(y);G = 0;disp(A);
for k = 1:size(A)
    f(k) = ilaplace(A(k),s,t);
    G = G + f(k);
end

disp(G)