#C27a.txt: May 2, 2013, Traffic flow (under construction)

Help:=proc(): print(` PotF(J,R,DelayJ,y,x) `): end:

JR1:=Optimization[Minimize]:


G1:=[{2},{1,3},{1,2}];
J1:=[[1,2],[2,1],[2,3],[3,2],[3,1]];
S1:=[[1,2],[1,3],[2,1],[2,3],[3,1],[3,2]];
R1:=[[1,2],[1,2,3],[2,1],[2,3],[3,1],[3,2,1],[3,2],[3,1,2]];
DelayJ1:=[10*y+4,35,12*y+15,40*y+11,50*y+1]:

#IsEdgeRoute(Edge,Route): Is Edge part of Route
IsEdgeInRoute:=proc(Edge,Route) local i:

member(Edge, {seq([Route[i],Route[i+1]],i=1..nops(Route)-1)}):


end:


#PotF(J,R,DelayJ,y,x): the Traffic potential function with
#edges J1 and Routes R1
PotF:=proc(J,R,DelayJ,y, x) local W,j,z:
#W is for Wait (Wisecrack due to Pat Devlin)
W:=add( int(DelayJ[j], y=0..z[j])    , j=1..nops(J)):

end: