%?%%case4
GSF=[];
B=[-2050/3010/3;50/3-155/62.520/3;02.5-75/1420/7;10/320/320/7-90/7];
p=[-1?0?0];
theta123=（B（1:31:3）^-1）*p‘;

???GSF（11）=（theta123（1）-theta123（2））/0.06;?
???GSF（21）=0;
???GSF（31）=（theta123（1）-0）/0.3;
???GSF（41）=（theta123（2）-0）/0.15;
???GSF（51）=（theta123（3）-0）/0.35;
???GSF（61）=（theta123（2）-theta123（3））/0.4;
???
p=[0?-1?0];
theta123=（B（1:31:3）^-1）*p‘;

???GSF（12）=（theta123（1）-theta123（2））/0.06;?
???GSF（22）=0;
???GSF（32）=（theta123（1）-0）/0.3;
???GSF（42）=（theta123（2）-0）/0.15;
???GSF（52）=（theta123（3）-0）/0.35;
???GSF（62）=（theta123（2）-theta123（3））/0.4;

???p=[0?0?-1];
theta123=（B（1:31:3）^-1）*p‘;

???GSF（13）=（theta123（1）-theta123（2））/0.06;?
???GSF（23）=0;
???GSF（33）=（theta123（1）-0）/0.3;
???GSF（43）=（theta123（2）-0）/0.15;
???GSF（53）=（theta123（3）-0）/0.35;
???GSF（63）=（theta123（2）-theta123（3））/0.4;
???
???
???f=[10?15?20?30];
lb=[0?0?0?0];
ub=[300?50?60?200];
b=[135.94];
A=[GSF（11）GSF（12）GSF（12）GSF（13）];
Aeq=[1?1?1?1];
beq=[450];
[xfvalexitflagoutputlamda]=linprog（fAbAeqbeqlbub）;
MP=lamda.eqlin;
LMP=-lamda.eqlin-lamda.ineqlin*GSF（11:3）;
LMP=[LMP-lamda.eqlin];
gencost1=10*x（1）+15*x（2）+20*x（3）+30*x（4）;
genrev=LMP（1）*x（1）+LMP（2）*（x（2）+x（3））+LMP（3）*x（4）;


%%
%case1