solvesos(sos(const2),[],options,[vec(P)',K1,K2,K3])
solvesos(sos(-expr1(1:3,1:3)),[],options,[vec(P)',K1,K2,K3])
solvesos(sos(-expr1(1:4,1:4)),[],options,[vec(P)',K1,K2,K3])
solvesos(sos(-expr1(3:4,3:4)),[],options,[vec(P)',K1,K2,K3])