配流04—全有全无配流算法（2）

说明：指定两点间的客流需求总量，建立费用函数（阻抗函数），一般情况下费用是流量的函数，就可以运用全有全无算法进行配流。

step1:建立费用函数m文件：feiyong.m

function [C,n]=feiyong
syms x1 x2 x3 real;
c1=8+0.1*x1;
c2=10+0.25*x2;
c3=15+0.15*x3;
C=vpa([c1;c2;c3],3);
n=length(C);

step2:主程序m文件：qyqw_feiyong.m

function qyqw_feiyong
syms x1 x2 x3 real
x=[x1;x2;x3];%为了方便使用矩阵进行符号变量替换
clc
disp('========================================================================');
disp('                           《全有全无算法》');
disp('运行环境：MATL