x = sdpvar(2,1);
f = [-x(2);x(1)+(x(1)^2-1)*x(2)];
v = 1.5*x(1)^2-x(1)*x(2)+x(2)^2; gam = 1.0; ep = 1.0e-6;
s1 = polynomial(x,2,2);
co = -jacobian(v,x)*f-s1*(gam-v)-ep*x'*x;
Fs = sos(s1)+sos(co);
params = recover(setdiff(depends(Fs),depends(x)));
[sol,monos,Q] = solvesos(Fs,[],[],params)