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

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

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

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

# data[].attributes.body
--> g3= syslin('c',(5*%s^2+3*%s+6)/(%s^3+11*%s^2+14*%s+12) )
--> num= numer (g3)
--> den= denom(g3)
--> zeros_g3= roots(num)
--> polos_g3= roots(den)
--> plzr(g3)

# 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
--> 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);
--> t=0:0.1:100;
--> u=ones(t);
--> y= dsimul(s1,u);
--> plot(t,y)

# data[].attributes.body
G(s) = C*(sI−A)^1*B

--> A= [0 1;-3 -4]
--> B= [0; 1]
--> C= [10 0]
--> s=poly(0,'s')
--> eye(2,2)*s
--> M=(eye(2,2)*s)-A
--> inv(M)
-->G= C*inv(M)*B