Obtaining the solution of a polynomial in SOS

[Qu Cao]

This might be a simple question, but it really annoys me...

I solved a SOS problem that involves a polynomial, i.e.,

[s1,c1] = polynomial([q1 q2 q3 q4 q5 q6],2);

...

prog1 = [sos(J-gmama), sos(s1)];

solvesos(prog1,[],ops,[c1])

So how can I get access to the solution of s1?

I tried to use

sdisplay(s1),

but it gives

c1(1)+q1*c1(2)+q2*c1(3)+q3*c1(4)+q4*c1(5)+q5*c1(6)+q6*c1(7)+q1^2*c1(8)+q1*q2*c1(9)+q2^2*c1(10)+q1*q3*c1(11)+q2*q3*c1(12)+q3^2*c1(13)+q1*q4*c1(14)+q2*q4*c1(15)+q3*q4*c1(16)+q4^2*c1(17)+q1*q5*c1(18)+q2*q5*c1(19)+q3*q5*c1(20)+q4*q5*c1(21)+q5^2*c1(22)+q1*q6*c1(23)+q2*q6*c1(24)+q3*q6*c1(25)+q4*q6*c1(26)+q5*q6*c1(27)+q6^2*c1(28)

which is still the initial value.

[Johan Löfberg]

replace(s1,c1,value(c1))

or

[s1,c1,v1] = polynomial([q1 q2 q3 q4 q5 q6],2);
...
value(c1)'*v1