diff --git a/OPF.m b/OPF.m
index a42e3f1..367410a 100644
--- a/OPF.m
+++ b/OPF.m
@@ -75,27 +75,6 @@ while(k<=kmax && maxD> EPS)
     iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
     iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
     iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
-%    %全部转换为负荷电流
-%    CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
-%    CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
-%    CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
-%    %转换为序电流
-%    f012=Tp2f*conj([CurpA';CurpB';CurpC']);
-%    %把三序电流分离出来
-%    If0=conj(f012(1,:)');
-%    If1=conj(f012(2,:)');
-%    If2=conj(f012(3,:)');
-%    %试着算一下正序电流
-%     fsY11*V1;
-%    %形成负荷序电流的测量值
-%    mIf0=If0;
-%    mIf1=If1;
-%    mIf1(3)=-mIf1(2);
-%    mIf2=If2;
-%    %计算
-%     fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus
-%     Num);%这里要置0，置1，否则是奇异的
-    %%做最小二乘法
     [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
         PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
     maxD=max(abs([dP;dQ;]));
@@ -217,7 +196,6 @@ while(abs(Gap)>Precision)
     if KK>kmax
         break;
     end
-%     plotGap(KK+1)=Gap;
     Init_u=Gap/2/RestraintCount*CenterA;
     AngleIJMat=0;
     %% 开始计算OPF