LMP for PJM-5 bus

% %%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