clear;
 clc;
sdpvar x n
sdpvar pc0 pc1 pc2 pc3 pc4 pc5
sdpvar v1 v2 v3 v4 v5 v6
sdpvar u1
FeasibilityConstraints=[sos(pc0+pc1*(x)+pc2*(n)+pc3*(x)*(x)+pc4*(x)*(n)+pc5*(n)*(n)-(0.25*n^2+0.25*n)+v1*(x-0)+v2*(x-0)*x+v3*(x-0)*n),...
sos(x-(pc0+pc1*(x)+pc2*(n)+pc3*(x)*(x)+pc4*(x)*(n)+pc5*(n)*(n))+v4*(0-n)+v5*(0-n)*x+v6*(0-n)*n),...
sos(0.5*(pc0+pc1*(x+n)+pc2*(n-1)+pc3*(x+n)*(x+n)+pc4*(x+n)*(n-1)+pc5*(n-1)*(n-1))+ (1-0.5)*(pc0+pc1*(x)+pc2*(n-1)+pc3*(x)*(x)+pc4*(x)*(n-1)+pc5*(n-1)*(n-1))-(pc0+pc1*(x)+pc2*(n)+pc3*(x)*(x)+pc4*(x)*(n)+pc5*(n)*(n))-u1*(n-0)),...
u1>=0];
options=sdpsettings('verbose',1,'solver','sedumi');
model = solvesos(FeasibilityConstraints,[ ],options,[pc0;pc1;pc2;pc3;pc4;pc5;v1;v2;v3;v4;v5;v6;u1]); 
 vpc0=value(pc0);  
 vpc1=value(pc1);  
 vpc2=value(pc2);  
 vpc3=value(pc3);  
 vpc4=value(pc4);  
 vpc5=value(pc5); 
cof_list = [vpc0,vpc1,vpc2,vpc3,vpc4,vpc5];
ref_list=(cof_list>10^-6).*cof_list;
L=ref_list(1)*1+ref_list(2)*(x)+ref_list(3)*(n)+ref_list(4)*(x)*(x)+ref_list(5)*(x)*(n)+ref_list(6)*(n)*(n)