1 c=[3 8 2 10 3;
 2    8 7 2 9 7;
 3    6 4 2 7 5
 4    8 4 2 3 5;
 5    9 10 6 9 10];
 6 c=c(:);%将矩阵C按列拉直，然后赋给C，例如矩阵C=[1,2,3;4,5,6]，操作完后就是列向量1,4,2,5,3,6
 7 a=zeros(10,25);
 8 for i=1:5
 9     a(i,(i-1)*5+1:5*i)=1;
10     %通过循环，将a的第一行的1到5，第二行的6到10一直到第五行的21到25赋值为1
11     
12     a(5+i,i:5:25)=1;
13     %将a的第六行的1、41026、11、16、21
14     %第七行的2、7、12、17、22
15     %第八1653行的3、8、13、18、23
16     %第九行的4、9、14、19、24
17     %第十行的5、10、15、20、25赋值为1
18 end
19 b=ones(10,1);
20 intcon=[1:25];
21 lb=zeros(25,1);
22 ub=ones(25,1)
23 [x,y]=intlinprog(c,intcon,[],[],a,b,lb,ub);
24 x=reshape(x,[5,5]),y