function dx = PoincareBendixson(t,x)

dx = zeros(2,1);
%state equation
b = -2; % no limit cycle
b = 1; % limit cycle
dx(1) = x(2);
dx(2) = x(1)+b*x(2)-x(1)^2*x(2)-x(1)^3;