diff --git a/OPF.m b/OPF.m
index e7321e4..ece53b5 100644
--- a/OPF.m
+++ b/OPF.m
@@ -39,7 +39,7 @@ for badDataNode=1:1
 %         mVolt(badDataNode)=rVolt(badDataNode)*(1-sigma*6);
         %% 画Case A的图
         % figure('Color',[1 1 1]);
-        [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([22],PD0,QD0,mVolt,sigma);%全部有
+        [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([19:22],PD0,QD0,mVolt,sigma);%全部有
         %         badDataResult(I,badDataNode)=sum(Vbi);
         %         badDataLocation(1:33,I)=Vbi;
         %         badDataLocation(34,I)=sum(abs((rVolt-Volt)./rVolt./length(rVolt)))+sum(abs( (UAngel(2:33)-rUAngel(2:33))./rUAngel(2:33)./length(rUAngel(2:33))));