diff --git a/OPF.m b/OPF.m
index df33f88..d448ad4 100644
--- a/OPF.m
+++ b/OPF.m
@@ -216,7 +216,7 @@ mVoltABCV=rVoltABCV.*VoltSigma;
 PD3PSigma=(1+normrnd(0,sigma,length(rPD3P),1));
 mPD3P=rPD3P.*PD3PSigma;
 QD3PSigma=(1+normrnd(0,sigma,length(rQD3P),1));
-mQD3P=rPD3P.*QD3PSigma;
+mQD3P=rQD3P.*QD3PSigma;
 %Á¿²â·½²î
 wVolt=1./(abs(mVoltABCV*sigma).^2);
 wPD=1./(abs(mPD3P*sigma).^2);