function xx1 = dynsu1(t,x)
global a b c d p
xx1(1,1) = p*(a*x(1)+b*x(2));xx1(2,1) =p*(c*x(1)+d*x(2));
