function xdot=php_pol(t,x)
xdot=[x(2);(1-x(1)^2)*x(2)-x(1)];