function Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi)

Mat_G=[
    sparse(PD(Loadi));
    sparse(PD(Loadi))-0.001;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
   % PDbi;
    ];
end