diff --git a/checkSSatisfied.m b/checkSSatisfied.m
index 860a9a6..7769ef9 100644
--- a/checkSSatisfied.m
+++ b/checkSSatisfied.m
@@ -19,7 +19,7 @@ output_args=Sp3+pLoadABC;
 ck=output_args;
 ck(3*(Balance-1)+1:3*(Balance-1)+3)=0;
 if any(abs(ck)>1e-5)
-    fprintf('反推回的功率不匹配。潮流方程约束不满足。');
+    fprintf('反推回的功率不匹配。潮流方程约束不满足。\n');
 end
 end
 
diff --git a/dataRead.m b/dataRead.m
index 02994f2..25d5cd3 100644
--- a/dataRead.m
+++ b/dataRead.m
@@ -1,5 +1,5 @@
 function [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
-    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,nodeNum,Balance,phaseABCY] = dataRead(lineZ,dataFile )
+    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY] = dataRead(lineZ,dataFile )
 data=dlmread(dataFile);
 zeroEntries=find(data(:,1)==0);
 lines=data(zeroEntries(1)+1:zeroEntries(2)-1,:);
diff --git a/lineWithConfig.m b/lineWithConfig.m
index 49d78b6..ef63625 100644
--- a/lineWithConfig.m
+++ b/lineWithConfig.m
@@ -15,6 +15,15 @@ phase3Y(3,2)=phase3Y(3,1);
 phase3Y(2,3)=phase3Y(1,3);
 phase3Y(2,2)=phase3Y(1,1);
 phase3Y(3,3)=phase3Y(1,1);
+% phase3Y(1,3)=10;
+% phase3Y(3,2)=10;
+% phase3Y(2,3)=10;
+% phase3Y(2,2)=phase3Y(1,1);
+% phase3Y(3,3)=phase3Y(1,1);
+% phase3Y(1,2)=10;
+% phase3Y(2,1)=10;
+% phase3Y(3,1)=10;
+% phase3Y(1,3)=10;
 %标幺值
 baseHighU=4.16;
 baseLowU=0.48;
diff --git a/run.m b/run.m
index b385c6d..4396098 100644
--- a/run.m
+++ b/run.m
@@ -2,8 +2,7 @@ clc
 clear
 lineZ=readLineZ('.\feeder13\lineParameter.txt');
 [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
-    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data2.txt');
-%fsY1(1,1)=fsY1(1,1)+1j*1e-10;
+    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data1.txt');
 a=exp(1j*2*pi/3);
 Tp2f=1/3*[1 1 1;
     1 a a^2;
@@ -52,7 +51,7 @@ while(k<=kmax && maxD> EPS)
     k=k+1;
     [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
         PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
-    maxD=max([dP;dQ;]);
+    maxD=max(abs([dP;dQ;]));
     jaco=Jacobi(Balance,busNum,QGi,Vmf1,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵
     [dV, dVangle]=Solv(busNum,jaco,dP,dQ);%解出修正量
     [Vmf1, Vaf1]=Modify(Vmf1,Vaf1,dV,dVangle,1);
@@ -76,6 +75,15 @@ while(k<=kmax && maxD> EPS)
     Vf2=fsY2\If2;
 end
 Vf1=Vmf1.*exp(1j*Vaf1);
+%%
+% V=[1.01*exp(1j*0);1.01*exp(-1j*pi/360*120);1.01*exp(-1j*pi/360*240)];
+% abs(Tp2f*V)
+% Tf2p*Tp2f*V
+% Vf0=Vf0+refV0;
+% Vf2=Vf2+refV2;
+% Vf0=Vf0*exp(-1j*(-1.0471975511966));
+% Vf1=Vf1*exp(-1j*(1.0471975511966));
+%%
 (Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3;
 conj(Tf2p*[If0(2);If1(2);If2(2)]).*(Tf2p*[Vf0(2);Vf1(2);Vf2(2)]);
 IpABC=Tf2p*conj([If0';If1';If2']);
@@ -83,6 +91,8 @@ IpABC=Tf2p*conj([If0';If1';If2']);
 VoltpABC=Tf2p*conj([ Vf0'; Vf1'; Vf2']);
 disp(['                  A                        B                          C'])
 full(abs(VoltpABC'))
+fprintf('节点号对应\n');
+disp([setIJ,nodeNum ])
 checkSSatisfied(Balance,phaseABCY,VoltpABC, ...
     phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP, ...
     phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ );