apis/examples_matlab_control/transport2_8m_source.html

function transport2(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


    data = {

    'Sets                                                   '

    '  i   canning plants   / seattle, san-diego /          '

    '  j   markets          / new-york, chicago, topeka / ; '

    'Parameters                                             '

    '                                                       '

    '  a(i)  capacity of plant i in cases                   '

    '                     /    seattle     350              '

    '                          san-diego   600  /           '

    '                                                       '

    '  b(j)  demand at market j in cases                    '

    '                     /    new-york    325              '

    '                          chicago     300              '

    '                          topeka      275  / ;         '

    '                                                       '

    'Table d(i,j)  distance in thousands of miles           '

    '               new-york       chicago      topeka      '

    '  seattle          2.5           1.7          1.8      '

    '  san-diego        2.5           1.8          1.4  ;   '

    '                                                       '

    'Scalar f  freight in dollars per case per thousand miles  /90/'};

    data = sprintf('%s\n', data{:});


    model = {

    'Sets                                                                    '

    '      i   canning plants                                                '

    '      j   markets                                                       '

    '                                                                        '

    'Parameters                                                              '

    '      a(i)   capacity of plant i in cases                               '

    '      b(j)   demand at market j in cases                                '

    '      d(i,j) distance in thousands of miles                             '

    'Scalar f  freight in dollars per case per thousand miles;               '

    '                                                                        '

    '$if not set incname $abort ''no include file name for data file provided'''

    '$include %incname%                                                      '

    '                                                                        '

    ' Parameter c(i,j)  transport cost in thousands of dollars per case ;    '

    '                                                                        '

    '            c(i,j) = f * d(i,j) / 1000 ;                                '

    '                                                                        '

    ' Variables                                                              '

    '       x(i,j)  shipment quantities in cases                             '

    '       z       total transportation costs in thousands of dollars ;     '

    '                                                                        '

    ' Positive Variable x ;                                                  '

    '                                                                        '

    ' Equations                                                              '

    '                                                                        '

    '      cost        define objective function                             '

    '      supply(i)   observe supply limit at plant i                       '

    '       demand(j)   satisfy demand at market j ;                         '

    '                                                                        '

    '  cost ..        z  =e=  sum((i,j), c(i,j)*x(i,j)) ;                    '

    '                                                                        '

    '  supply(i) ..   sum(j, x(i,j))  =l=  a(i) ;                            '

    '                                                                        '

    '  demand(j) ..   sum(i, x(i,j))  =g=  b(j) ;                            '

    '                                                                        '

    ' Model transport /all/ ;                                                '

    '                                                                        '

    ' Solve transport using lp minimizing z ;                                '

    '                                                                        '

    'Display x.l, x.m ;                                                      '

    '                                                                        '};

    model = sprintf('%s\n', model{:});


    % write 'data' into file 'tdata.gms' under Workspace's working directory

    fid = fopen(fullfile(ws.workingDirectory, 'tdata.gms'), 'w');

    fprintf(fid, data);

    fclose(fid);


    % create Job 't2' from the 'model' string variable

    t2 = ws.addJobFromString(model);

    % create Option 'opt' and define 'incname' as 'tdata'

    opt = ws.addOptions();

    opt.defines('incname', 'tdata');

    % run Job 't2' with Options 'opt'

    t2.run(opt);


    % retrieve Variable 'x' from Job's output databases

    for x = t2.outDB.getVariable('x').records

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

    end


    % clear option and database

    opt.dispose();

    t2.outDB.dispose();


    % remove working directory

    rmdir(ws.workingDirectory, 's');


end