function dy = thrusys(t,y)
% Edy=Ay+Bu ODE formulation
% dy=E^{-1}[Ay+Bu]
dy = zeros(30,1); u = 1;
global Am; global Bm;
dy = Am*y+Bm*u;