variational inequality problem

function flow_choice（）

%   global variation
global ts;
global con;
%*********************************************************
%   initial map : 5 nodes 5 links destination is node 5!
%   init OD with origin demand: 5 nodes and 30 time span
od=zeros（530）;
od（11）=14; od（12）=9;              od（14）=13;
od（21）=14; od（22）=10;             od（24）=12;
od（41）=5;  od（42）=4;  od（43）=2;  od（44）=1;

%    initial map connection
con=[0 1 0 0 0
    0 0 2 3 0
    0 0 0 4 0
    0 0 0 0 5
    0 0 0 0 0];

%   init time-space: 5 links and 30 time span each cell has struct:{influx flow cost}
%   init link state structure
linkstate=struct（‘influx‘0‘flow‘0‘cost‘1）;
ts=cell（530）;
for i_1=1:5
    for j_1=1:30
        ts{i_1j_1}=linkstate;
    end
end

%Test
%sp =