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

Mat_G=[
    sparse(PD(Loadi));
    sparse(QD(Loadi));
    Volt'-bigM*Vbi-mVolt'-sigma*RealVolt';
    Volt'+bigM*Vbi+mVolt'+sigma*RealVolt';
    Vbi;
    ];
end