function dxdt=relayfb(t,x)
%Differential equation for P4 with relay feedback
dxdt=[-x(1)+x(2);-x(2)+x(3);-x(3)+x(4);-x(4)-sign(x(1))];