function [xprime] = test(t,x) global f11 f12; A = [0 -1 1 0]; B = [0 1]'; % Exosystem state w = exp(-t); % Asymptotic control u = -2*w + f11*x(1) + f12*x(2) - (f11+f12)*w; xprime=A*x+B*u;