1
function [p,cap]=max_cap_path(i,j,g)
4
if rhs<>3 then error(39), end
6
if prod(size(i))<>1 then
7
error('First argument must be a scalar')
9
if prod(size(j))<>1 then
10
error('Second argument must be a scalar')
14
// compute lp, la and ls
16
ma=prod(size(g('tail')))
17
if g('directed')==1 then
18
[lp,la,ls]=m6ta2lpd(g('tail'),g('head'),n+1,n)
20
[lp,la,ls]=m6ta2lpu(g('tail'),g('head'),n+1,n,2*ma)
23
if g('edge_max_cap')==[] then
28
// compute max capacity path
29
[c,v]=m6chcm(i,la,lp,ls,n,g('edge_max_cap'))
30
p=m6prevn2p(i,j,v,la,lp,ls,g('directed'))