apis/examples_matlab_control/transport13_8m_source.html

function transport13(varargin)


    % check workspace info from arguments

    if nargin > 0

        wsInfo = gams.control.WorkspaceInfo();

        wsInfo.systemDirectory = varargin{1};

        ws = gams.control.Workspace(wsInfo);

    else

        ws = gams.control.Workspace();

    end


    % prepare input data

    plants = {'Seattle', 'San-Diego'};

    markets = {'New-York', 'Chicago', 'Topeka'};

    capacity = containers.Map();

    capacity('Seattle') = 350;

    capacity('San-Diego') = 600;

    demand = containers.Map();

    demand('New-York') = 325;

    demand('Chicago') = 300;

    demand('Topeka') = 275;

    distance = containers.Map();

    distance('Seattle.New-York') = 2.5;

    distance('Seattle.Chicago') = 1.7;

    distance('Seattle.Topeka') = 1.8;

    distance('San-Diego.New-York') = 2.5;

    distance('San-Diego.Chicago') = 1.8;

    distance('San-Diego.Topeka') = 1.4;


    t = TransportModel(ws);


    for p = plants

        t.i.addRecord(p{1});

    end


    for m = markets

        t.j.addRecord(m{1});

    end


    for p = plants

        rec = t.a.addRecord(p{1});

        rec.value = capacity(p{1});

    end


    for m = markets

        rec = t.b.addRecord(m{1});

        rec.value = demand(m{1});

    end


    for p = plants

        for m = markets

            rec = t.d.addRecord(p{1}, m{1});

            rec.value = distance([p{1}, '.', m{1}]);

        end

    end


    rec = t.f.addRecord();

    rec.value = 90;


    t.opt.setAllModelTypes('cplex');

    t.run([], gams.control.PrintStream());


    fprintf('Objective: %g\n', t.z.record.level);


    for x = t.x.records

        fprintf('x(%s,%s): level=%g marginal=%g\n', x{1}.keys{:}, x{1}.level, x{1}.marginal);

    end

end