257/CH7/EX7.33/example_7_33.sce


1
2
3
4
5
6
7
8
9
10

syms s t
p=poly([8 1],'s','coeff');
q=poly([0 4 1],'s','coeff');
G=2*p/q //gain FACTOR=2
H=1
T=G/(1+G*H) 
R=1/s;
C=T*R
c=ilaplace(C,s,t)
disp(c,"c = ")