function y = hopf(x,p)
%
y=x;
y(1) = p*x(1)+ x(2)-x(1)*x(1)*x(1);
y(2)=-x(1);