# data[].attributes.body
--> t=[0:0.1:30];
--> y=sin(3*t);
--> plot(t,y)

# data[].attributes.body
--> polinomio = poly ([7 2 3 0 1], 's', 'coeff');
--> roots(polinomio)

# data[].attributes.body
ans  =

0.7816583 + 1.65851i
0.7816583 - 1.65851i
-0.7816583 + 1.2129815i
-0.7816583 - 1.2129815i

# data[].attributes.body
--> p=poly ([-1 -5 -2 -8], 's', 'roots')

# data[].attributes.body
--> G1 = 1/(%s+1)
--> G2=(%s+2)/(%s+3)
--> MA=G1*G2;
--> H = 1
--> MF=MA/.H

# data[].attributes.body
--> G1 = 1/(%s+1)
--> G2=(%s+2)/(%s+3)
--> MA=G1*G2;
--> H = 1
--> MF=MA/.H
--> t=0.1:0.1:100;
--> u=ones(t);
--> y=csim(u,t,MF);
--> plot(t,y)

# data[].attributes.body
--> G3= syslin('c',(5*%s^2+3*%s+6)/(%s^3+11*%s^2+14*%s+12));
--> H=1;
--> MF=G3/.H;
--> s1=tf2ss(MF)

# data[].attributes.body
Continuação dos comandos da questão anterior

-> t=0:0.1:100;
--> u=ones(t);
--> y= dsimul(s1,u);
--> plot(t,y)

# data[].attributes.body
--> A= [0 1;-3 -4]
A  =

0.   1.
-3.  -4.
--> B= [0; 1]
B  =

0.
1.
--> C= [10 0]
C  =

10.   0.
--> s=poly(0,'s')
s  =

s
--> eye(2,2)*s
ans  =

s    0

0    s
--> M=(eye(2,2)*s)-A
M  =

s  -1

3   4 +s
--> inv(M)
ans  =

4 + s          1
-----------  -----------
2            2
3 + 4s + s   3 + 4s + s

-3           s
-----------  -----------
2            2
3 + 4s + s   3 + 4s + s
--> G= C*inv(M)*B