function [y,phi] = Ksi(K,m,n)
[Ainf,Binf] = lipab(m,n);

[phi,y] = linprog(K,Ainf,Binf,[],[],0*K+1,0*K+1e5);
phi = phi';
return