function guihua01(f)
n=length(f);
f=f(:);
f=f';
a=zeros(2*n,n*n);
for i=1:n;
a(i,(i-1)*n+1:n*i)=1;
a(n+i,i:n:n*n)=1;
end
b=ones(2*n,1);
[x,fval]=linprog(f,[],[],a,b,zeros(n*n,1),ones(n*n,1));
X=reshape(x,n,n);
opt=X;
disp('最优方案为:'),disp(opt);
disp('结果为:'),disp(fval);