m=2
n=3

[Ainf,Binf] = lipab(m,n);

global cc

A0=pnormdim(rand(m,m,n-1),1);
C0 = pnormdim(rand(m,1),1);
%C0 = zeros(m,1); C0(2)=1;

phi = optmasklp(A0,C0);
phi = phi-min(phi)

K1 = getK1(cc,m,n);

[Ainf,Binf] = lipab(m,n-1);

[x,fval] = linprog(K1,Ainf,Binf,[],[],[],[]);
phi1 = x';
phi1 = phi1-min(phi1)