diff --git a/feeder13/data1.txt b/feeder13/data1.txt
new file mode 100644
index 0000000..3a1a90a
--- /dev/null
+++ b/feeder13/data1.txt
@@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 70 90 70
+631 60 30 60 30 70 30
+632 60 30 70 30 60 30
+0
+
+
diff --git a/run.m b/run.m
index 9a5ef51..d746994 100644
--- a/run.m
+++ b/run.m
@@ -269,30 +269,30 @@ for II=1:3
     
     if II==1
         fprintf('ÕýÐò\n');
-        wV1r(1:2)=1e10;
-        wV1i(1:2)=1e10;
-        wI1r(1:2)=1e10;
-        wI1i(1:2)=1e10;
+%         wV1r(1:2)=1e10;
+%         wV1i(1:2)=1e10;
+%         wI1r(1:2)=1e10;
+%         wI1i(1:2)=1e10;
         tic
         [ V1r,V1i,I1r,I1i,isConverged1 ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 );
         toc
     end
     if II==2
         fprintf('¸ºÐò\n');
-        wV2r(1:2)=1e10;
-        wV2i(1:2)=1e10;
-        wI2r(1:2)=1e10;
-        wI2i(1:2)=1e10;
+%         wV2r(1:2)=1e10;
+%         wV2i(1:2)=1e10;
+%         wI2r(1:2)=1e10;
+%         wI2i(1:2)=1e10;
         tic
         [ V2r,V2i,I2r,I2i,isConverged2 ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 );
         toc
     end
     if II==3
         fprintf('ÁãÐò\n');
-        wV0r(1:2)=1e10;
-        wV0i(1:2)=1e10;
-        wI0r(1:2)=1e10;
-        wI0i(1:2)=1e10;
+%         wV0r(1:2)=1e10;
+%         wV0i(1:2)=1e10;
+%         wI0r(1:2)=1e10;
+%         wI0i(1:2)=1e10;
         tic
         [ V0r,V0i,I0r,I0i,isConverged0 ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
         toc