function xdot = compesti(t,x)
%This file is called by the script ESTI_EX
% to produce the unit-step response of Example 3.12

A=[-25 0 -0.5;0 0 1;46.296 0 0];
c=[0 0.02 0];
b=[5;0;0];
k=[-2.6 0.0864 -0.005];
u=-k*x+1;
xdot=A*x+b*u;