function xdot=P4_1fun(t,x)
A=[-3 1;200 -100];
xdot=A*x;