Q=1500;??????????????????????????????????????????????????%OD需求
W=[30?33?25?30?30?35?20?23?20];??????????????????????????%路段零流出行時間
C=[600?500?400?500?500?600?300?700?300];
Mxf=[0?1?0?0?0?0?1?1?0;0?1?1?0?0?1?0?0?0;0?1?1?0?0?0?0?1?1;1?0?0?0?1?0?0?1?1;1?0?0?0?1?1?0?0?0;0?0?0?1?0?0?0?1?1;0?0?0?1?0?1?0?0?0];
syms?lambda?real
for?i=1:numx
????syms?x（i）?real;
end
cont=0;
e=inf;
x=x（1:7）;
X0=zeros（17）;
t=zeros（17）;
tt=t;
while?e>1e-3???????????????????????????%精度判斷
????cont=cont+1;???????????????????????%迭代次數更新
????D（11）=W（11）*[1+0.15*（（x1（14）+x1（15））/C（11））^4];?
D（12）=W（12）*[1+0.15*（（x1（11）+x1（12）+x1（13））/C（12））^4];?
D（13）=W（13）*[1+0.15*（（x1（13）+x1（12））/C（13））^4];?
D（14）=W（13）*[1+0.15*（（x1（16）+x1（17））/C（14））^4];?
D（15）=W（13）*[1+0.15*（（x1（14）+x1（15））/C（15））^4];?
D（16）=W（13）*[1+0.15*（（x1（12）+x1（15）+x1（17））/C（16））^4];?
D（17）=W（13）*[1+0.15*（x1（11）/C（17））^4];?
D（18）=W（13）*[1+0.15*（（x1（11）+x1（13）+x1（14）+x1（16））/C（18））^4];?
D（19）=W（13）*[1+0.15*（（x1（13）+x1（14）+x1（16））/C（19））^4];?
R（11）=D（12）+D（17）+D（18）;
R（12）=D（12）+D（13）+D（16）;
R（13）=D（12）+D（13）+D（19）+D（18）;
R（14）=D（11）+D（15）+D（19）+D（18）;
R（15）=D（11）+D（15）+D（16）;
R（16）=D（14）+D（19）+D（18）;
R（17）=D（14）+D（16）;
????Y1=Mxf（index:）.*Q;????????????????%全有全無法求輔助流量
????S=Y1-X1;???????????????????????????%搜索方向
????X2=X1+lambda*S;????????????????????%先將X2用X1和lambda進行表達
????t=（W）.*（1+0.15*（X2./Cmax）.^4）;?????%含lambd