function dx=hopf(t,x)
u=0.5;
dx=[x(2)+x(1)*(u-x(1)*x(1)-x(2)*x(2));
    -x(1)+x(2)*(u-x(1)*x(1)-x(2)*x(2))];

