function dxdt=windup(t,x)
%Differential equation for simulation of windup
r=1;a=0.2;
u=r-x(1)+x(2);v=sat(u,a);e=r-x(1);
dxdt=[v;e];