fileID = fopen('prob6.lp','w');
fprintf(fileID,'min\n');
fprintf(fileID,'obj: x2');
for i=3:contf
    fprintf(fileID,' + x%d', i);
end
for m= 2:n-1
    for i=2:contf
        fprintf(fileID,' + x%d', i+(m-1)*contf);
    end
end
fprintf(fileID,'\n\n');
fprintf(fileID,'subject to\n\n');
for m=1:n-1
    fprintf(fileID,'c%d: x%d = 1\n\n', m,(m-1)*contf+1);
end
for m=1:n-2
    for i=1:8*(l-3)+4*(l-2-k)
        fprintf(fileID,'c%d: %dx%d', i+n-1+(m-1)*(8*(l-3)+4*(l-2-k)), final((m-1)*(8*(l-3)+4*(l-2-k))+i,1,1),1+(m-1)*contf);
        for j=2:contf
            if (final((m-1)*(8*(l-3)+4*(l-2-k))+i,j,1)==1)
                fprintf(fileID,' + x%d', j+(m-1)*contf);
            end 
        end
        fprintf(fileID,' >= 1\n\n');
    end
end
m=n-1;
for i=1:4*(l-1)+(l-4)*(l-4)+4*(l-2-k)
    fprintf(fileID,'c%d: %dx%d', i+n-1+(m-1)*(8*(l-3)+4*(l-2-k)), final((m-1)*(8*(l-3)+4*(l-2-k))+i,1,1),1+(m-1)*contf);
    for j=2:contf
        if (final((m-1)*(8*(l-3)+4*(l-2-k))+i,j,1)==1)
            fprintf(fileID,' + x%d', j+(m-1)*contf);
        end 
    end
    fprintf(fileID,' >= 1\n\n');
end
fprintf(fileID,'bounds\n');
fprintf(fileID,'binary\n');
for i=1:contf*(n-1)
    fprintf(fileID,'x%d\n', i);
end
fclose(fileID);