diff --git a/Run_YALMIP2.m b/Run_YALMIP2.m
index 7486365..d5b4310 100644
--- a/Run_YALMIP2.m
+++ b/Run_YALMIP2.m
@@ -112,8 +112,11 @@ for I=1:1
     SEVolt=x(1:length(Volt));
     SEVAngel=x(length(Volt)+1:end);
     fprintf('×î´óÆ«²î\n')
-    PD=x(1:length(PDi));
-    QD=x(length(PDi)+1:length(PDi)+length(QDi));
+    AngleIJ=sparse(r,c,SEVAngel(r)-SEVAngel(c)-Angle,Busnum,Busnum);
+    zeroP=diag(SEVolt)*Y.*cos(AngleIJ)*SEVolt;
+    zeroQ=diag(SEVolt)*Y.*sin(AngleIJ)*SEVolt;
+    PD=-zeroP(PDi);
+    QD=-zeroQ(QDi);
     maxDev1=MaxDeviation(sigma,mVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD)
     fprintf('ͳ¼ÆÆ«²î\n')
     statDev1=StatDeviation(sigma,rVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD)