s=%s
L=polly([6 9 3],'s','c')
M=polly([32 56 28 4],'s','c')
a=roots(M)
b=roots(L)
//postac iloczynowa
//licznik
polfact(L)
//mianownik
polfact(M)

g=syslin('c',L/M)
plzr(g)
t=[0:0.001:10]
figure()
Y1=csim('impulse',t,g)
plot(t, Y1)
Y2=csim('step',t,g)
figure()
plot(t, Y2)
gz=g/(1+g)
e=1-gz
ge=syslin('c',e)
Ye=csim('impulse', t,ge)
plot(t,Ye)