clc,clear
a=struct2array(load('22.mat'))
a(:,1)=[]
a=a(:)
b=ceil(a0.2)
xy1=zeros(104,1)
xy2=zeros(104,1)
z=zeros(104,1)
by1=zeros(104,1)
by2=zeros(104,1)
for i=1:104
c=a(i)
sl=ceil(c/10)
by1(i)=@(xy1,a,b,sl)50+sum(xy1(1:i-1,:),1)-4a(i)-4sum(b(1:i-1,:),1)-sl
by2(i)=@(xy2)13+sum(xy2(1:i-1,:),1)-a(i)-sum(b(1:i-1,:),1)
xl=xy1(i)+sl
z(i)=@(xy1,xy2)100xy1(i)+200xy2(i)+by15+by210+xl10
end
by1f=matlabFunction(by1)
prob=optimproblem
x=optimvar('xy',2,'LowerBound',4);
prob.Objective=sum(z)
prob.Constraints.con1=13+sum(xy2(1:i,:),1)-sum(b(1:i-1,:),1)>=a(i+1)
prob.Constraints.con2=50+sum(xy1(1:i,:),1)-4sum(b(1:i-1,:),1)-4a(i)>=4a(i+1)
prob.Constraints.con3=52+4sum(xy2(1:i,:),1)<=sum(xy1(1:i,:),1)+50
[sol,fval,flag,out]=solve(prob),sol.x