clear;
 clc;
sdpvar b n
sdpvar pc0 pc1 pc2 pc3 pc4 pc5
sdpvar v1 v2 v3 v4 v5 v6
sdpvar u1 u2 u3 u4 u5
FeasibilityConstraints=[sos(pc0+pc1*(b)+pc2*(n)+pc3*(b)*(b)+pc4*(b)*(n)+pc5*(n)*(n)-(-2*b)+v1*(n-0)-u1*(0-b)+v2*(n-0)*b+v3*(n-0)*n),...
sos(n-(pc0+pc1*(b)+pc2*(n)+pc3*(b)*(b)+pc4*(b)*(n)+pc5*(n)*(n))-u2*b-u3*b*n),...
sos(0.25*(pc0+pc1*(b+2)+pc2*(n+1)+pc3*(b+2)*(b+2)+pc4*(b+2)*(n+1)+pc5*(n+1)*(n+1))+ (1-0.25)*(pc0+pc1*(b)+pc2*(n+1)+pc3*(b)*(b)+pc4*(b)*(n+1)+pc5*(n+1)*(n+1))-(pc0+pc1*(b)+pc2*(n)+pc3*(b)*(b)+pc4*(b)*(n)+pc5*(n)*(n))+u4*b+u5*b*n),...
u1>=0, u2>=0, u3>=0, u4>=0, u5>=0];
options=sdpsettings('verbose',1,'solver','sedumi');
model = solvesos(FeasibilityConstraints,[ ],options,[pc0;pc1;pc2;pc3;pc4;pc5;v1;v2;v3;v4;v5;v6;u1;u2;u3;u4;u5]); 
 vpc0=value(pc0);  
 vpc1=value(pc1);  
 vpc2=value(pc2);  
 vpc3=value(pc3);  
 vpc4=value(pc4);  
 vpc5=value(pc5); 