MATLAB利用yamlip求解优化问题

clear all
close all
Cstart = [5000 1600 2400 1200]‘;
Len = [6 3 3 2 4 4 2]‘;
Cmin = [2250 1800 3750 4800]‘;
Cadd = [2.7 2.2 1.8 3.8]‘;
Dem = [12000 32000 25000 36000 25000 30000 18000]‘;
Pmax = [1750 1500 2000 3500]‘;
Pmin = [750 1000 1200 1800]‘;
Avail = [10 4 8 3]‘;
start = intvar（47）;
work = intvar（47）;
padd = intvar（47）;
f = sum（Cstart‘ * start） + Len‘ * （Cmin‘ * work + Cadd‘ * padd）‘;
F = [];
F = [F start>=0 padd>=0];
for i = 1:4
    F = [F 0<=work（i:）<=Avail（i）];
end
for i = 1:4
    for j = 1:7
        F = [F padd（ij