c=[zeros(1,49),-1,zeros(1,100)];
A=[];
b=[];
Aeq=[4,zeros(1,49),-4,zeros(1,49),zeros(1,50);1,zeros(1,49),zeros(1,50),1,zeros(1,49);zeros(1,48),-1,4,zeros(1,49),-4,zeros(1,50);zeros(1,48),1,4,zeros(1,50),zeros(1,49),4];beq=[1;1;0;4];
vlb=[zeros(1,150)];
vub=[];
[x,fval,ef,out,lag]=linprog(c,A,b,Aeq,beq,vlb,vub)