function test_sos_unc


sdpvar x y a t 
gamma = sdpvar(1)

p = a*x^4 + y^4+x*y + 1+gamma
sol = solvesos(set(uncertain(gamma))+set(-1/2 < gamma < 1/2)+set(sos(p - t))+set(ismember(a,[3 4 5]))+set(4.6 > abs(a) > 3.5),-t,[],[a t])

mbg_asserttolequal(sol.problem,0);
mbg_asserttolequal(double(t),0.4375,1e-5);