%參數a
a=1;
%目標函數系數
f=[0 0 0 0 0 0 0 0 0 0 0 -6 0 0]';
%等式約束矩陣系數
Aeq=[0 0 0 0 0 0 0 0 0 0 0 6 2 1;
      1 1 0 0 0 0 0 0 0 0 0 0 0 0;
      1 0 -1 -1 -1 0 0 0 0 0 0 0 0 0 ;
      0 1 1 -1 0 0 -1 0 0 0 0 0 0 0;
      0 0 0 0 1 -1 1 0 0 0 0 0 0 0 ;
      0 0 0 0 0 1 0 0 1 1 0 0 0 0 ;
      0 0 1 0 0 0 0 -1 0 0 0 0 -1 0;
      0 0 0 0 0 0 -1 1 1 0 0 0 0 0 ;
      0 0 0 0 0 0 0 0 -1 -1 -1 0 0 0;
      0 0 0 1 0 0 0 0 0 0 -1 -1 0 0 ;
      -1 1 0 0 0 0 0 0 0 2 0 0 0 0 ];
beq=[3*a a 0 0 0 0 0 0 0 0 0 ]';
%無不等式約束
A=[];
b=[];
%邊界約束
lb=[0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
ub=[];
%初始點
x0=[0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ]';
%調用函數linprog
[x,fval]=linprog(f,A,b,Aeq,beq,lb,ub,x0)