diff --git a/OPF.m b/OPF.m
index cc068d5..3b3c2c2 100644
--- a/OPF.m
+++ b/OPF.m
@@ -56,9 +56,9 @@ mPD=PD0;
 mQD=QD0;
 mVolt=rVolt.*(1+normrnd(0,sigma,length(rVolt),1))';
 bigM=10;
-Vbi=sparse(ones(Busnum,1));
-PDbi=sparse(ones(length(Loadi),1));
-QDbi=sparse(ones(length(Loadi),1));
+Vbi=sparse(0.5*ones(Busnum,1));
+PDbi=sparse(0.5*ones(length(Loadi),1));
+QDbi=sparse(0.5*ones(length(Loadi),1));
 while(abs(Gap)>Precision)
     if KK>kmax
         break;