function xp = ex4(t,x,flag,p)

xp=[-x(1)+p*(1-x(1))*exp(x(2));
   -3*x(2)+14*p*(1-x(1))*exp(x(2))];
