 LIB "control.lib";
// a WindTunnel example
ring A = (0,a, omega, zeta, k),(D1, delta),dp;
module R;
R = [D1+a, k*a*delta, 0, 0],
[0, D1, 1, 0],
[0, omega^2, D1+2*zeta*omega, omega^2];
R=transpose(R);
view(R);
==> D1+(a),(a*k)*delta,0 ,0 ,
==> 0 ,D1 ,1 ,0 ,
==> 0 ,(omega^2) ,D1+(2*omega*zeta),(omega^2)
view(control(R));
==> number of first nonzero Ext:
==>
==> 2
==>
==> controllable, not reflexive, image representation:
==>
==> (a*omega^2*k)*delta ,
==> (omega^2)*D1+(a*omega^2) ,
==> (omega^2)*D1^2+(a*omega^2)*D1 ,
==> D1^3+(a+2*omega*zeta)*D1^2+(2*a*omega*zeta+omega^2)*D1+(a*omega^2)
==>
==> dimension of the system:
==>
==> 2
==>
==> Parameter constellations which might lead to a noncontrollable system:
==>
==> a,k,omega
==>
==>
