Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function [sys,x0]=model(t,x,u,flag)
- if (flag==0)
- x0=[1;-2;0;0];
- sys=[4;0;1;1;0;0];
- else if (flag==1)
- e=u-4*x(4);
- p=-1 .* (x(1) < -1) + (3 * x(1) + 2) .* (x(1) >= -1 & x(1) < 0) + 2 .* (x(1) >= 0);
- sys=[x(2)
- 2*e - 3*x(2)*x(1) - 2*x(1)
- x(3)
- 0.5*(0.2*p-x(3))];
- else if (flag==3)
- sys=x(4);
- else
- sys=[];
- end
Add Comment
Please, Sign In to add comment